R Merging Map ShapeFiles Together

I am attempting to merge two shapefiles together in R to plot.

I have downloaded the ShapeFile format of the political boundaries of the Carribbean from here and the ShapeFile of U.S. county borders from this link.

I have read in both files with

US_counties <- shapefile("~/us_county_shp_files/tl_2017_us_county.shp")

carribbean <-  shapefile("~/pol6bg/pol6bg.shp")

From looking at other posts, I have tried using both the union and the rbind methods.

With rbind (shown below), I get the error: identicalCRS(dots) is not TRUE.

rbind(US_counties, carribbean, makeUniqueIDs = TRUE)

With union, I get Error in as.vector(x) : no method for coercing this S4 class to a vector.

union(US_counties, carribbean)

Thanks in advance!

edited Nov 19 '18 at 17:53

asked Nov 19 '18 at 16:24

Alan

9110

add a comment |

I am attempting to merge two shapefiles together in R to plot.

I have downloaded the ShapeFile format of the political boundaries of the Carribbean from here and the ShapeFile of U.S. county borders from this link.

I have read in both files with

US_counties <- shapefile("~/us_county_shp_files/tl_2017_us_county.shp")

carribbean <-  shapefile("~/pol6bg/pol6bg.shp")

From looking at other posts, I have tried using both the union and the rbind methods.

With rbind (shown below), I get the error: identicalCRS(dots) is not TRUE.

rbind(US_counties, carribbean, makeUniqueIDs = TRUE)

With union, I get Error in as.vector(x) : no method for coercing this S4 class to a vector.

union(US_counties, carribbean)

Thanks in advance!

edited Nov 19 '18 at 17:53

asked Nov 19 '18 at 16:24

Alan

9110

add a comment |

I am attempting to merge two shapefiles together in R to plot.

I have downloaded the ShapeFile format of the political boundaries of the Carribbean from here and the ShapeFile of U.S. county borders from this link.

I have read in both files with

US_counties <- shapefile("~/us_county_shp_files/tl_2017_us_county.shp")

carribbean <-  shapefile("~/pol6bg/pol6bg.shp")

From looking at other posts, I have tried using both the union and the rbind methods.

With rbind (shown below), I get the error: identicalCRS(dots) is not TRUE.

rbind(US_counties, carribbean, makeUniqueIDs = TRUE)

With union, I get Error in as.vector(x) : no method for coercing this S4 class to a vector.

union(US_counties, carribbean)

Thanks in advance!

edited Nov 19 '18 at 17:53

asked Nov 19 '18 at 16:24

Alan

9110

I am attempting to merge two shapefiles together in R to plot.

I have downloaded the ShapeFile format of the political boundaries of the Carribbean from here and the ShapeFile of U.S. county borders from this link.

I have read in both files with

US_counties <- shapefile("~/us_county_shp_files/tl_2017_us_county.shp")

carribbean <-  shapefile("~/pol6bg/pol6bg.shp")

From looking at other posts, I have tried using both the union and the rbind methods.

With rbind (shown below), I get the error: identicalCRS(dots) is not TRUE.

rbind(US_counties, carribbean, makeUniqueIDs = TRUE)

With union, I get Error in as.vector(x) : no method for coercing this S4 class to a vector.

union(US_counties, carribbean)

Thanks in advance!

r shapefile

edited Nov 19 '18 at 17:53

asked Nov 19 '18 at 16:24

Alan

9110

edited Nov 19 '18 at 17:53

asked Nov 19 '18 at 16:24

Alan

9110

edited Nov 19 '18 at 17:53

asked Nov 19 '18 at 16:24

Alan

9110

asked Nov 19 '18 at 16:24

Alan

9110

asked Nov 19 '18 at 16:24

Alan

9110

add a comment |

1 Answer
1

active

oldest

votes

Firstly, the error when using union seems to arise due to conflicting packages. In particular, not the right function union is being called. Using

raster:::union(US_counties, carribbean)

gives

Warning message:

In union(US_counties, carribbean) : non identical CRS

