I am doing a real-time SAR spot simulation where the SAR orbits a
target point and saves a series of the georeferenced spot image
files. My initial implementation generated non-GDAL GeoTIFF or JPEG
files. I use a least-squares algorithm to develop the 6 affine
geotransform values. This all works fine, and the files display
nicely in my GIS viewer. I have since brought in GDAL 2.2.3 to
create NITF files. Unfortunately, over 30% of the NITF files lose the
geotransform info. Rather, they show generic GCPs. The “good” NITF
files display OK in my GIS viewer. As a sanity test, I took a series
of non-GDAL GeoTIFF files and ran them through gdal_translate to
NITF. Interestingly, around 30% of the NITF files again dropped the
geotransform data. I’ve attached below gdalinfo reports of a
gdal_translate input GeoTIFF file and the output NITF file for a
failing example. My geotransform data appears to be well-behaved.
Question: What could cause geotransform data to be rejected? Thanks…..
Non-GDAL gdal_translate input GeoTIFF file:
$ gdalinfo sar2023071722525322.tif
Driver: GTiff/GeoTIFF
Files: sar2023071722525322.tif
Size is 1024, 1024
Coordinate System is:
PROJCS["VITAL SAR Simulation",
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4326"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",-81],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]]]
GeoTransform =
586659.1336352756, -0.2027675805429681, -0.2274441341933198
2969181.096499263, -0.2274441341933198, 0.2027675805429681
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left ( 586659.134, 2969181.096) ( 80d 7'40.06"W, 26d50'29.50"N)
Lower Left ( 586426.231, 2969388.731) ( 80d 7'48.44"W, 26d50'36.30"N)
Upper Right ( 586451.500, 2968948.194) ( 80d 7'47.64"W, 26d50'21.98"N)
Lower Right ( 586218.597, 2969155.828) ( 80d 7'56.02"W, 26d50'28.78"N)
Center ( 586438.865, 2969168.462) ( 80d 7'48.04"W, 26d50'29.14"N)
Band 1 Block=1024x8 Type=Byte, ColorInterp=Gray
gdal_translate output NITF file:
$ gdalinfo sar2023071722525322.nitf
Driver: NITF/National Imagery Transmission Format
Files: sar2023071722525322.nitf
Size is 1024, 1024
Coordinate System is `'
GCP Projection =
PROJCS["UTM Zone 17, Northern Hemisphere",
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0,
AUTHORITY["EPSG","8901"]],
UNIT["degree",0.0174532925199433,
AUTHORITY["EPSG","9122"]],
AUTHORITY["EPSG","4326"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",-81],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",0],
UNIT["Meter",1]]
GCP[ 0]: Id=UpperLeft, Info=
(0.5,0.5) -> (586659,2969181,0)
GCP[ 1]: Id=UpperRight, Info=
(1023.5,0.5) -> (586451,2968948,0)
GCP[ 2]: Id=LowerRight, Info=
(1023.5,1023.5) -> (586219,2969156,0)
GCP[ 3]: Id=LowerLeft, Info=
(0.5,1023.5) -> (586426,2969389,0)
Metadata:
NITF_ABPP=08
NITF_CCS_COLUMN=0
NITF_CCS_ROW=0
NITF_CLEVEL=03
NITF_ENCRYP=0
NITF_FBKGC= 0, 0, 0
NITF_FDT=20021216151629
NITF_FHDR=NITF02.10
NITF_FSCATP=
NITF_FSCAUT=
NITF_FSCLAS=U
NITF_FSCLSY=
NITF_FSCLTX=
NITF_FSCODE=
NITF_FSCOP=00000
NITF_FSCPYS=00000
NITF_FSCRSN=
NITF_FSCTLH=
NITF_FSCTLN=
NITF_FSDCDT=
NITF_FSDCTP=
NITF_FSDCXM=
NITF_FSDG=
NITF_FSDGDT=
NITF_FSREL=
NITF_FSSRDT=
NITF_FTITLE=
NITF_IALVL=0
NITF_IC=NC
NITF_ICAT=VIS
NITF_ICORDS=N
NITF_IDATIM=20021216151629
NITF_IDLVL=1
NITF_IGEOLO=175866592969181175864512968948175862192969156175864262969389
NITF_IID1=Missing
NITF_IID2=
NITF_ILOC_COLUMN=0
NITF_ILOC_ROW=0
NITF_IMAG=1.0
NITF_IMODE=B
NITF_IREP=MONO
NITF_ISCATP=
NITF_ISCAUT=
NITF_ISCLAS=U
NITF_ISCLSY=
NITF_ISCLTX=
NITF_ISCODE=
NITF_ISCRSN=
NITF_ISCTLH=
NITF_ISCTLN=
NITF_ISDCDT=
NITF_ISDCTP=
NITF_ISDCXM=
NITF_ISDG=
NITF_ISDGDT=
NITF_ISORCE=Unknown
NITF_ISREL=
NITF_ISSRDT=
NITF_ONAME=
NITF_OPHONE=
NITF_OSTAID=GDAL
NITF_PJUST=R
NITF_PVTYPE=INT
NITF_STYPE=BF01
NITF_TGTID=
Corner Coordinates:
Upper Left ( 0.0, 0.0)
Lower Left ( 0.0, 1024.0)
Upper Right ( 1024.0, 0.0)
Lower Right ( 1024.0, 1024.0)
Center ( 512.0, 512.0)
Band 1 Block=1024x1 Type=Byte, ColorInterp=Gray
_______________________________________________
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev
