Run R script from Python -

-

i want run r-script python script. r-script required projection of lat lon coordinates in different coordinate system. have considered 2 options this. in first option parse lat , lon coordinates r-script shown below. r-script returns x , y python script, can't figure out how this.

project<-function(lat,lon){  library(sp) library(rgdal)  xy <- cbind(x = lon, y = lat) s <- spatialpoints(xy) proj4string(s) <- crs("+proj=longlat +ellps=wgs84 +datum=wgs84 +no_defs") snew <- sptransform(s, crs("+proj=sterea +lat_0=52.15616055555555 +lon_0=5.38763888888889 +k=0.9999079 +x_0=155000 +y_0=463000 +ellps=bessel +units=m +no_defs")) x <- coordinates(snew)[1] y <- coordinates(snew)[2]  return(x,y) }

for second option i've considered using r-script @ bottom lat , lon in it. try run python subprocess.popen('rscript project.r', shell=true).wait() not seem work. not write xy.txt file. if run cmd line, however, r-script job. can me out 1 of these 2 options?

library(sp) library(rgdal)  lat <- 52.29999924 lon <- 4.76999998  xy <- cbind(x = lon, y = lat) s <- spatialpoints(xy) proj4string(s) <- crs("+proj=longlat +ellps=wgs84 +datum=wgs84 +no_defs") snew <- sptransform(s, crs("+proj=sterea +lat_0=52.15616055555555 +lon_0=5.38763888888889 +k=0.9999079 +x_0=155000 +y_0=463000 +ellps=bessel +units=m +no_defs")) x <- coordinates(snew)[1] y <- coordinates(snew)[2]  cat(x, file="xy.txt",sep="") cat(y,file="xy.txt",append=true)

for first solution -- r output coordinates console (using print or cat or whatever in r), capture in python (see here: running shell command python , capturing output)

that give string teh lat/lon coords in it, have parse them out using appropraite python functions.


Comments

Post a Comment

Popular posts from this blog

java - Jmockit String final length method mocking Issue -

-
using jmockit 1.2 jar, trying mock string's length method getting unexpected invocation exception: failed: test java.lang.illegalstateexception: missing invocation mocked type @ point; please make sure such invocations appear after declaration of suitable mock field or parameter @ stringdemo.testa$1.<init>(testa.java:17) @ stringdemo.testa.test(testa.java:13) i using testng: @test public void test() throws exception { new expectations() { @mocked("length") string astring; { astring.length(); result = 2; } }; system.out.println(a.showa("test")); } } actual class a: public class { public static int showa(string str){ int a= str.length(); return a; } } this wrong way of recording expected results. shouldn't mock , record string's length() method, record showa() instead. here solution @suppresswarnings("unused&qu
Read more

asp.net - Razor Page Hosted on IIS 6 Fails Every Morning -

-
i have simple page in intranet uses razor/asp fetch single record table , display it, plus few graphic, dashboard display. meaning, no user intervention other first time open ie internal url. the problem having every morning ie displays "page not found or network error" message. refresh page, , it's stuck , doesn't display anything. try different pc, open ie internal url , it's stuck... ...until following: login server hosting page run inetmgr go web sites, etc right click on page that's giving me problem , select browse at moment, error message: server error in '/' application. -------------------------------------------------------------------------------- type of page not served. description: type of page have requested not served because has been explicitly forbidden. extension '.cshtml' may incorrect. please review url below , make sure spelled correctly. requested url: /application/dashboard.cshtml ---------------
Read more

c++ - wxwidget compiling on windows command prompt -

-
can please guide me on method,as how compile wxwidget using windows command prompt or msys. [i can't seem find on wxwiki/wxwidget official book/anywhere else] i have compiled wxwidget instructions have provided (using msys). from searching on internet seems 1 can (maybe?) through use of makefile (is correct?) if yes: directory & how should link wx libraries in makefile if no: way besides makefile? i know can use ide's code::block , make life simpler prefer compile using command prompt/msys. thanks in advance.
Read more