which possibly is an important warning and says the same as when trying to use `rbind. Perhaps you know more about CRS (Coordinate Reference Systems) than me, but the problem is that

proj4string(carribbean)

# [1] "+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0"

proj4string(US_counties)

# [1] "+proj=longlat +datum=NAD83 +no_defs +ellps=GRS80 +towgs84=0,0,0"

So, not very surprisingly, the two shape files are of somewhat different formats. With the help of this useful read we see that it's not hard to unify those formats, after all. E.g.,

library(sp)

carribbean <- spTransform(carribbean, proj4string(US_counties))

Then

rbind(US_counties, carribbean, makeUniqueIDs = TRUE)

# Error in rbind(deparse.level, ...) : 

#   numbers of columns of arguments do not match

still doesn't work, but the problem is clear given that we are trying to bind different variables

names(US_counties)

#  [1] "STATEFP"  "COUNTYFP" "COUNTYNS" "GEOID"    "NAME"     "NAMELSAD" "LSAD"     "CLASSFP" 

#  [9] "MTFCC"    "CSAFP"    "CBSAFP"   "METDIVFP" "FUNCSTAT" "ALAND"    "AWATER"   "INTPTLAT"

# [17] "INTPTLON"

names(carribbean)

# [1] "F_CODE"     "NAM"        "NA2"        "NA2_DESCRI" "NA3"        "NA3_DESCRI" "WFBPATHLN" 

# [8] "AREA"       "PERIMETER"

However,

raster:::union(US_counties, carribbean)

now works and shows no warnings.

edited Nov 20 '18 at 2:21

answered Nov 19 '18 at 18:33

Julius Vainora

36.8k76481

Thanks for the answer Julius! When I try to replicate your code, I get this error from the union method call: Error in as.vector(x) : no method for coercing this S4 class to a vector. Is it because of a package I am missing?

– Alan
Nov 20 '18 at 1:51

1

@Alan, see the update. It was actually more about too many packages rather than missing something.

– Julius Vainora
Nov 20 '18 at 2:22

add a comment |

Your Answer

StackExchange.ifUsing("editor", function () {
StackExchange.using("externalEditor", function () {
StackExchange.using("snippets", function () {
StackExchange.snippets.init();
});
});
}, "code-snippets");

StackExchange.ready(function() {
var channelOptions = {
tags: "".split(" "),
id: "1"
};
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: true,
noModals: true,
showLowRepImageUploadWarning: true,
reputationToPostImages: 10,
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

Sign up or log in

StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});

Post as a guest

Name

Required, but never shown

StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fstackoverflow.com%2fquestions%2f53378845%2fr-merging-map-shapefiles-together%23new-answer', 'question_page');
}
);

Post as a guest

Name

Required, but never shown

1 Answer
1

active

oldest

votes

1 Answer
1

active

oldest

votes

Firstly, the error when using union seems to arise due to conflicting packages. In particular, not the right function union is being called. Using

raster:::union(US_counties, carribbean)

gives

Warning message:

In union(US_counties, carribbean) : non identical CRS

which possibly is an important warning and says the same as when trying to use `rbind. Perhaps you know more about CRS (Coordinate Reference Systems) than me, but the problem is that

proj4string(carribbean)

# [1] "+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0"

proj4string(US_counties)

# [1] "+proj=longlat +datum=NAD83 +no_defs +ellps=GRS80 +towgs84=0,0,0"

So, not very surprisingly, the two shape files are of somewhat different formats. With the help of this useful read we see that it's not hard to unify those formats, after all. E.g.,

library(sp)

carribbean <- spTransform(carribbean, proj4string(US_counties))

Then

rbind(US_counties, carribbean, makeUniqueIDs = TRUE)

# Error in rbind(deparse.level, ...) : 

#   numbers of columns of arguments do not match

still doesn't work, but the problem is clear given that we are trying to bind different variables

names(US_counties)

#  [1] "STATEFP"  "COUNTYFP" "COUNTYNS" "GEOID"    "NAME"     "NAMELSAD" "LSAD"     "CLASSFP" 

#  [9] "MTFCC"    "CSAFP"    "CBSAFP"   "METDIVFP" "FUNCSTAT" "ALAND"    "AWATER"   "INTPTLAT"

# [17] "INTPTLON"

names(carribbean)

# [1] "F_CODE"     "NAM"        "NA2"        "NA2_DESCRI" "NA3"        "NA3_DESCRI" "WFBPATHLN" 

# [8] "AREA"       "PERIMETER"

However,

raster:::union(US_counties, carribbean)

