haodawei123的博客

原创 surf绘制表面图应用

ff=zeros(40,40,40);[x,y,z]=distancep(10,500,-150,-10,40,‘1.csv’);ff(:,:,10)=z;[x,y,z]=distancep(10,500,-150,-10,40,‘2.csv’);ff(:,:,20)=z;x=linspace(10,500,40);y=linspace(-150,-10,40);z=linspace(1,40,40);[xx,yy,zz]=meshgrid(x,y,z);hold on;%h=slice

原创 slice贴片图去网格

slice绘制的贴片图带网格，使用shading interp命令可以去掉网格不加shading interp 命令是这样的加了以后变成这样：ff=zeros(40,40,40);[x,y,z]=distancep(10,500,-150,-10,40,‘1.csv’);ff(:,:,10)=z;[x,y,z]=distancep(10,500,-150,-10,40,‘2.csv’);ff(:,:,20)=z;x=linspace(10,500,40);y=linspace(-150,

原创 matlab获取等值线的位置

使用连续的两行来存储坐标信息，如下图 等值线为1的点有5个，等值线为2的有13个，一直往下遍历1、做出图的话用，[X,Y] = meshgrid(-2:.2:2,-2:.2:3);Z = X.*exp(-X.2-Y.2);[C,h] = contour(X,Y,Z);C就是你需要的。或者（多此一举）C=get(h,‘ContourMatrix’) %一样的2、不作图，直接用一个输出得到等高矩阵[X,Y] = meshgrid(-2:.2:2,-2:.2:3);Z = X.*exp

原创 matlab生成exe

一、生成独立可执行的程序（exe文件）步骤需要注意一下目录，如我的m文件在MatFor3D目录下面，则需要双击进入该目录后再执行相关的命令。1、安装编译器。可有多种选择，matlab自带了一个LCC，推荐使用VC++6.0，我基于VS2010实现。2、设置编译器。在matlab命令行输入mbuild –setup以及mex –setup，选择安装的c编译器（可能会选1也可能选0，这个要看你显示的是什么）。选择y 一直到done3、调用编译器。当然，也可以输入mcc -m filaname，

原创 matlab 逻辑操作函数

原创 matlab listbox使用

1、几个重要的属性Tag属性设置控件的名称，m文件编辑获取控件时使用Value属性指示当前选中的是第几行String属性设置显示在listbox中的内容，可以在设计界面通过“ENTER”换行输入多行内容获取listbox内容function listbox2_Callback(hObject, eventdata, handles)% hObject handle to listbox2 (see GCBO)% eventdata reserved - to be define

原创 Matlab显示打开文件对话框

1、matlab7 输入guide命令打开对话框，插入一个按钮2、找到按钮的回调程序添加下面的内容使用了 uigetfile这个函数，这个函数会打开打开文件对话框，返回文件名称、路径、index索引% — Executes on button press in pushbutton4.function pushbutton4_Callback(hObject, eventdata, handles)% hObject handle to pushbutton4 (see GCBO)% ev

原创 matlab按钮不能显示中文处理

matlab7 按钮的string属性设置为中文不能显示，双击按钮打开属性面板，找到FontName，修改为“宋体”或者是其它的中文字体即可。

原创 matlab绘制平面等值线图

%首先使用surfer将高密度数据网格化，注意网格化的行与列要一致%即沿测线方向与往地下方向要分割相等，生成grd之后，将grd的数据%另存为csv,这里剖分为50x50a=load(‘1.csv’);x=linspace(7.5,487.5,50);y=linspace(5,165,50);z=a(:,3,:);%获取视电阻率值zz=reshape(z,50,50);%将z变形为50*50的矩阵zzz=zz’;%需要再取一个转置的操作[X,Y]=meshgrid(x,y);contou

原创 matlab实现3维测井曲线绘制

a=load(‘cj.csv’);z=a(:,1,:);%将第一列给zx=a(:,2,:);%将第二列给xy=a(:,3,:);%将第三列给yplot3(x,y,z,‘r’);set(gca,‘ZDir’,‘reverse’);%设置Z轴倒序box on;%显示盒子grid on;%开启网格xlim([-5,5]);xlabel(‘东西方向’);ylim([-5,5]);ylabel(‘南边方向’);line([0 0 0],[0 0 0],[0,0,180]);set(gcf,‘

转载 MapObejcts组件应用设计(地图数据的创建与编辑)

https://blog.csdn.net/lgmawei/article/details/1285718

转载 mapobject绘图

https://www.cnblogs.com/gispeng/articles/755823.html

原创 mapobject学习

1、下载安装mapobject2，并破解2、实例及代码2.1添加控件在窗体上添加两个map控件，一个legend控件，一个commondialog控件2.2在map1控件显示地图（*.shp）Private Sub Command1_Click()Dim dc As New DataConnectionDim gs As GeoDatasetDim name As StringDim layer As MapObjects2.MapLayerCommonDialog1.Filter =

原创 VC++在内存中绘图

HDC hdc=::GetDC(m_hWnd);RECT rect;GetClientRect(hwnd,&rect);int cxClient = rect.right;int cyClient = rect.bottom;//为显示设备创建一块相应的内存HDC hdcmem = CreateCompatibleDC(NULL);//创建位图HBITMAP hBitmap = CreateCompatibleBitmap(hdc,cxClient,cyClient);//将内存句

转载 python绘制三维曲线（转）

-- coding: UTF-8 --import matplotlib as mplfrom mpl_toolkits.mplot3d import Axes3Dimport numpy as npimport matplotlib.pyplot as plt设置图例字号mpl.rcParams[‘legend.fontsize’] = 10fig = plt.figure()设置三维图形模式ax = fig.gca(projection=‘3d’)测试数据theta = np.li

原创 pandas的两种定位方式loc与iloc

复制代码import pandas as pddf1 = pd.DataFrame([[‘Snow’,‘M’,22],[‘Tyrion’,‘M’,32],[‘Sansa’,‘F’,18],[‘Arya’,‘F’,14]], columns=[‘name’,‘gender’,‘age’])print("--------更换单个值----------")loc和iloc 可以更换单行、单列、多行、多列的值df1.loc[0,‘age’]=25 # 思路：先用loc找到要更改的值，再用赋值（=

原创 python 使用代理下载图片

使用urllib.request.urlretrieve方式：from urllib.request import urlretrieveurlretrieve(“https://www.awaimai.com/wp-content/uploads/2017/09/phpinfo.png”, “phpinfo.png”)如果要使用代理，需要先导入：socket和httplib2.socks模块：from urllib.request import urlretrieveimport config

原创 pandas操作技巧

创建DataFrame,转置操作由字典创建一个DateFramed={}#创建一个空的字典d[‘000001’]=[1.2,1.3,-1.3,1.4]##字典里面添加键值对d[‘000002’]=[2.2,2.3,-2.3,1.4]d[‘000003’]=[3.2,3.3,-3.3,3.4]df=pd.DataFrame(d)##创建DataFramedf2=pd.DataFrame(df.values.T, index=df.columns, columns=df.index)##计算转置

原创 Python-Pandas 打乱数据（shuffle）

在Python里面，使用Pandas里面的DataFrame来存放数据的时候想要把数据集进行shuffle会许多的方法，本文介绍两种比较常用而且简单的方法。应用情景：我们有下面以个DataFrame处理完之后的效果如下：实现方法：最简单的方法就是采用pandas中自带的 sample这个方法。假设df是这个DataFramedf.sample(frac=1)这样对可以对df进行shuffle。其中参数frac是要返回的比例，比如df中有10行数据，我只想返回其中的30%,那么frac=0

原创 tushare复权数据下载设置

首先打开下面网站登陆tushare，然后注册登陆https://tushare.pro/register?reg=403626切换到下面界面获取token在idle中输入一下代码即可：token=‘你复制下来的token’ts.set_token(token)pro=ts.pro_api()

原创 python list列表特殊使用

下面的一个二维列表tt[0:5]array([[ 34.19 , 34.89 , 34.81 , 34.01 , 23601.24 , 34.81 ,34.81 , 34.81 , 23601.24 , 23601.24 , 23601.24 ],[ 34.82 , 34.82 , 34.59 , 34.37 , 12093. , 34.7 ,34.7 , 34.7 , 17847.12 , 17847.12 ,

原创 py2exe总结

将py脚本编译成exe需要py2exe、py脚本、setup.py脚本（设置编译参数）1 安装py2execmd下输入：pip install py2exe注意通过我的测试发现3.6、3.7、3.8版本的py2exe转换py到exe的时候都会出错，只有3.4版本的没有问题，这里我用的是x86的3.4版测试的2转换py到exe2.1 setup.py部分#mysetup.pyfrom distutils.core import setupimport py2exesetup(console=

原创 vc调用python函数返回列表的处理

1、python代码输入两个变量，返回一个列表xdef main(a,b): with open("1.txt","w") as fw: fw.write(a) fw.write(b) x=[1,2,3,4] return x,len(x) #//return x,len(x)#两个值就是元组了，也可以处理2、vc代码 Py_Initialize(); PyObject * pModule = NULL; PyObject * p

原创 vc调用python返回值处理

1、c代码*通过测试发现只有 char a="";这样定义的字符型指针能通过PyArg_ParseTuple返回正确的字符串，使用char s1[100]={0};这样的形式定义的字符串将s1传入PyArg_ParseTuple函数会返回错误的字符串。Py_Initialize();PyObject * pModule = NULL;PyObject * pFunc = NULL;PyObject * ret = NULL;pModule = PyImport_ImportModule("1")

原创 vc调用python函数传递可变参数

在绝大多数的文章中vc调用python的函数传递的参数都是确定的值，通过Py_BuildValue(“s”,“12334”)，这样的形式实现的，有时候我们想传递一个可变的参数如有个vc的界面输入一段内容，将这段内容通过python的函数显示出来。下面是完整代码：// test.cpp : 定义控制台应用程序的入口点。//#include "stdafx.h"#include <Python.h>#include <windows.h>#include <stri

原创 北英武测水

定了119号点、189号点、245号点3处位置，打在189号点，20多米见破碎带水量大，估计245号点会更好。

原创 python np.linspace函数使用

linspace(start, stop, num=50, endpoint=True, retstep=False, dtype=None)x3=np.linspace(1,5,5,dtype=int) #我们注意观察会发现，默认输出的是浮点形，我们这里改为整形。print(x3)

原创 python threading.Thread实现多线程

参考下面的网址https://www.cnblogs.com/csnd/p/11807840.html#%E4%BD%BF%E7%94%A8%C2%A0threading%E6%A8%A1%E5%9D%97%E5%88%9B%E5%BB%BA%E7%BA%BF%E7%A8%8B自己实现了一下运行效果如下import threadingimport timefrom datetime import datetimedef Test(name): for i in range(3):

原创 python response.iter_content()的使用

使用requests模块下载文件时使用流氏下载，或者要查看文件的下载进度的时候会用到 response.iter_content()def gettsfile(i):#传入待下载的ts文件序号，找到对应的下载地址下载 url="http://y.downya.com/soft/MATLABchengxuyuyan_downyi.com.zip" user_agent="Mozilla/5.0 (Windows NT 6.1; WOW64) AppleWebKit/537.36 (KHTML

原创 python求解线性方程组的解

>>>import numpy as np>>> a=[[4,6,2],[3,4,1],[2,8,13]]>>> s=[9,7,2]>>> x=np.linalg.solve(a,s)>>> xarray([ 3. , -0.5, 0. ])

原创 第二类修正贝塞尔函数

1 递推算式摘自《MATLAB语音常用程序集》第二类修正贝塞尔函数需要用第一类修正贝塞尔函数来计算2 vc编程实现#include <stdio.h>#include <math.h>double first_xzbsl(int n,double x)//第一类修正的贝塞尔函数{ double a[7]={1.0,3.5156229,3.0899424,1.2067492,0.2659732,0.0360768,0.0045813}; double c[7]

原创 python中的复数操作（实部、虚部、求模、求共轭）

Python 3.6.2 (v3.6.2:5fd33b5, Jul 8 2017, 04:57:36) [MSC v.1900 64 bit (AMD64)] on win32Type "copyright", "credits" or "license()" for more information.>>> aa=123-12j>>> aa.real #（取实部）123.0>>> aa.imag #（取虚部）-12.0>>&

原创 栖霞杨础找水总结

栖霞杨础测水，本次工作是在下雨之后的第二天立马去做的。从下图上看，3个图的对应都很好，并且3号异常位于已知井位的正北小于20m的位置，所以认为本次测量的是可靠的。一开始根据下图圈定了1，2号异常，先打的1号位置打到150m，在80m附近有一层水0.5方左右，基本是废井。去现场确定1号异常位置与实际打井位置正确，后又布设2号异常位置，打了150m没水。第三次定在3号异常位置，地表浅部就有水，大水在120m深度，最终是10方水。从MNB图上看，OB在120m-180m等值线有个区域数值不变。120mX0

原创 dx 生成dex文件实现

1、新建一个java项目并编译新开一个cmd，javac hello.java 生成hello.class文件2、使用android sdk build-tools 的dx工具生成dex文件将android sdk build-tools目录添加到用户环境变量，然后到项目的src目录新开一个cmd窗口执行下面的命令即可在src目录下面生成hello.dex文件dx --dex --output=hello.dex wei\hello.class...

原创 eclipse JNI 开发总结

经过将近两天的努力终于在eclipse上实现了JNI编程，这里记录一下！！！！1、JNI开发的组成两部分一是java或android部分（函数方法的声明），二是so或者dll的实现部分（函数方法的具体实现）2、javah 生成.h头文件将 javac 将java转换成class，然后利用javah命令生成 so 或者 dll调用的 .h头文件这里容易出错的是生成.h头文件会出现找不到类文件的错误，参见上一篇文章解决文章地址就是使用 javah -classpath 指定src的目录，多了这么一步

原创 javah -jni 找不到类的解决

被这个问题困扰了1天1、建立JNI项目eclipse新建一个项目JNIDemo,添加一个wei的包，创建一个HelloWorld的类添加如下内容：添加一个native的sayHello方法package wei;public class HelloWorld { public native void sayHello(); public static void main(String[] arsc) { System.out.print("Hello!!!!!!!"); System

原创 android studio获取IMEI码

获取到的就是*#06#的码，public class MainActivity extends AppCompatActivity { public String GetPhoneNumber(Context paramContext)//#手机IMEI码号码 { TelephonyManager paramTelephonyManager = (TelephonyManager) paramContext.getSystemService(Context.TELEPHON

原创 快速傅里叶变换及其逆变换的python的实现

代码实现了快速傅里叶变换及其逆变换的计算，幅值、初始相位的计算import numpy as npimport matplotlib.pyplot as pltFs=500#采样频率ts=1/Fs#采样时间间隔L=500#信号长度t=np.linspace(0,L-1,L)*tsx=2+0.7*np.cos(2*np.pi*50*t+np.pi/2)+2*np.cos(2*np.pi*200*t+np.pi/6)y=x#+2*np.random.rand()fft_y=np.fft.ff

转载 Android反调试方法总结以及源码实现之检测篇（一）（转）

https://blog.csdn.net/feibabeibei_beibei/article/details/60956307

原创 windows脱壳总结

1、单步跟踪法原理：单步跟踪就是有向上的跳转就在下面的代码上面F4，执行的跳转下面的代码，一直找直到找到OEP2、ESP定律对UPX ASPack的壳查找OEP使用标志是第一行代码是 pushad单步pushad然后在ESP寄存器上右键数据窗口中跟随数据窗口中ESP寄存器的地址处下硬件访问断点，F9让程序跑起来断下之后JMP xxx 基本就是ESP了3、2次镜像法od中ALT+M打开内存，在程序的第一个.rsrc区段上F2断点，F9运行，断下之后再ALT+M找到00401000这个位置下