FITS driver improved: new functions for Scale Offset and FITS World C…

…oordinate System read and write (#1298)

.gitignore

@@ -13,3 +13,4 @@
 .cproject
 .idea
 .vscode
+

autotest/gdrivers/fits.py

@@ -38,7 +38,7 @@
 
 @pytest.mark.parametrize(
     'filename',
-    ['byte', 'int16', 'int32', 'float32', 'float64']
+    ['byte', 'int16', 'uint16', 'int32', 'uint32', 'float32', 'float64']
 )
 def test_fits(filename):
     driver = gdal.GetDriverByName('FITS')
@@ -54,8 +54,7 @@ def test_fits(filename):
     ds2 = None
     driver.Delete('tmp/' + filename + '.fits')
 
-
-def fits_metadata():
+def test_fits_metadata():
     driver = gdal.GetDriverByName('FITS')
 
     ds = gdal.Open('../gcore/data/byte.tif')
@@ -82,3 +81,36 @@ def fits_metadata():
     ds2 = None
 
     assert md['TEST2'] == 'test_value2'
+
+def test_fits_nodata():
+    driver = gdal.GetDriverByName('FITS')
+
+    ds = gdal.Open('../gcore/data/nodata_byte.tif')
+    ds2 = driver.CreateCopy('tmp/nodata_byte.fits', ds)
+    ds2 = None
+    gdal.Unlink('tmp/nodata_byte.fits.aux.xml')
+
+    ds2 = gdal.Open('tmp/nodata_byte.fits')
+    nd = ds2.GetRasterBand(1).GetNoDataValue()
+    ds2 = None
+    driver.Delete('tmp/nodata_byte.fits')
+
+    assert nd == 0
+
+def test_fits_offscale():
+    driver = gdal.GetDriverByName('FITS')
+
+    ds = gdal.Open('../gdrivers/data/offscale_byte.tif')
+    ds2 = driver.CreateCopy('tmp/offscale_byte.fits', ds)
+    ds2 = None
+    gdal.Unlink('tmp/offscale_byte.fits.aux.xml')
+
+    ds2 = gdal.Open('tmp/offscale_byte.fits')
+    offset = ds2.GetRasterBand(1).GetOffset()
+    scale = ds2.GetRasterBand(1).GetScale()
+    ds2 = None
+    driver.Delete('tmp/offscale_byte.fits')
+
+    assert offset == -0.0039525691699605
+    assert scale == 1.00395256917
+