UTM Projection of DEM-data with gdalwarpgdalwarp with -srcwinHow can I fix badly specified geographic coordinate system with units in arc-seconds?Cropping Geotiff with gdalwarpProperly Converting from UTM to Lat/Long with gdalwarpgdalwarp reprojection causes missing datagdalwarp malfunctioning for GEOS projectionClipping raster takes forever in QGISgdalwarp with DSM as DEM leaves shifted pixels in their original unshifted locationTransforming Lambert Conformal Conic to EPSG:3857GDALwarp HDF5 Image With Unknown Projection
Can town administrative "code" overule state laws like those forbidding trespassing?
Circuitry of TV splitters
How do you conduct xenoanthropology after first contact?
Infinite past with a beginning?
Can an x86 CPU running in real mode be considered to be basically an 8086 CPU?
What is the offset in a seaplane's hull?
What are these boxed doors outside store fronts in New York?
Calculus Optimization - Point on graph closest to given point
How can I fix this gap between bookcases I made?
Why did the Germans forbid the possession of pet pigeons in Rostov-on-Don in 1941?
A Journey Through Space and Time
Can a German sentence have two subjects?
Why has Russell's definition of numbers using equivalence classes been finally abandoned? ( If it has actually been abandoned).
Banach space and Hilbert space topology
Why is the design of haulage companies so “special”?
Can I make popcorn with any corn?
Motorized valve interfering with button?
Prevent a directory in /tmp from being deleted
Patience, young "Padovan"
Is Social Media Science Fiction?
Email Account under attack (really) - anything I can do?
Is there really no realistic way for a skeleton monster to move around without magic?
How to type dʒ symbol (IPA) on Mac?
least quadratic residue under GRH: an EXPLICIT bound
UTM Projection of DEM-data with gdalwarp
gdalwarp with -srcwinHow can I fix badly specified geographic coordinate system with units in arc-seconds?Cropping Geotiff with gdalwarpProperly Converting from UTM to Lat/Long with gdalwarpgdalwarp reprojection causes missing datagdalwarp malfunctioning for GEOS projectionClipping raster takes forever in QGISgdalwarp with DSM as DEM leaves shifted pixels in their original unshifted locationTransforming Lambert Conformal Conic to EPSG:3857GDALwarp HDF5 Image With Unknown Projection
.everyoneloves__top-leaderboard:empty,.everyoneloves__mid-leaderboard:empty,.everyoneloves__bot-mid-leaderboard:empty margin-bottom:0;
I have downloaded a certain tile from https://www.eorc.jaxa.jp/ALOS/en/aw3d30/data/index.htm, and after the requisite logistics (unzipping, putting all *.AVE_DSM.tif
files into a single folder), I built a virtual
dataset with gdalbuildvrt *.tif asia.vrt
.
Now gdalinfo asia.vrt
gives this:
Size is 108000, 90000
Coordinate System is:
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"]]
Origin = (89.999999999999957,50.000000000000028)
Pixel Size = (0.000277777777778,-0.000277777777778)
Corner Coordinates:
Upper Left ( 90.0000000, 50.0000000) ( 90d 0' 0.00"E, 50d 0' 0.00"N)
Lower Left ( 90.0000000, 25.0000000) ( 90d 0' 0.00"E, 25d 0' 0.00"N)
Upper Right ( 120.0000000, 50.0000000) (120d 0' 0.00"E, 50d 0' 0.00"N)
Lower Right ( 120.0000000, 25.0000000) (120d 0' 0.00"E, 25d 0' 0.00"N)
Center ( 105.0000000, 37.5000000) (105d 0' 0.00"E, 37d30' 0.00"N)
Band 1 Block=128x128 Type=Int16, ColorInterp=Gray
I can query the data as follows:gdallocationinfo -geoloc asia.vrt 110.00 30.00
, meaning "Find the elevation
of the point at 110 degrees East and 30 degrees North",
and get a proper
Report:
Location: (72000P,72000L)
Band 1:
<LocationInfo><File>N029E110_AVE_DSM.tif</File></LocationInfo>
Value: 894
Now I would like to project the entire region I've got here
to a rectangle P = [a,b]x[c,d]
, so that I can query not with lat/lon pairs but
rather with a point from P
. Naturally, I need the distance
between the two points be somewhat reflective of their "original" distance.
I do not care much about the specific values of a,b,c,d
, and anyway I think
they depend on the units
option we would supply to gdalwarp
.
Now my question is,
what is the exact gdalwarp
-command I need to issue to accomplish this?
I know it has to do with proj
-string, but I am lost to its multitude of parameters.
dem gdalwarp
add a comment |
I have downloaded a certain tile from https://www.eorc.jaxa.jp/ALOS/en/aw3d30/data/index.htm, and after the requisite logistics (unzipping, putting all *.AVE_DSM.tif
files into a single folder), I built a virtual
dataset with gdalbuildvrt *.tif asia.vrt
.
Now gdalinfo asia.vrt
gives this:
Size is 108000, 90000
Coordinate System is:
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"]]
Origin = (89.999999999999957,50.000000000000028)
Pixel Size = (0.000277777777778,-0.000277777777778)
Corner Coordinates:
Upper Left ( 90.0000000, 50.0000000) ( 90d 0' 0.00"E, 50d 0' 0.00"N)
Lower Left ( 90.0000000, 25.0000000) ( 90d 0' 0.00"E, 25d 0' 0.00"N)
Upper Right ( 120.0000000, 50.0000000) (120d 0' 0.00"E, 50d 0' 0.00"N)
Lower Right ( 120.0000000, 25.0000000) (120d 0' 0.00"E, 25d 0' 0.00"N)
Center ( 105.0000000, 37.5000000) (105d 0' 0.00"E, 37d30' 0.00"N)
Band 1 Block=128x128 Type=Int16, ColorInterp=Gray
I can query the data as follows:gdallocationinfo -geoloc asia.vrt 110.00 30.00
, meaning "Find the elevation
of the point at 110 degrees East and 30 degrees North",
and get a proper
Report:
Location: (72000P,72000L)
Band 1:
<LocationInfo><File>N029E110_AVE_DSM.tif</File></LocationInfo>
Value: 894
Now I would like to project the entire region I've got here
to a rectangle P = [a,b]x[c,d]
, so that I can query not with lat/lon pairs but
rather with a point from P
. Naturally, I need the distance
between the two points be somewhat reflective of their "original" distance.
I do not care much about the specific values of a,b,c,d
, and anyway I think
they depend on the units
option we would supply to gdalwarp
.
Now my question is,
what is the exact gdalwarp
-command I need to issue to accomplish this?
I know it has to do with proj
-string, but I am lost to its multitude of parameters.
dem gdalwarp
2
Look at-te
and-te_srs
. Those will let you specify the target extents of your rectangle in whatever spatial reference you want
– mikewatt
Apr 3 at 23:32
add a comment |
I have downloaded a certain tile from https://www.eorc.jaxa.jp/ALOS/en/aw3d30/data/index.htm, and after the requisite logistics (unzipping, putting all *.AVE_DSM.tif
files into a single folder), I built a virtual
dataset with gdalbuildvrt *.tif asia.vrt
.
Now gdalinfo asia.vrt
gives this:
Size is 108000, 90000
Coordinate System is:
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"]]
Origin = (89.999999999999957,50.000000000000028)
Pixel Size = (0.000277777777778,-0.000277777777778)
Corner Coordinates:
Upper Left ( 90.0000000, 50.0000000) ( 90d 0' 0.00"E, 50d 0' 0.00"N)
Lower Left ( 90.0000000, 25.0000000) ( 90d 0' 0.00"E, 25d 0' 0.00"N)
Upper Right ( 120.0000000, 50.0000000) (120d 0' 0.00"E, 50d 0' 0.00"N)
Lower Right ( 120.0000000, 25.0000000) (120d 0' 0.00"E, 25d 0' 0.00"N)
Center ( 105.0000000, 37.5000000) (105d 0' 0.00"E, 37d30' 0.00"N)
Band 1 Block=128x128 Type=Int16, ColorInterp=Gray
I can query the data as follows:gdallocationinfo -geoloc asia.vrt 110.00 30.00
, meaning "Find the elevation
of the point at 110 degrees East and 30 degrees North",
and get a proper
Report:
Location: (72000P,72000L)
Band 1:
<LocationInfo><File>N029E110_AVE_DSM.tif</File></LocationInfo>
Value: 894
Now I would like to project the entire region I've got here
to a rectangle P = [a,b]x[c,d]
, so that I can query not with lat/lon pairs but
rather with a point from P
. Naturally, I need the distance
between the two points be somewhat reflective of their "original" distance.
I do not care much about the specific values of a,b,c,d
, and anyway I think
they depend on the units
option we would supply to gdalwarp
.
Now my question is,
what is the exact gdalwarp
-command I need to issue to accomplish this?
I know it has to do with proj
-string, but I am lost to its multitude of parameters.
dem gdalwarp
I have downloaded a certain tile from https://www.eorc.jaxa.jp/ALOS/en/aw3d30/data/index.htm, and after the requisite logistics (unzipping, putting all *.AVE_DSM.tif
files into a single folder), I built a virtual
dataset with gdalbuildvrt *.tif asia.vrt
.
Now gdalinfo asia.vrt
gives this:
Size is 108000, 90000
Coordinate System is:
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"]]
Origin = (89.999999999999957,50.000000000000028)
Pixel Size = (0.000277777777778,-0.000277777777778)
Corner Coordinates:
Upper Left ( 90.0000000, 50.0000000) ( 90d 0' 0.00"E, 50d 0' 0.00"N)
Lower Left ( 90.0000000, 25.0000000) ( 90d 0' 0.00"E, 25d 0' 0.00"N)
Upper Right ( 120.0000000, 50.0000000) (120d 0' 0.00"E, 50d 0' 0.00"N)
Lower Right ( 120.0000000, 25.0000000) (120d 0' 0.00"E, 25d 0' 0.00"N)
Center ( 105.0000000, 37.5000000) (105d 0' 0.00"E, 37d30' 0.00"N)
Band 1 Block=128x128 Type=Int16, ColorInterp=Gray
I can query the data as follows:gdallocationinfo -geoloc asia.vrt 110.00 30.00
, meaning "Find the elevation
of the point at 110 degrees East and 30 degrees North",
and get a proper
Report:
Location: (72000P,72000L)
Band 1:
<LocationInfo><File>N029E110_AVE_DSM.tif</File></LocationInfo>
Value: 894
Now I would like to project the entire region I've got here
to a rectangle P = [a,b]x[c,d]
, so that I can query not with lat/lon pairs but
rather with a point from P
. Naturally, I need the distance
between the two points be somewhat reflective of their "original" distance.
I do not care much about the specific values of a,b,c,d
, and anyway I think
they depend on the units
option we would supply to gdalwarp
.
Now my question is,
what is the exact gdalwarp
-command I need to issue to accomplish this?
I know it has to do with proj
-string, but I am lost to its multitude of parameters.
dem gdalwarp
dem gdalwarp
asked Apr 3 at 18:09
IlonpilaajaIlonpilaaja
1366
1366
2
Look at-te
and-te_srs
. Those will let you specify the target extents of your rectangle in whatever spatial reference you want
– mikewatt
Apr 3 at 23:32
add a comment |
2
Look at-te
and-te_srs
. Those will let you specify the target extents of your rectangle in whatever spatial reference you want
– mikewatt
Apr 3 at 23:32
2
2
Look at
-te
and -te_srs
. Those will let you specify the target extents of your rectangle in whatever spatial reference you want– mikewatt
Apr 3 at 23:32
Look at
-te
and -te_srs
. Those will let you specify the target extents of your rectangle in whatever spatial reference you want– mikewatt
Apr 3 at 23:32
add a comment |
1 Answer
1
active
oldest
votes
It appears that the answer to my question is about this:gdalwarp -t_srs '+proj=eqc +lat_0=37.5 +lon_0=105.0 +units=m' *.tif result.tif
add a comment |
Your Answer
StackExchange.ready(function()
var channelOptions =
tags: "".split(" "),
id: "79"
;
initTagRenderer("".split(" "), "".split(" "), channelOptions);
StackExchange.using("externalEditor", function()
// Have to fire editor after snippets, if snippets enabled
if (StackExchange.settings.snippets.snippetsEnabled)
StackExchange.using("snippets", function()
createEditor();
);
else
createEditor();
);
function createEditor()
StackExchange.prepareEditor(
heartbeatType: 'answer',
autoActivateHeartbeat: false,
convertImagesToLinks: false,
noModals: true,
showLowRepImageUploadWarning: true,
reputationToPostImages: null,
bindNavPrevention: true,
postfix: "",
imageUploader:
brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
allowUrls: true
,
onDemand: true,
discardSelector: ".discard-answer"
,immediatelyShowMarkdownHelp:true
);
);
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function ()
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fgis.stackexchange.com%2fquestions%2f317688%2futm-projection-of-dem-data-with-gdalwarp%23new-answer', 'question_page');
);
Post as a guest
Required, but never shown
1 Answer
1
active
oldest
votes
1 Answer
1
active
oldest
votes
active
oldest
votes
active
oldest
votes
It appears that the answer to my question is about this:gdalwarp -t_srs '+proj=eqc +lat_0=37.5 +lon_0=105.0 +units=m' *.tif result.tif
add a comment |
It appears that the answer to my question is about this:gdalwarp -t_srs '+proj=eqc +lat_0=37.5 +lon_0=105.0 +units=m' *.tif result.tif
add a comment |
It appears that the answer to my question is about this:gdalwarp -t_srs '+proj=eqc +lat_0=37.5 +lon_0=105.0 +units=m' *.tif result.tif
It appears that the answer to my question is about this:gdalwarp -t_srs '+proj=eqc +lat_0=37.5 +lon_0=105.0 +units=m' *.tif result.tif
answered Apr 4 at 3:36
IlonpilaajaIlonpilaaja
1366
1366
add a comment |
add a comment |
Thanks for contributing an answer to Geographic Information Systems Stack Exchange!
- Please be sure to answer the question. Provide details and share your research!
But avoid …
- Asking for help, clarification, or responding to other answers.
- Making statements based on opinion; back them up with references or personal experience.
To learn more, see our tips on writing great answers.
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function ()
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fgis.stackexchange.com%2fquestions%2f317688%2futm-projection-of-dem-data-with-gdalwarp%23new-answer', 'question_page');
);
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
2
Look at
-te
and-te_srs
. Those will let you specify the target extents of your rectangle in whatever spatial reference you want– mikewatt
Apr 3 at 23:32