[gdal-dev] RADARSAT2 NITF RPC georeferencing
Gong, Shawn (Contractor)
Shawn.Gong at drdc-rddc.gc.ca
Thu Jun 25 16:04:58 EDT 2009
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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/gdal-dev/attachments/20090625/61cdd2b9/attachment-0001.html
More information about the gdal-dev
mailing list