now works and shows no warnings.

edited Nov 20 '18 at 2:21

answered Nov 19 '18 at 18:33

Julius Vainora

36.8k76481

Thanks for the answer Julius! When I try to replicate your code, I get this error from the union method call: Error in as.vector(x) : no method for coercing this S4 class to a vector. Is it because of a package I am missing?

– Alan
Nov 20 '18 at 1:51

1

@Alan, see the update. It was actually more about too many packages rather than missing something.

– Julius Vainora
Nov 20 '18 at 2:22

add a comment |

Firstly, the error when using union seems to arise due to conflicting packages. In particular, not the right function union is being called. Using

raster:::union(US_counties, carribbean)

gives

Warning message:

In union(US_counties, carribbean) : non identical CRS

which possibly is an important warning and says the same as when trying to use `rbind. Perhaps you know more about CRS (Coordinate Reference Systems) than me, but the problem is that

proj4string(carribbean)

# [1] "+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0"

proj4string(US_counties)

# [1] "+proj=longlat +datum=NAD83 +no_defs +ellps=GRS80 +towgs84=0,0,0"

So, not very surprisingly, the two shape files are of somewhat different formats. With the help of this useful read we see that it's not hard to unify those formats, after all. E.g.,

library(sp)

carribbean <- spTransform(carribbean, proj4string(US_counties))

Then

rbind(US_counties, carribbean, makeUniqueIDs = TRUE)

# Error in rbind(deparse.level, ...) : 

#   numbers of columns of arguments do not match

still doesn't work, but the problem is clear given that we are trying to bind different variables

names(US_counties)

#  [1] "STATEFP"  "COUNTYFP" "COUNTYNS" "GEOID"    "NAME"     "NAMELSAD" "LSAD"     "CLASSFP" 

#  [9] "MTFCC"    "CSAFP"    "CBSAFP"   "METDIVFP" "FUNCSTAT" "ALAND"    "AWATER"   "INTPTLAT"

# [17] "INTPTLON"

names(carribbean)

# [1] "F_CODE"     "NAM"        "NA2"        "NA2_DESCRI" "NA3"        "NA3_DESCRI" "WFBPATHLN" 

# [8] "AREA"       "PERIMETER"

However,

raster:::union(US_counties, carribbean)

now works and shows no warnings.

edited Nov 20 '18 at 2:21

answered Nov 19 '18 at 18:33

Julius Vainora

36.8k76481

Thanks for the answer Julius! When I try to replicate your code, I get this error from the union method call: Error in as.vector(x) : no method for coercing this S4 class to a vector. Is it because of a package I am missing?

– Alan
Nov 20 '18 at 1:51

1

@Alan, see the update. It was actually more about too many packages rather than missing something.

– Julius Vainora
Nov 20 '18 at 2:22

add a comment |

Firstly, the error when using union seems to arise due to conflicting packages. In particular, not the right function union is being called. Using

raster:::union(US_counties, carribbean)

gives

Warning message:

In union(US_counties, carribbean) : non identical CRS

which possibly is an important warning and says the same as when trying to use `rbind. Perhaps you know more about CRS (Coordinate Reference Systems) than me, but the problem is that

proj4string(carribbean)

# [1] "+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0"

proj4string(US_counties)

# [1] "+proj=longlat +datum=NAD83 +no_defs +ellps=GRS80 +towgs84=0,0,0"

So, not very surprisingly, the two shape files are of somewhat different formats. With the help of this useful read we see that it's not hard to unify those formats, after all. E.g.,

library(sp)

carribbean <- spTransform(carribbean, proj4string(US_counties))

Then

rbind(US_counties, carribbean, makeUniqueIDs = TRUE)

# Error in rbind(deparse.level, ...) : 

#   numbers of columns of arguments do not match

still doesn't work, but the problem is clear given that we are trying to bind different variables

names(US_counties)

#  [1] "STATEFP"  "COUNTYFP" "COUNTYNS" "GEOID"    "NAME"     "NAMELSAD" "LSAD"     "CLASSFP" 

#  [9] "MTFCC"    "CSAFP"    "CBSAFP"   "METDIVFP" "FUNCSTAT" "ALAND"    "AWATER"   "INTPTLAT"

# [17] "INTPTLON"

names(carribbean)

