Hfrhyu

How can bays and straits be determined in a procedurally generated map?

Is it legal for company to use my work email to pretend I still work there?

Do I have a twin with permutated remainders?

Alternative to sending password over mail?

What is the word for reserving something for yourself before others do?

Does detail obscure or enhance action?

Why is consensus so controversial in Britain?

"You are your self first supporter", a more proper way to say it

What's the output of a record needle playing an out-of-speed record

Rock identification in KY

Perform and show arithmetic with LuaLaTeX

Convert two switches to a dual stack, and add outlet - possible here?

Cross compiling for RPi - error while loading shared libraries

How is the claim "I am in New York only if I am in America" the same as "If I am in New York, then I am in America?

Are astronomers waiting to see something in an image from a gravitational lens that they've already seen in an adjacent image?

NMaximize is not converging to a solution

Has there ever been an airliner design involving reducing generator load by installing solar panels?

Today is the Center

tikz convert color string to hex value

Maximum likelihood parameters deviate from posterior distributions

Why can't we play rap on piano?

Was any UN Security Council vote triple-vetoed?

Is it possible to run Internet Explorer on OS X El Capitan?

Important Resources for Dark Age Civilizations?

Calculating tile coordinates based on location and crs for ArcGIS maps?

How to georeference a web mercator tile correctly using gdal?What is wrong with this Proj4 transformation?proj.4 definition for Moscow GCS_Bessel_1841 to WGS84 convertionWarping raster to new domain/resolution/projection causes data lossConverting long, lat to tile of NASA GIBS - WGS 84 / Lat-lon / Geographic (EPSG:4326)?Numerical Vertex Edit Error: latitude or longitude exceeded limitsWhy is my reprojection from ESPG4326,WGS84 to ESPG3005,NAD83/BCAlbers failing?Is SRID 4326 Lon/Lat or Lat/Lon?Using project / layer CRS in QGIS?Google Maps - get all tiles locations around a specific latitude longitude

.everyoneloves__top-leaderboard:empty,.everyoneloves__mid-leaderboard:empty,.everyoneloves__bot-mid-leaderboard:empty margin-bottom:0;

0

I am working on a solution to download map tiles based on given location:(latitude, longitude) and zoom level

Followed https://wiki.openstreetmap.org/wiki/Slippy_map_tilenames

The key concept is to convert (latitude, longitude) to tile reference, eg. javascript code:

function long2tileX(lon,zoom) return (Math.floor((lon+180)/360*Math.pow(2,zoom)));

function lat2tileY(lat,zoom) return (Math.floor((1-Math.log(Math.tan(lat*Math.PI/180) + 1/Math.cos(lat*Math.PI/180))/Math.PI)/2 *Math.pow(2,zoom)));

For sample location (63.433772, 10.393456) and zoom level (17) this formulas return (69320, 35421) which is exactly what's needed to get the tile:

https://opencache.statkart.no/gatekeeper/gk/gk.open_gmaps?layers=topo4&zoom=17&x=69320&y=35421

Now I am struggling to achieve the same with ArcGIS maps.
The following CRS object is supplied:


"Name":"Geocache Basis",
"Crs":
"Resolutions":[
21674.7100160867,
10837.355008043351,
5418.6775040216753,
2709.3387520108377,
1354.6693760054188,
677.33468800270941,
338.66734400135471,
169.33367200067735,
84.666836000338677,
42.333418000169338,
21.166709000084669,
10.583354500042335,
5.2916772500211673,
2.6458386250105836,
1.3229193125052918,
0.66145965625264591,
0.33072982812632296,
0.16536491406316148
],
"OriginX":-2500000.0,
"OriginY":9045984.0,
"Srid":"EPSG:32633",
"Definition":"+proj=utm +zone=33 +ellps=WGS84 +datum=WGS84 +units=m +no_defs"
,
"Attribution":"Kartverket, Geovekst og kommuner - Geodata AS",
"MapProvider":"arcgis"


Tried to use proj4 and define EPSG:32633 with 'Definition' field value, then do a projection, but that didn't get me anywhere close to calculating
the respective tile coordinates, which I experimentally determined to be (47333, 65439)
https://services.geodataonline.no/arcgis/rest/services/Geocache_UTM33_WGS84/GeocacheBasis/MapServer/tile/17/47333/65439

What mathematical formula can I apply to follow the same scenario and get an accurate result?

coordinate-system tiles proj

share|improve this question

edited Apr 2 at 21:50

PolyGeo♦

53.9k1781245

asked Apr 2 at 12:56

Irek ZuchowskiIrek Zuchowski

1

New contributor

Irek Zuchowski is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.

add a comment | 

0

I am working on a solution to download map tiles based on given location:(latitude, longitude) and zoom level

Followed https://wiki.openstreetmap.org/wiki/Slippy_map_tilenames

The key concept is to convert (latitude, longitude) to tile reference, eg. javascript code:

function long2tileX(lon,zoom) return (Math.floor((lon+180)/360*Math.pow(2,zoom)));

function lat2tileY(lat,zoom) return (Math.floor((1-Math.log(Math.tan(lat*Math.PI/180) + 1/Math.cos(lat*Math.PI/180))/Math.PI)/2 *Math.pow(2,zoom)));

