hi list, RSAT2's NITF product comes in "imagery_##.ntf" and "product.xml" gdalinfo output for these two are listed below.
Georeferencing is off between the two files (their corners don't match using OpenEV). Q1: Is it true that imagery_##.ntf uses affine transform (see GeoTransform tag), and product.xml uses 143 GCPs? Q2: I see that NITF RPC georeferencing has been implemented, http://trac.osgeo.org/gdal/wiki/rfc22_rpc How do I use the RPC in RADARSAT2's NITF file? Will the RPC provide better georeferencing? thanks, Shawn ========================================= imagery_##.ntf Driver: NITF/National Imagery Transmission Format Files: imagery_HH.ntf Size is 8792, 10070 Coordinate System is: GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS 84",6378137,298.257223563, AUTHORITY["EPSG","7030"]], TOWGS84[0,0,0,0,0,0,0], AUTHORITY["EPSG","6326"]], PRIMEM["Greenwich",0, AUTHORITY["EPSG","8901"]], UNIT["degree",0.0174532925199433, AUTHORITY["EPSG","9108"]], AUTHORITY["EPSG","4326"]] GeoTransform = 50.03718817626492, 5.656037108659124e-005, 1.153154346122985e-005 13.76523037807838, 1.126467725830654e-005, -5.535416735635244e-005 Metadata: NITF_FHDR=NITF02.10 NITF_CLEVEL=06 NITF_STYPE=BF01 NITF_OSTAID=RSI NITF_FDT=20081019203042 NITF_FTITLE=RS2_RECTIFIED:2008-10-19T14:51:25.306614Z ... ... RPC Metadata: LINE_OFF=5034 LINE_SCALE=5036 SAMP_OFF=4388 SAMP_SCALE=4538 LONG_OFF=50.3435 LONG_SCALE=0.3068 LAT_OFF=13.5359 LAT_SCALE=0.3283 HEIGHT_OFF=0 HEIGHT_SCALE=700 LINE_NUM_COEFF=0.0006728793 0.2105967 -1.130644 -1.223845e-006 -0.0003698973 -1.002738e-007 -3.118422e-007 -0.0006155271 -0.000504439 0 -1.174964e-008 2.143919e-007 -4.856652e-006 7.659651e-008 -6.164503e-006 -8.210502000000001e-007 -4.267789e-007 -4.321238e-008 0 -2.272615e-008 LINE_DEN_COEFF=1 -0.0004301605 0.0002676128 -3.302012e-007 -3.970334e-006 0 1.027667e-008 -6.566892e-006 2.541771e-006 3.784513e-007 0 0 -1.749914e-006 0 3.37938e-007 3.157953e-006 0 0 0 0 SAMP_NUM_COEFF=0.0003633695 1.147341 0.2557103 -0.03114342 -0.0006508754 0.0008477953 0.0001658026 -5.388591e-005 -0.0002000257 2.006948e-005 -1.95455e-005 -2.907722e-006 -9.035467999999999e-006 -1.276494e-006 1.372415e-006 -1.699765e-006 -5.031646e-007 -6.214799e-005 -1.609981e-006 3.722058e-008 SAMP_DEN_COEFF=1 -0.0001622142 0.0007102185 -0.001001374 -8.652518e-007 4.728022e-005 1.732678e-005 -2.196369e-006 8.936224e-006 -1.514486e-006 -2.812931e-006 0 1.072109e-007 1.4604e-007 -1.900576e-008 -2.564915e-007 0 -6.090468e-006 -3.221172e-007 0 MIN_LONG=50.1901 MAX_LONG=50.4969 MIN_LAT=13.37175 MAX_LAT=13.70005 Corner Coordinates: Upper Left ( 50.0371882, 13.7652304) ( 50d 2'13.88"E, 13d45'54.83"N) Lower Left ( 50.1533108, 13.2078139) ( 50d 9'11.92"E, 13d12'28.13"N) Upper Right ( 50.5344670, 13.8642694) ( 50d32'4.08"E, 13d51'51.37"N) Lower Right ( 50.6505896, 13.3068530) ( 50d39'2.12"E, 13d18'24.67"N) Center ( 50.3438889, 13.5360417) ( 50d20'38.00"E, 13d32'9.75"N) Band 1 Block=1024x1024 Type=UInt16, ColorInterp=Gray ========================================= product.xml Driver: RS2/RadarSat 2 XML Product Files: product.xml lutBeta.xml lutSigma.xml lutGamma.xml imagery_HH.ntf imagery_HV.ntf Size is 8792, 10070 Coordinate System is `' GCP Projection = GEOGCS[,DATUM[,SPHEROID["WGS 1984",6378137,298.2572235629972]],PRIMEM["Greenwich",0],UNIT["degree",0. 0174532925199433]] GCP[ 0]: Id=1, Info= (0,0) -> (50.0368241447173,13.7652233850433,0) GCP[ 1]: Id=2, Info= (879.099976,0) -> (50.0866031135979,13.775165475226,0) GCP[ 2]: Id=3, Info= (1758.19995,0) -> (50.1363862930026,13.7850974953681,0) GCP[ 3]: Id=4, Info= (2637.30005,0) -> (50.1861736889333,13.7950194384848,0) GCP[ 4]: Id=5, Info= (3516.3999,0) -> (50.2359652797982,13.8049312921043,0) GCP[ 5]: Id=6, Info= (4395.5,0) -> (50.2857610923955,13.8148330533982,0) ... ... GCP[139]: Id=140, Info= (6153.7002,10069) -> (50.5012171146149,13.2772541435474,0) GCP[140]: Id=141, Info= (7032.7998,10069) -> (50.5509029539819,13.2871595011868,0) GCP[141]: Id=142, Info= (7911.8999,10069) -> (50.6005928456264,13.297055157266,0) GCP[142]: Id=143, Info= (8791,10069) -> (50.6502867615454,13.3069410980946,0) Metadata: PRODUCT_TYPE=SGF PIXEL_SPACING=6.25000000e+00 LINE_SPACING=6.25000000e+00 BETA_NOUGHT_LUT=lutBeta.xml SIGMA_NOUGHT_LUT=lutSigma.xml GAMMA_LUT=lutGamma.xml SATELLITE_IDENTIFIER=RADARSAT-2 SENSOR_IDENTIFIER=SAR BEAM_MODE=F1 ACQUISITION_START_TIME=2008-10-19T14:51:24.966569Z NEAR_RANGE_INCIDENCE_ANGLE=3.68685226e+01 FAR_RANGE_INCIDENCE_ANGLE=3.98800888e+01 SLANT_RANGE_NEAR_EDGE=9.631139343954010e+05 Subdatasets: SUBDATASET_3_NAME=RADARSAT_2_CALIB:BETA0:product.xml SUBDATASET_3_DESC=Beta Nought calibrated SUBDATASET_2_NAME=RADARSAT_2_CALIB:SIGMA0:product.xml SUBDATASET_2_DESC=Sigma Nought calibrated SUBDATASET_4_NAME=RADARSAT_2_CALIB:GAMMA:product.xml SUBDATASET_4_DESC=Gamma calibrated SUBDATASET_1_NAME=RADARSAT_2_CALIB:UNCALIB:product.xml SUBDATASET_1_DESC=Uncalibrated digital numbers Corner Coordinates: Upper Left ( 0.0, 0.0) Lower Left ( 0.0,10070.0) Upper Right ( 8792.0, 0.0) Lower Right ( 8792.0,10070.0) Center ( 4396.0, 5035.0) Band 1 Block=1024x1024 Type=UInt16, ColorInterp=Undefined Metadata: POLARIMETRIC_INTERP=HH Band 2 Block=1024x1024 Type=UInt16, ColorInterp=Undefined Metadata: POLARIMETRIC_INTERP=HV
_______________________________________________ gdal-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/gdal-dev
