nav[aria-label="Primary Navigation"] { padding: 0; & ul { list-style: none; width: 100%; display: flex; flex-direction: row; justify-content: start; align-items: start; gap: 30px; padding: 0; & li { margin: 0; } & ul li { list-style: none; } } }

Projecting multiple images from one coordinate system to another,

Anonymous User

Original User: Jamal432@gmail.com

Projecting multiple images from one coordinate system to another,

I wanted to project multiple Jpeg images form Israel to Palestine coordinate system using the batch of the ???project raster??? tool but couldn???t figure out how the input raster are populated automatically in the ???output raster dataset???

[ATTACH=CONFIG]27794[/ATTACH], [ATTACH=CONFIG]27795[/ATTACH]

Is there a way to populate the images automatically in the ???output raster dataset????

Thank you

Best

Jamal

Clip_38.jpg

Clip_39.jpg

タグ付けされた投稿をもっと見る

Accepted answers

Anonymous User

Original User: JSkinn3

Hi Jamal,

I don't believe there is a way to auto populate this parameter. You can easily accomplish this in python, though. Here is an example:

import arcpy from arcpy import env env.workspace = r"D:\Practical1_AddingData\Data"  for raster in arcpy.ListRasters("*.jpg"):     arcpy.ProjectRaster_management(raster, raster + "_project.jpg", "GEOGCS['GCS_Palestine_1923',DATUM['D_Palestine_1923',SPHEROID['Clarke_1880_Benoit',6378300.789,293.4663155389802]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]","NEAREST","0.0000025","Palestine_1923_To_Israel_1","#","GEOGCS['GCS_Israel',DATUM['D_Israel',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]")

The above code will iterate through the 'D:\Practical1_AddingData\Data' directory and reproject each JPG to a new JPG with the same name, just with "_project" added to it.

All comments

Anonymous User

Original User: GSCUser85

Just had a look at the site and downloaded one of the daily files.
No spatial reference info....
Even the extent :
Top 40.05
Bot -40.05
Left -20.05
Right 55.05
Doesn't look like anything...
Is there any other info available on the site telling you where this stuff is in the world?
Cheers,
N

PrabinRokaya

Dear Neil,

Thank you for the kind assistance. I downloaded this precipitation images from Famine Early Warning Systems Network of US Government http://earlywarning.usgs.gov/fews/africa/web/dwndailyrfe.php . When I open the file in ArcGIS and see spatial reference in properties, I so no reference. Is there any way out?

Prabin,
errr...
If the original bil images are not georeferenced then no amount of assigning a coordinate system / projecting and what not is going to help.
They are simply not georeferenced. You will have to go through a georeferencing step which involves picking control points from the image and assigning to these points their correct coordinates.
But that is rather odd, where did this raster data come from?
Neil

Anonymous User

Original User: Jamal432@gmail.com

Thank you Neil for the assistance. It was really kind of you. But yes, the bil files are not correctly geo-referenced. I also tried to do it via model builder but because file are not properly geo-referenced, it was difficult. Actually I have precipitation images which I want to first geo-referenced (usually I do it via arccatalog, going to spatial reference, edit and selecting GCS WGS 1984) and then through data management tool, project raster, I project them into WGS_1984_UTM_Zone_37N. Then in ArcGIS, I use extract by mask tool to get precipitation map for only my watershed. And calculate value using raster calculator. But the problem I have is, I have thousands of images and I can not do it manually. I wanted to do it via model builder or python and I am looking for assistance as all my attempts to figure it out have failed. It would be really really grateful if you could help me.

Thank you for your time and consideration.
Prabin Rokaya
MSC, UNESCO-IHE Institute for Water Education

Hi Neil,

I�??m not sure if this code can be developed in form of geoprocessing tool. This would be much more easier for the end users. This command is a need for everybody

Also, I would like to point out that the code originally provided by Jake includes all parameters required for the projection.

Best

Jamal

NeilAyres

Prabin,
errr...
If the original bil images are not georeferenced then no amount of assigning a coordinate system / projecting and what not is going to help.
They are simply not georeferenced. You will have to go through a georeferencing step which involves picking control points from the image and assigning to these points their correct coordinates.
But that is rather odd, where did this raster data come from?
Neil

Anonymous User

Original User: prabinrokaya@gmail.com

Thank you Neil for the assistance. It was really kind of you. But yes, the bil files are not correctly geo-referenced. I also tried to do it via model builder but because file are not properly geo-referenced, it was difficult. Actually I have precipitation images which I want to first geo-referenced (usually I do it via arccatalog, going to spatial reference, edit and selecting GCS WGS 1984) and then through data management tool, project raster, I project them into WGS_1984_UTM_Zone_37N. Then in ArcGIS, I use extract by mask tool to get precipitation map for only my watershed. And calculate value using raster calculator. But the problem I have is, I have thousands of images and I can not do it manually. I wanted to do it via model builder or python and I am looking for assistance as all my attempts to figure it out have failed. It would be really really grateful if you could help me.

Thank you for your time and consideration.
Prabin Rokaya
MSC, UNESCO-IHE Institute for Water Education

Are all of your bil files correctly georeferenced and on the WGS84 datum?
If so, follow the code template given by Jake above.