For sample location (63.433772, 10.393456) and zoom level (17) this formulas return (69320, 35421) which is exactly what's needed to get the tile:

https://opencache.statkart.no/gatekeeper/gk/gk.open_gmaps?layers=topo4&zoom=17&x=69320&y=35421

Now I am struggling to achieve the same with ArcGIS maps.
The following CRS object is supplied:


"Name":"Geocache Basis",
"Crs":
"Resolutions":[
21674.7100160867,
10837.355008043351,
5418.6775040216753,
2709.3387520108377,
1354.6693760054188,
677.33468800270941,
338.66734400135471,
169.33367200067735,
84.666836000338677,
42.333418000169338,
21.166709000084669,
10.583354500042335,
5.2916772500211673,
2.6458386250105836,
1.3229193125052918,
0.66145965625264591,
0.33072982812632296,
0.16536491406316148
],
"OriginX":-2500000.0,
"OriginY":9045984.0,
"Srid":"EPSG:32633",
"Definition":"+proj=utm +zone=33 +ellps=WGS84 +datum=WGS84 +units=m +no_defs"
,
"Attribution":"Kartverket, Geovekst og kommuner - Geodata AS",
"MapProvider":"arcgis"


Tried to use proj4 and define EPSG:32633 with 'Definition' field value, then do a projection, but that didn't get me anywhere close to calculating
the respective tile coordinates, which I experimentally determined to be (47333, 65439)
https://services.geodataonline.no/arcgis/rest/services/Geocache_UTM33_WGS84/GeocacheBasis/MapServer/tile/17/47333/65439

What mathematical formula can I apply to follow the same scenario and get an accurate result?

coordinate-system tiles proj

share|improve this question

edited Apr 2 at 21:50

PolyGeo♦

53.9k1781245

asked Apr 2 at 12:56

Irek ZuchowskiIrek Zuchowski

1

New contributor

Irek Zuchowski is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.

add a comment | 

I am working on a solution to download map tiles based on given location:(latitude, longitude) and zoom level

Followed https://wiki.openstreetmap.org/wiki/Slippy_map_tilenames

The key concept is to convert (latitude, longitude) to tile reference, eg. javascript code:

function long2tileX(lon,zoom) return (Math.floor((lon+180)/360*Math.pow(2,zoom)));

function lat2tileY(lat,zoom) return (Math.floor((1-Math.log(Math.tan(lat*Math.PI/180) + 1/Math.cos(lat*Math.PI/180))/Math.PI)/2 *Math.pow(2,zoom)));

For sample location (63.433772, 10.393456) and zoom level (17) this formulas return (69320, 35421) which is exactly what's needed to get the tile:

https://opencache.statkart.no/gatekeeper/gk/gk.open_gmaps?layers=topo4&zoom=17&x=69320&y=35421

Now I am struggling to achieve the same with ArcGIS maps.
The following CRS object is supplied:


"Name":"Geocache Basis",
"Crs":
"Resolutions":[
21674.7100160867,
10837.355008043351,
5418.6775040216753,
2709.3387520108377,
1354.6693760054188,
677.33468800270941,
338.66734400135471,
169.33367200067735,
84.666836000338677,
42.333418000169338,
21.166709000084669,
10.583354500042335,
5.2916772500211673,
2.6458386250105836,
1.3229193125052918,
0.66145965625264591,
0.33072982812632296,
0.16536491406316148
],
"OriginX":-2500000.0,
"OriginY":9045984.0,
"Srid":"EPSG:32633",
"Definition":"+proj=utm +zone=33 +ellps=WGS84 +datum=WGS84 +units=m +no_defs"
,
"Attribution":"Kartverket, Geovekst og kommuner - Geodata AS",
"MapProvider":"arcgis"


Tried to use proj4 and define EPSG:32633 with 'Definition' field value, then do a projection, but that didn't get me anywhere close to calculating
the respective tile coordinates, which I experimentally determined to be (47333, 65439)
https://services.geodataonline.no/arcgis/rest/services/Geocache_UTM33_WGS84/GeocacheBasis/MapServer/tile/17/47333/65439

What mathematical formula can I apply to follow the same scenario and get an accurate result?

coordinate-system tiles proj

share|improve this question

edited Apr 2 at 21:50

PolyGeo♦

53.9k1781245

asked Apr 2 at 12:56

Irek ZuchowskiIrek Zuchowski

1

New contributor

Irek Zuchowski is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.

share|improve this question

edited Apr 2 at 21:50

PolyGeo♦

53.9k1781245

asked Apr 2 at 12:56

Irek ZuchowskiIrek Zuchowski

1

New contributor

Irek Zuchowski is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.

share|improve this question

edited Apr 2 at 21:50

PolyGeo♦

53.9k1781245

asked Apr 2 at 12:56

Irek ZuchowskiIrek Zuchowski

1

New contributor

Irek Zuchowski is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.

edited Apr 2 at 21:50

PolyGeo♦

53.9k1781245

edited Apr 2 at 21:50

PolyGeo♦

53.9k1781245

asked Apr 2 at 12:56

Irek ZuchowskiIrek Zuchowski

1

New contributor

Irek Zuchowski is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.

asked Apr 2 at 12:56

Irek ZuchowskiIrek Zuchowski

1