# [1] "F_CODE"     "NAM"        "NA2"        "NA2_DESCRI" "NA3"        "NA3_DESCRI" "WFBPATHLN" 

# [8] "AREA"       "PERIMETER"

However,

raster:::union(US_counties, carribbean)

now works and shows no warnings.

edited Nov 20 '18 at 2:21

answered Nov 19 '18 at 18:33

Julius Vainora

36.8k76481

Firstly, the error when using union seems to arise due to conflicting packages. In particular, not the right function union is being called. Using

raster:::union(US_counties, carribbean)

gives

Warning message:

In union(US_counties, carribbean) : non identical CRS

which possibly is an important warning and says the same as when trying to use `rbind. Perhaps you know more about CRS (Coordinate Reference Systems) than me, but the problem is that

proj4string(carribbean)

# [1] "+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0"

proj4string(US_counties)

# [1] "+proj=longlat +datum=NAD83 +no_defs +ellps=GRS80 +towgs84=0,0,0"

So, not very surprisingly, the two shape files are of somewhat different formats. With the help of this useful read we see that it's not hard to unify those formats, after all. E.g.,

library(sp)

carribbean <- spTransform(carribbean, proj4string(US_counties))

Then

rbind(US_counties, carribbean, makeUniqueIDs = TRUE)

# Error in rbind(deparse.level, ...) : 

#   numbers of columns of arguments do not match

still doesn't work, but the problem is clear given that we are trying to bind different variables

names(US_counties)

#  [1] "STATEFP"  "COUNTYFP" "COUNTYNS" "GEOID"    "NAME"     "NAMELSAD" "LSAD"     "CLASSFP" 

#  [9] "MTFCC"    "CSAFP"    "CBSAFP"   "METDIVFP" "FUNCSTAT" "ALAND"    "AWATER"   "INTPTLAT"

# [17] "INTPTLON"

names(carribbean)

# [1] "F_CODE"     "NAM"        "NA2"        "NA2_DESCRI" "NA3"        "NA3_DESCRI" "WFBPATHLN" 

# [8] "AREA"       "PERIMETER"

However,

raster:::union(US_counties, carribbean)

now works and shows no warnings.

edited Nov 20 '18 at 2:21

answered Nov 19 '18 at 18:33

Julius Vainora

36.8k76481

edited Nov 20 '18 at 2:21

answered Nov 19 '18 at 18:33

Julius Vainora

36.8k76481

answered Nov 19 '18 at 18:33

Julius Vainora

36.8k76481

answered Nov 19 '18 at 18:33

Julius Vainora

36.8k76481

Thanks for the answer Julius! When I try to replicate your code, I get this error from the union method call: Error in as.vector(x) : no method for coercing this S4 class to a vector. Is it because of a package I am missing?

– Alan
Nov 20 '18 at 1:51

1

@Alan, see the update. It was actually more about too many packages rather than missing something.

– Julius Vainora
Nov 20 '18 at 2:22

add a comment |

Thanks for the answer Julius! When I try to replicate your code, I get this error from the union method call: Error in as.vector(x) : no method for coercing this S4 class to a vector. Is it because of a package I am missing?

– Alan
Nov 20 '18 at 1:51

1

@Alan, see the update. It was actually more about too many packages rather than missing something.

– Julius Vainora
Nov 20 '18 at 2:22

Thanks for the answer Julius! When I try to replicate your code, I get this error from the union method call: Error in as.vector(x) : no method for coercing this S4 class to a vector. Is it because of a package I am missing?

– Alan
Nov 20 '18 at 1:51

@Alan, see the update. It was actually more about too many packages rather than missing something.

– Julius Vainora
Nov 20 '18 at 2:22

add a comment |

draft saved

draft discarded

Thanks for contributing an answer to Stack Overflow!

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

Sign up or log in

StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});

Post as a guest

Name

Required, but never shown

Post as a guest

Name

Required, but never shown

Sign up or log in

StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});

Post as a guest

Name

Required, but never shown

Sign up or log in

StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});

Post as a guest

Name

Required, but never shown

Sign up or log in

StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});

Post as a guest

Name

Required, but never shown

Name

Required, but never shown

Name

Required, but never shown

This page is only for reference, If you need detailed information, please check here

搜尋此網誌

Agfdhyk