I have a shapefile of the Bering Sea in the Pacific Ocean I've been working with in R. It was originally split down the middle, so I reprojected it to Mollweide to connect them (image example attached). The only problem is that there's still the boundary line down the middle. How would I get rid of it? They aren't two polygons that can be combined with Dissolve for some reason. Any advice would be greatly appreciatined. Here is my code so far:

Bering_Sea<- readOGR("Bering_Sea.shp")

Bering_Transform<- spTransform(Bering_Sea, CRS("+proj=moll +lon_0=150 +ellps=WGS84"))


enter image description here

enter image description here

  • I would try making a new feature class, then copying and pasting the two into it, while you're in the projection you like, then merging them via editor... – Michael Mar 21 '18 at 19:55
  • Why can't they be dissolved? – Simbamangu Mar 21 '18 at 20:15
  • What does dim or length tell you? If there are, in fact, two polygons, create a new column, in the @data slot, that contains identical values (eg. new.id = c(1,1) ) and then use rgeos::gUnaryUnion(Bering_Transform, id = Bering_Transform@data$new.id) – Jeffrey Evans Mar 21 '18 at 20:23

Browse other questions tagged or ask your own question.