I have multiple raster images without geoinformation. The geoinformation is stored seperately in two other rasters: one contains the Lat-coordinate of every pixel, the another the Long-coordinates. So my question is: What's the best way to bring the information together to update the raster images with geoinformation?
I'm working with IDL 8.4/ENVI5.2 and format of the rasters is binary (.dat). But I could save it in any format.
What did I try so far:I wrote lists of GCPs based on the Lat-and Long-rasters (I only used every 100th pixel and not every pixel because of the size of the image). A list looks like this:
; ENVI Image to Map GCP File; projection info = {Geographic Lat/Lon, WGS-84, units=Degrees}; warp file: test.dat; Map (x,y), Image (x,y);13.239339 53.374766 1 10013.238695 53.372731 1 20013.238004 53.370701 1 30013.237349 53.368668 1 40013.236707 53.366633 1 500... ... ... ...Based on the GCP-lists, I geocoded few image rasters manually using ENVI > Map > Registration > Warp from GCPs: Image to Map. That worked well but needs a lot of time.
Then, I tried to batch this processing chain in IDL, but I got already stuck at the first step:The import of the GCPs via
e=envi()GCP.pts via GCPs = ENVIGCPSet('GCP.pts') does not work. I got no error message, but the variable GCPs is a Null object.
Then I tried to create the GCP in IDL:
GCPs = ENVIGCPSet() ;open empty arrayGCPs.Add, X=13.239339, Y=53.374766, Z=0.0D, COLUMN=1, ROW=100, NAME='Geocod' ;add GCP 1GCPs.Add, X=13.238695, Y=53.372731, Z=0.0D, COLUMN=1, ROW=200, NAME='Geocod' ;add GCP 2GCPs.Add, X=13.238004, Y=53.370701, Z=0.0D, COLUMN=1, ROW=300, NAME='Geocod' ;add GCP 2GCPs.Export, 'D:\GCP_IDL.pts' However, no .pts file was created and I got no error message:
Do you have any hints concerning my problems with ENVIGCPSet? So far I couldn't find a IDL function for Warp from GCPs: Image to Map. Any hints there?Or is there maybe an much easier way to bring geoinformation and image together?
I would prefere to stay in IDL, but solutions in R are welcome, too.
أكثر...
I'm working with IDL 8.4/ENVI5.2 and format of the rasters is binary (.dat). But I could save it in any format.
What did I try so far:I wrote lists of GCPs based on the Lat-and Long-rasters (I only used every 100th pixel and not every pixel because of the size of the image). A list looks like this:
; ENVI Image to Map GCP File; projection info = {Geographic Lat/Lon, WGS-84, units=Degrees}; warp file: test.dat; Map (x,y), Image (x,y);13.239339 53.374766 1 10013.238695 53.372731 1 20013.238004 53.370701 1 30013.237349 53.368668 1 40013.236707 53.366633 1 500... ... ... ...Based on the GCP-lists, I geocoded few image rasters manually using ENVI > Map > Registration > Warp from GCPs: Image to Map. That worked well but needs a lot of time.
Then, I tried to batch this processing chain in IDL, but I got already stuck at the first step:The import of the GCPs via
e=envi()GCP.pts via GCPs = ENVIGCPSet('GCP.pts') does not work. I got no error message, but the variable GCPs is a Null object.
Then I tried to create the GCP in IDL:
GCPs = ENVIGCPSet() ;open empty arrayGCPs.Add, X=13.239339, Y=53.374766, Z=0.0D, COLUMN=1, ROW=100, NAME='Geocod' ;add GCP 1GCPs.Add, X=13.238695, Y=53.372731, Z=0.0D, COLUMN=1, ROW=200, NAME='Geocod' ;add GCP 2GCPs.Add, X=13.238004, Y=53.370701, Z=0.0D, COLUMN=1, ROW=300, NAME='Geocod' ;add GCP 2GCPs.Export, 'D:\GCP_IDL.pts' However, no .pts file was created and I got no error message:
Do you have any hints concerning my problems with ENVIGCPSet? So far I couldn't find a IDL function for Warp from GCPs: Image to Map. Any hints there?Or is there maybe an much easier way to bring geoinformation and image together?
I would prefere to stay in IDL, but solutions in R are welcome, too.
أكثر...