def loadPainting3Data(imageFile): fileName = "/Users/cengelsm/Test/pgsinterp3.p3" fis = FileInputStream(fileName) ois = PythonObjectInputStream(fis) p3 = ois.readObject() fis.close() image = readImage(imageFile) n1, n2, n3 = s1.count, s2.count, s3.count d1, d2, d3 = s1.delta, s2.delta, s3.delta f1, f2, f3 = s1.first, s2.first, s3.first l1, l2, l3 = s1.last, s2.last, s3.last ipg = ImagePanelGroup2(s1, s2, s3, image, p3.getPaint()) ipg.setColorModel1(ColorMap.getGray()) ipg.setColorModel2(ColorMap.getJet(0.5)) ipg.setSlices(700, 50, 80) sf = SimpleFrame() world = sf.getWorld() world.addChild(ipg) sf.setSize(1250, 900) view = sf.getOrbitView() zscale = 0.75 * max(n2 * d2, n3 * d3) / (n1 * d1) view.setAxesScale(1.0, 1.0, zscale) view.setScale(1.3) view.setAzimuth(-65.0) view.setWorldSphere(BoundingSphere(BoundingBox(f3, f2, f1, l3, l2, l1)))
def readTensors(): tensorsFile = "st"+str(int(sigmaTensor)) fis = FileInputStream(dataDir+dataPref+tensorsFile+".dat") ois = PythonObjectInputStream(fis) tensors = ois.readObject() fis.close() return tensors
def loadPainting3Data(imageFile): fileName = "/Users/cengelsm/Test/pgsinterp3.p3" fis = FileInputStream(fileName) ois = PythonObjectInputStream(fis) p3 = ois.readObject() fis.close() image = readImage(imageFile) n1,n2,n3 = s1.count,s2.count,s3.count d1,d2,d3 = s1.delta,s2.delta,s3.delta f1,f2,f3 = s1.first,s2.first,s3.first l1,l2,l3 = s1.last,s2.last,s3.last ipg = ImagePanelGroup2(s1,s2,s3,image,p3.getPaint()) ipg.setColorModel1(ColorMap.getGray()) ipg.setColorModel2(ColorMap.getJet(0.5)) ipg.setSlices(700,50,80) sf = SimpleFrame() world = sf.getWorld() world.addChild(ipg) sf.setSize(1250,900) view = sf.getOrbitView() zscale = 0.75*max(n2*d2,n3*d3)/(n1*d1) view.setAxesScale(1.0,1.0,zscale) view.setScale(1.3) view.setAzimuth(-65.0) view.setWorldSphere(BoundingSphere(BoundingBox(f3,f2,f1,l3,l2,l1)))
def readTensors(name): """ Reads tensors from file with specified basename; e.g., "tpet". """ fis = FileInputStream(seismicDir+name+".dat") ois = PythonObjectInputStream(fis) tensors = ois.readObject() fis.close() return tensors
def readTensors(name): """ Reads tensors from file with specified basename; e.g., "tpet". """ fis = FileInputStream(PATH+name) ois = PythonObjectInputStream(fis) tensors = ois.readObject() fis.close() return tensors
def readTensors(name): """ Reads tensors from file with specified basename; e.g., "f3et". """ fis = FileInputStream(f3dDataSetDir+name+".dat") ois = PythonObjectInputStream(fis) tensors = ois.readObject() fis.close() return tensors
def loadPaintingSlice(imageFile): fileName = "/Users/cengelsm/Test/paintings/pgs_02.p3" fis = FileInputStream(fileName) ois = PythonObjectInputStream(fis) p3 = ois.readObject() fis.close() image = readImage(imageFile) print s1.count, s2.count sl = 30 xz = image[sl] pxz = p3.getPaint()[sl] panel = PlotPanel(PlotPanel.Orientation.X1DOWN_X2RIGHT) pv1 = panel.addPixels(s1, s2, xz) pv1.setPercentiles(1, 99) pv2 = panel.addPixels(s1, s2, pxz) pv2.setColorModel(ColorMap.getJet(0.5)) panel.setVLabel("Depth (km)") panel.setHLabel("Distance (km)") frame = PlotFrame(panel) frame.setDefaultCloseOperation(PlotFrame.EXIT_ON_CLOSE) frame.setSize(1800, 1000) frame.setVisible(True) frame.setFontSizeForSlide(1, 1) frame.paintToPng(300, 6, "/Users/cengelsm/Desktop/test2.png")
def loadPaintingSlice(imageFile): fileName = "/Users/cengelsm/Test/paintings/pgs_02.p3" fis = FileInputStream(fileName) ois = PythonObjectInputStream(fis) p3 = ois.readObject() fis.close() image = readImage(imageFile) print s1.count,s2.count sl = 30 xz = image[sl] pxz = p3.getPaint()[sl] panel = PlotPanel(PlotPanel.Orientation.X1DOWN_X2RIGHT) pv1 = panel.addPixels(s1,s2,xz) pv1.setPercentiles(1,99); pv2 = panel.addPixels(s1,s2,pxz) pv2.setColorModel(ColorMap.getJet(0.5)) panel.setVLabel("Depth (km)") panel.setHLabel("Distance (km)") frame = PlotFrame(panel) frame.setDefaultCloseOperation(PlotFrame.EXIT_ON_CLOSE) frame.setSize(1800,1000) frame.setVisible(True) frame.setFontSizeForSlide(1,1) frame.paintToPng(300,6,"/Users/cengelsm/Desktop/test2.png")
def readObject(name): fis = FileInputStream(seismicDir+name+".dat") ois = PythonObjectInputStream(fis) obj = ois.readObject() ois.close() return obj
def readObject(name): fis = FileInputStream(seismicDir + name + ".dat") ois = PythonObjectInputStream(fis) obj = ois.readObject() ois.close() return obj
def readTensors(tensorsFile): fis = FileInputStream(tensorsFile) ois = PythonObjectInputStream(fis) tensors = ois.readObject() fis.close() return tensors
def readTensors(name): fis = FileInputStream(seismicDir+name+".dat") ois = PythonObjectInputStream(fis) tensors = ois.readObject() fis.close() return tensors
def readTensors(name): fis = FileInputStream(seismicDir + name + ".dat") ois = PythonObjectInputStream(fis) tensors = ois.readObject() fis.close() return tensors