The WKID for WGS_1984_UTM_Zone_37N is 32637, so, rather than listing the entire text for the projection, I would use the WKID to create a spatial reference object for use in the Project function.
Thus :
import arcpy
from arcpy import env
env.workspace = r"Your path to the data here"

srOut = arcpy.SpatialReference(32637)


for raster in arcpy.ListRasters("*.bil"):
    arcpy.ProjectRaster_management(raster, raster + "_project.bil", srOut)
See the examples here :
http://resources.arcgis.com/en/help/main/10.2/#/Project_Raster/00170000007q000000/

Cheers and good luck,
Neil

NeilAyres

Are all of your bil files correctly georeferenced and on the WGS84 datum?
If so, follow the code template given by Jake above.

The WKID for WGS_1984_UTM_Zone_37N is 32637, so, rather than listing the entire text for the projection, I would use the WKID to create a spatial reference object for use in the Project function.
Thus :

import arcpy
from arcpy import env
env.workspace = r"Your path to the data here"

srOut = arcpy.SpatialReference(32637)


for raster in arcpy.ListRasters("*.bil"):
    arcpy.ProjectRaster_management(raster, raster + "_project.bil", srOut)

See the examples here :
http://resources.arcgis.com/en/help/main/10.2/#/Project_Raster/00170000007q000000/

Cheers and good luck,
Neil

Anonymous User

Original User: prabinrokaya@gmail.com

Can anyone tell me how can I project thousands of bil images into WGS_1984_UTM_Zone_37N?

JamalNUMAN

Open your python window in ArcMap and you can copy/paste the code in there.

The code will iterate through your 'D:\Practical1_Adding\Data directory and reproject each JPG. It will write the new projected raster to the same directory with the same name plus "_project" added to it. For example, the image 'yosh_09-06-12_77.jpg' will be reprojected and named 'yosh_09-06-12_77_project.jpg'.

Thank you Jake very much for the prompt help. It works fantastically.

While the code you have already sent is great in its current form, I still couldn�??t figure out why the behavior of the �??project raster�?� can�??t be developed to be the same as �??raster to other format�?� such that multiple projected raster can be populated automatically in the specified folder.

[ATTACH=CONFIG]27818[/ATTACH]

Meanwhile, how the coded you have already sent can be further developed to allow the user to specify the input and output of raster (on the hard drive)?
Is it also possible to be developed in form of script tool?

[ATTACH=CONFIG]27819[/ATTACH]

Clip_50.jpg

Clip_51.jpg

Anonymous User

Original User: JSkinn3

Open your python window in ArcMap and you can copy/paste the code in there.

The code will iterate through your 'D:\Practical1_Adding\Data directory and reproject each JPG. It will write the new projected raster to the same directory with the same name plus "_project" added to it. For example, the image 'yosh_09-06-12_77.jpg' will be reprojected and named 'yosh_09-06-12_77_project.jpg'.

JamalNUMAN

Hi Jamal,

I don't believe there is a way to auto populate this parameter. You can easily accomplish this in python, though. Here is an example:
import arcpy
from arcpy import env
env.workspace = r"D:\Practical1_AddingData\Data"

for raster in arcpy.ListRasters("*.jpg"):
    arcpy.ProjectRaster_management(raster, raster + "_project.jpg", "GEOGCS['GCS_Palestine_1923',DATUM['D_Palestine_1923',SPHEROID['Clarke_1880_Benoit',6378300.789,293.4663155389802]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]","NEAREST","0.0000025","Palestine_1923_To_Israel_1","#","GEOGCS['GCS_Israel',DATUM['D_Israel',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]")
The above code will iterate through the 'D:\Practical1_AddingData\Data' directory and reproject each JPG to a new JPG with the same name, just with "_project" added to it.

Many thanks Jake for the help,

Sorry but I couldn�??t know how to integrate the code you have provided with the �??project raster�?� tool. Could you please elaborate a bit more.

[ATTACH=CONFIG]27804[/ATTACH]

By the way, I like the behavior of the �??raster to other format�?� tool where the output is populated automatically in the specified folder. Why the design of the �??project raster�?� is not developed to be the same as �??raster to other format�?�?

[ATTACH=CONFIG]27805[/ATTACH]

Best

Jamal

Clip_42.jpg

Clip_41.jpg

Anonymous User

Original User: JSkinn3

Hi Jamal,

I don't believe there is a way to auto populate this parameter. You can easily accomplish this in python, though. Here is an example:

import arcpy from arcpy import env env.workspace = r"D:\Practical1_AddingData\Data"  for raster in arcpy.ListRasters("*.jpg"):     arcpy.ProjectRaster_management(raster, raster + "_project.jpg", "GEOGCS['GCS_Palestine_1923',DATUM['D_Palestine_1923',SPHEROID['Clarke_1880_Benoit',6378300.789,293.4663155389802]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]","NEAREST","0.0000025","Palestine_1923_To_Israel_1","#","GEOGCS['GCS_Israel',DATUM['D_Israel',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]")

The above code will iterate through the 'D:\Practical1_AddingData\Data' directory and reproject each JPG to a new JPG with the same name, just with "_project" added to it.

会話に参加する

サインインしたメンバーは投稿、更新のフォローなどができます。初めてですか？無料アカウントを登録してください。

クイックリンク

Community Help
Find useful guides, FAQs, and documents to help you navigate and make the most of Esri Community.