<div dir="ltr">Mark,<div><br></div><div>Rather than use the proj.db to look up an ESRI or EPSG code you could try explicitly defining your source coordinate system. There's an example here:</div><div><a href="https://mapserver.org/mapfile/projection.html">https://mapserver.org/mapfile/projection.html</a><br></div><div><br></div><div>You can also edit, or create your own epsg file. In fact, that used to be a recommended optimization so that proj didn't have to parse a big long epsg file. It's just an ASCII text file. On a Linux computer the default location is in /usr/share/proj/. </div><div><br></div><div>And as a test, on the server where you're running MapServer, can you do something like:</div><div> ogr2ogr out.shp source.shp -t_srs epsg:3857</div><div>This would tell you if the installed gdal and proj libraries recognise the coordinate system of the source.shp file. Or simpler still, ogrinfo source.shp and see if the coordinate system is identified.</div><div><br></div><div>If you don't have a .prj file accompanying your source.shp you can declare the source coordinate system with -s_srs esri:103749 or possibly -s_srs epsg:103749 which might give you some clues as to what's working and not working on your computer.</div><div><br></div><div>Rich</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Jan 7, 2022 at 8:26 AM Mark Volz <<a href="mailto:MarkVolz@co.lyon.mn.us">MarkVolz@co.lyon.mn.us</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div lang="EN-US">
<div class="gmail-m_-8207557529120665152WordSection1">
<p class="MsoNormal"><span style="color:rgb(31,73,125)">Hello,<u></u><u></u></span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125)"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125)">I apologize please ignore my mention of EPSG 4326. My end result is a WMS map in Web Mercator 3857. My underlying data is stored in a shapefile that is in WKID 103749 which is a projection has already been
included in proj.db as ESRI as the authority . The underlying map application is using javascript. I know that my map is finding proj.db because it MapServer does not complain when I convert UTM 14 to Web Mercator… but that does not help because the underlying
data is not in UTM 14.<u></u><u></u></span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125)"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125)">Here is how I have the projections of my map is set up:<u></u><u></u></span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125)">MAP<u></u><u></u></span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125)"> PROJECTION<u></u><u></u></span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125)"> "init=epsg:3857"<u></u><u></u></span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125)"> END #end map projection<u></u><u></u></span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125)"> LAYER<u></u><u></u></span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125)"> DATA 'ParcelsLyon.shp'<u></u><u></u></span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125)"> PROJECTION<u></u><u></u></span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125)"> "init=epsg:103749" #Error: proj_create: cannot expand init=espg:103749 type=crs<u></u><u></u></span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125)"> "init=crs:103749" #</span>
<span style="color:rgb(31,73,125)">Error: proj_create: cannot expand init=crs:103749 type=crs<u></u><u></u></span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125)"> "init=ESRI:103749" #Error: proj_create: cannot expand init=esri:103749 type=crs<u></u><u></u></span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125)"> "init=epsg:26914" #MapServer does not warn about any errors, but of course nothing draws because my data is in a different projection<u></u><u></u></span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125)"> END #layer projection<u></u><u></u></span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125)"> END #end layer<u></u><u></u></span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125)">END #end map<u></u><u></u></span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125)"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125)">Example request<u></u><u></u></span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125)"><a href="http://localhost/cgi-bin/mapserv.exe?SERVICE=WMS&VERSION=1.1.1&REQUEST=GetMap&FORMAT=image%2Fpng&TRANSPARENT=true&LAYERS=parcels%2CCityLimits%2CTownships&MAP=%2Fms4w%2Fapps%2Flyondata%2F.%2Fmaps%2Falllayers%2FAllLayers.map&SRS=EPSG%3A3857&STYLES=&WIDTH=1569&HEIGHT=912&BBOX=-10678090.52943604%2C5544408.3738651555%2C-10668156.653791208%2C5550182.55779638" target="_blank">http://localhost/cgi-bin/mapserv.exe?SERVICE=WMS&VERSION=1.1.1&REQUEST=GetMap&FORMAT=image%2Fpng&TRANSPARENT=true&LAYERS=parcels%2CCityLimits%2CTownships&MAP=%2Fms4w%2Fapps%2Flyondata%2F.%2Fmaps%2Falllayers%2FAllLayers.map&SRS=EPSG%3A3857&STYLES=&WIDTH=1569&HEIGHT=912&BBOX=-10678090.52943604%2C5544408.3738651555%2C-10668156.653791208%2C5550182.55779638</a><u></u><u></u></span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125)"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125)"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125)">@Bob<u></u><u></u></span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125)">I could project my shapefile into any projection. However, I will hit a wall with my raster layers which I prefer not to project at this time.<u></u><u></u></span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125)"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125)">Thanks!<u></u><u></u></span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125)"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125)"><u></u> <u></u></span></p>
<div>
<div style="border-right:none;border-bottom:none;border-left:none;border-top:1pt solid rgb(225,225,225);padding:3pt 0in 0in">
<p class="MsoNormal"><b>From:</b> <a href="mailto:artclifford@gmail.com" target="_blank">artclifford@gmail.com</a> <<a href="mailto:artclifford@gmail.com" target="_blank">artclifford@gmail.com</a>> <br>
<b>Sent:</b> Thursday, January 6, 2022 7:19 PM<br>
<b>To:</b> Mark Volz <<a href="mailto:MarkVolz@co.lyon.mn.us" target="_blank">MarkVolz@co.lyon.mn.us</a>><br>
<b>Subject:</b> RE: [mapserver-users] working with ESRI coordinate systems in proj.db<u></u><u></u></p>
</div>
</div>
<div>
<p class="MsoNormal"><span style="color:rgb(31,73,125)"><u></u> <u></u></span></p>
<p class="MsoNormal">It’s a little unclear to me what you are attempting.<u></u><u></u></p>
<p class="MsoNormal">A mapfile can be defined where the map’s projection is EPSG:3857 (web mercator) set in a projection block at the map level, and your layer definition can have a projection block set to whatever projection you are using. That tells the map
to render in web Mercator but says what projection your layer is stored in so that map server knows when to reproject stuff. I think you know that, but I’m just covering the “is it plugged in” issues and to suggest making sure the right projections are in
the right projection blocks. <u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Although not quite the problem you are having yet, epsg:4326 is wgs84 lat/lon. Web Mercator is: EPSG:3857.<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">If I remember correctly the projection can be defined for a mapfile projection block using the proj code definition if not the ESPG/numeric code. If that works, it might mean that MapServer isn’t using the proj.db file like you expect.
Or if it is using a proj.db file it may be located elsewhere. Whether putting the whole proj definition is viable for your ultimate use it would at least identify the problem.<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Then there’s a question of access. How are you “accessing” things? Are you using WFS? Or are you using a python/php/other script? Are things stored in a shp file or postgis or other database?
<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">If you are using wfs are you passing the desired projection on the request? Also, are you defining the allowed wms/wfs projections?<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">I’m not sure how much I personally can be of help, but any further context on your setup and what you are attempting that you are free and willing to provide will help folks here narrow down the possibilities.<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<div style="border-right:none;border-bottom:none;border-left:none;border-top:1pt solid rgb(225,225,225);padding:3pt 0in 0in">
<p class="MsoNormal"><b>From:</b> MapServer-users <<a href="mailto:mapserver-users-bounces@lists.osgeo.org" target="_blank">mapserver-users-bounces@lists.osgeo.org</a>>
<b>On Behalf Of </b>Mark Volz<br>
<b>Sent:</b> Thursday, January 6, 2022 9:19 AM<br>
<b>To:</b> <a href="mailto:mapserver-users@lists.osgeo.org" target="_blank">mapserver-users@lists.osgeo.org</a><br>
<b>Subject:</b> [mapserver-users] working with ESRI coordinate systems in proj.db<u></u><u></u></p>
</div>
</div>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Hello,<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Could someone help me with setting up MapServer to use a non ESRI projection that is already stored in proj.db? The coordinate system is for Lyon County, MN #103749 and is stored with ESRI as the authority. I know that we can set up MapServer
to use Web Mercator by setting the adding "init=epsg:4326" to a projection block. Unfortunately I cannot access my counties coordinate system by setting the projection with “init=esri:103749”.<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Thank you for any help!<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125)">Sincerely,<u></u><u></u></span></p>
<p class="MsoNormal"><b><span style="color:rgb(31,73,125)">Mark Volz, GISP</span></b><b><span style="font-size:12pt;color:rgb(31,73,125)"><u></u><u></u></span></b></p>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div>
</div>
_______________________________________________<br>
MapServer-users mailing list<br>
<a href="mailto:MapServer-users@lists.osgeo.org" target="_blank">MapServer-users@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/mapserver-users" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/mapserver-users</a><br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr">Richard W. Greenwood<br><a href="http://www.greenwoodmap.com" target="_blank">www.greenwoodmap.com</a></div></div>