Coordinate system for shapefile, no .prj file, using R












2















I am trying to use the shapefiles found on this page for oil & gas data from the Netherlands: https://www.nlog.nl/en/files-interactive-map from table at the bottom with "ARC_grid" and "Google Earth (WGS84)" columns and icons.



Specifically the "Oil and gas fields" ARC_grid files (https://www.nlog.nl/sites/default/files/nlog_velden_ed_1950_utm_31n_20170829.zip)



I am using R to download, read and plot the shapefile. I am able to do this OK, but the shapefile does not seem to have a coordinate reference system, I believe due to the absence of the .prj file.



The coordinates I get when I read the file in look like this:



 [1,] 683985.7 5931987
[2,] 684138.5 5931975


I would like them in Lat/Long, but do not know how to find out what conversion to use initially.



I am looking for an R-based solution.



Using R, this is how I have loaded the file:



temp <- tempfile(fileext = ".zip")
download.file("https://www.nlog.nl/sites/default/files/nlog_velden_ed_1950_utm_31n_20170829.zip", destfile = temp)
filepaths <- unzip(temp)

map <- readOGR(".")

map@proj4string
# CRS arguments: NA









share|improve this question



























    2















    I am trying to use the shapefiles found on this page for oil & gas data from the Netherlands: https://www.nlog.nl/en/files-interactive-map from table at the bottom with "ARC_grid" and "Google Earth (WGS84)" columns and icons.



    Specifically the "Oil and gas fields" ARC_grid files (https://www.nlog.nl/sites/default/files/nlog_velden_ed_1950_utm_31n_20170829.zip)



    I am using R to download, read and plot the shapefile. I am able to do this OK, but the shapefile does not seem to have a coordinate reference system, I believe due to the absence of the .prj file.



    The coordinates I get when I read the file in look like this:



     [1,] 683985.7 5931987
    [2,] 684138.5 5931975


    I would like them in Lat/Long, but do not know how to find out what conversion to use initially.



    I am looking for an R-based solution.



    Using R, this is how I have loaded the file:



    temp <- tempfile(fileext = ".zip")
    download.file("https://www.nlog.nl/sites/default/files/nlog_velden_ed_1950_utm_31n_20170829.zip", destfile = temp)
    filepaths <- unzip(temp)

    map <- readOGR(".")

    map@proj4string
    # CRS arguments: NA









    share|improve this question

























      2












      2








      2








      I am trying to use the shapefiles found on this page for oil & gas data from the Netherlands: https://www.nlog.nl/en/files-interactive-map from table at the bottom with "ARC_grid" and "Google Earth (WGS84)" columns and icons.



      Specifically the "Oil and gas fields" ARC_grid files (https://www.nlog.nl/sites/default/files/nlog_velden_ed_1950_utm_31n_20170829.zip)



      I am using R to download, read and plot the shapefile. I am able to do this OK, but the shapefile does not seem to have a coordinate reference system, I believe due to the absence of the .prj file.



      The coordinates I get when I read the file in look like this:



       [1,] 683985.7 5931987
      [2,] 684138.5 5931975


      I would like them in Lat/Long, but do not know how to find out what conversion to use initially.



      I am looking for an R-based solution.



      Using R, this is how I have loaded the file:



      temp <- tempfile(fileext = ".zip")
      download.file("https://www.nlog.nl/sites/default/files/nlog_velden_ed_1950_utm_31n_20170829.zip", destfile = temp)
      filepaths <- unzip(temp)

      map <- readOGR(".")

      map@proj4string
      # CRS arguments: NA









      share|improve this question














      I am trying to use the shapefiles found on this page for oil & gas data from the Netherlands: https://www.nlog.nl/en/files-interactive-map from table at the bottom with "ARC_grid" and "Google Earth (WGS84)" columns and icons.



      Specifically the "Oil and gas fields" ARC_grid files (https://www.nlog.nl/sites/default/files/nlog_velden_ed_1950_utm_31n_20170829.zip)



      I am using R to download, read and plot the shapefile. I am able to do this OK, but the shapefile does not seem to have a coordinate reference system, I believe due to the absence of the .prj file.



      The coordinates I get when I read the file in look like this:



       [1,] 683985.7 5931987
      [2,] 684138.5 5931975


      I would like them in Lat/Long, but do not know how to find out what conversion to use initially.



      I am looking for an R-based solution.



      Using R, this is how I have loaded the file:



      temp <- tempfile(fileext = ".zip")
      download.file("https://www.nlog.nl/sites/default/files/nlog_velden_ed_1950_utm_31n_20170829.zip", destfile = temp)
      filepaths <- unzip(temp)

      map <- readOGR(".")

      map@proj4string
      # CRS arguments: NA






      shapefile r






      share|improve this question













      share|improve this question











      share|improve this question




      share|improve this question










      asked Feb 8 at 19:43









      moman822moman822

      1153




      1153






















          3 Answers
          3






          active

          oldest

          votes


















          3














          Go to projfinder, type in one of your coordinate points and zoom the map to where you think the point should be:



          http://projfinder.com/



          projfinder will then try lots of coordinate systems and for the ones that map onto your map it will list them.



          enter image description here



          you then need to use your knowledge to figure out which one its most likely to be. In this case WGS84 UTM zone 31N unless its old data that might have used a previous geodetic reference like WGS 72 - but they're pretty close together.






          share|improve this answer































            3














            From a quick look, most other datasets on the website have a specific PRJ file - perhaps that particular one just wasn't copied over properly. I'd try copying one from a related dataset, put it in the same directory as the shapefile you need and rename it to: nlog_velden_ed_1950_utm_31n_20170829.prj.



            Then you can reproject as needed in ArcGIS. I'd double check to make sure everything seems right, but I'd bet that would do it.






            share|improve this answer

































              2














              The data appear to be projected in UTM Zone 31N, given the naming convention of the files and the coordinate values. That makes sense, as that UTM zone covers the western half of the Netherlands.






              share|improve this answer























                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
                });


                }
                });














                draft saved

                draft discarded


















                StackExchange.ready(
                function () {
                StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fgis.stackexchange.com%2fquestions%2f311563%2fcoordinate-system-for-shapefile-no-prj-file-using-r%23new-answer', 'question_page');
                }
                );

                Post as a guest















                Required, but never shown

























                3 Answers
                3






                active

                oldest

                votes








                3 Answers
                3






                active

                oldest

                votes









                active

                oldest

                votes






                active

                oldest

                votes









                3














                Go to projfinder, type in one of your coordinate points and zoom the map to where you think the point should be:



                http://projfinder.com/



                projfinder will then try lots of coordinate systems and for the ones that map onto your map it will list them.



                enter image description here



                you then need to use your knowledge to figure out which one its most likely to be. In this case WGS84 UTM zone 31N unless its old data that might have used a previous geodetic reference like WGS 72 - but they're pretty close together.






                share|improve this answer




























                  3














                  Go to projfinder, type in one of your coordinate points and zoom the map to where you think the point should be:



                  http://projfinder.com/



                  projfinder will then try lots of coordinate systems and for the ones that map onto your map it will list them.



                  enter image description here



                  you then need to use your knowledge to figure out which one its most likely to be. In this case WGS84 UTM zone 31N unless its old data that might have used a previous geodetic reference like WGS 72 - but they're pretty close together.






                  share|improve this answer


























                    3












                    3








                    3







                    Go to projfinder, type in one of your coordinate points and zoom the map to where you think the point should be:



                    http://projfinder.com/



                    projfinder will then try lots of coordinate systems and for the ones that map onto your map it will list them.



                    enter image description here



                    you then need to use your knowledge to figure out which one its most likely to be. In this case WGS84 UTM zone 31N unless its old data that might have used a previous geodetic reference like WGS 72 - but they're pretty close together.






                    share|improve this answer













                    Go to projfinder, type in one of your coordinate points and zoom the map to where you think the point should be:



                    http://projfinder.com/



                    projfinder will then try lots of coordinate systems and for the ones that map onto your map it will list them.



                    enter image description here



                    you then need to use your knowledge to figure out which one its most likely to be. In this case WGS84 UTM zone 31N unless its old data that might have used a previous geodetic reference like WGS 72 - but they're pretty close together.







                    share|improve this answer












                    share|improve this answer



                    share|improve this answer










                    answered Feb 8 at 20:59









                    SpacedmanSpacedman

                    24.1k23551




                    24.1k23551

























                        3














                        From a quick look, most other datasets on the website have a specific PRJ file - perhaps that particular one just wasn't copied over properly. I'd try copying one from a related dataset, put it in the same directory as the shapefile you need and rename it to: nlog_velden_ed_1950_utm_31n_20170829.prj.



                        Then you can reproject as needed in ArcGIS. I'd double check to make sure everything seems right, but I'd bet that would do it.






                        share|improve this answer






























                          3














                          From a quick look, most other datasets on the website have a specific PRJ file - perhaps that particular one just wasn't copied over properly. I'd try copying one from a related dataset, put it in the same directory as the shapefile you need and rename it to: nlog_velden_ed_1950_utm_31n_20170829.prj.



                          Then you can reproject as needed in ArcGIS. I'd double check to make sure everything seems right, but I'd bet that would do it.






                          share|improve this answer




























                            3












                            3








                            3







                            From a quick look, most other datasets on the website have a specific PRJ file - perhaps that particular one just wasn't copied over properly. I'd try copying one from a related dataset, put it in the same directory as the shapefile you need and rename it to: nlog_velden_ed_1950_utm_31n_20170829.prj.



                            Then you can reproject as needed in ArcGIS. I'd double check to make sure everything seems right, but I'd bet that would do it.






                            share|improve this answer















                            From a quick look, most other datasets on the website have a specific PRJ file - perhaps that particular one just wasn't copied over properly. I'd try copying one from a related dataset, put it in the same directory as the shapefile you need and rename it to: nlog_velden_ed_1950_utm_31n_20170829.prj.



                            Then you can reproject as needed in ArcGIS. I'd double check to make sure everything seems right, but I'd bet that would do it.







                            share|improve this answer














                            share|improve this answer



                            share|improve this answer








                            edited Feb 9 at 2:10









                            Vince

                            14.7k32749




                            14.7k32749










                            answered Feb 8 at 19:57









                            footloosietyfootloosiety

                            313




                            313























                                2














                                The data appear to be projected in UTM Zone 31N, given the naming convention of the files and the coordinate values. That makes sense, as that UTM zone covers the western half of the Netherlands.






                                share|improve this answer




























                                  2














                                  The data appear to be projected in UTM Zone 31N, given the naming convention of the files and the coordinate values. That makes sense, as that UTM zone covers the western half of the Netherlands.






                                  share|improve this answer


























                                    2












                                    2








                                    2







                                    The data appear to be projected in UTM Zone 31N, given the naming convention of the files and the coordinate values. That makes sense, as that UTM zone covers the western half of the Netherlands.






                                    share|improve this answer













                                    The data appear to be projected in UTM Zone 31N, given the naming convention of the files and the coordinate values. That makes sense, as that UTM zone covers the western half of the Netherlands.







                                    share|improve this answer












                                    share|improve this answer



                                    share|improve this answer










                                    answered Feb 8 at 19:57









                                    lambertjlambertj

                                    1,9031623




                                    1,9031623






























                                        draft saved

                                        draft discarded




















































                                        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.




                                        draft saved


                                        draft discarded














                                        StackExchange.ready(
                                        function () {
                                        StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fgis.stackexchange.com%2fquestions%2f311563%2fcoordinate-system-for-shapefile-no-prj-file-using-r%23new-answer', 'question_page');
                                        }
                                        );

                                        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







                                        Popular posts from this blog

                                        Probability when a professor distributes a quiz and homework assignment to a class of n students.

                                        Aardman Animations

                                        Are they similar matrix