編集履歴

質問編集履歴

分かっていない部分や怪しい部分を更新したのと太字にしました。

2019/11/02 04:31

投稿

littlestream

スコア37

test CHANGED Viewed

File without changes

test CHANGED Viewed

@@ -10,7 +10,7 @@
-この共有する頂点という意味がよく分からないです。
+**この共有する頂点という意味がよく分からないです。**
@@ -114,6 +114,14 @@
 (0.394055203,-0.472866243,0.788110406)
+光源ベクトルと頂点ベクトルの内積は、
+頂点ベクトルが(Nx,Ny,Nz),光源ベクトルが(Lx,Ly,Lz)だとすると、
+double DotProductResult=Nx*Lx+Ny*Ly+Nz*Lz;
 ちなみに、内積をこの法線ベクトルと視線ベクトル（カメラ？）と行う事で
@@ -126,9 +134,9 @@
-3.あと内積で求めた頂点カラーをどういうループにすれば
+**3.あと内積で求めた頂点カラーをどういうループにすれば
-頂点カラーを適用できるのか分からないです。
+頂点カラーを適用できるのか分からないです。**
@@ -146,6 +154,18 @@
+**なので、よく分かっていないのですが、面ごとの法線ベクトルを
+頂点ベクトルという変数をでっちあげて面ごとの法線ベクトルをその面ごとの頂点ベクトルに
+足すのと足した回数を調べておいて平均化したものに、ＲＧＢ値のＲ値、Ｇ値、Ｂ値
+それぞれに掛けてしまう方法なのかどうかすら分かってません。**
 ```c/c++言語

ソースの最新版です。

2019/11/02 04:31

投稿

littlestream

スコア37

test CHANGED Viewed

File without changes

test CHANGED Viewed

@@ -146,7 +146,7 @@
-```C言語
+```c/c++言語
@@ -592,10 +592,48 @@
+		//PSET(320,240,0x00ff0000);//点を打つ
+        //LINE(320,240,640,380,0x00ff0000);
+		int xy[4][2]={
+			320,5,
+			220,300,
+			500,430,
+			630,280
+		};
+		int Col[4]=
+		{
+			0x00ff0000,
+			0x0000ff00,
+			0x000000ff,
+			0x00ffffff
+		};
 		Projection2DAxisFrom3DAxis(hdc);//透視投影
+		//PSET(320,240,0x00FFFFFF);
 		XYZ vvv[12];//全ての頂点の法線をクリアする
@@ -610,12 +648,16 @@
 			vvv[i].z=0.0f;//Clear
-			//これが無いとどんどん加算されてしまう
 		}
+		XYZ a[13];
+        XYZ b[13];
+        XYZ n[13];
 		for(int i=0;i<12;i++)
@@ -690,102 +732,134 @@
-		//頂点ベクトルは求まった(はず...)から、頂点カラーを面にある全ての頂点で求める
+			//頂点ベクトルは求まった(はず...)から、頂点カラーを面にある全ての頂点で求める
 		double ViewCos[12];
+		//XYZ vpx=200,vpy=-120,vpz=100;
+		for(int i=0;i<12;i++)
+		{
+			double a=sqrt(vvv[i].x*vvv[i].x+vvv[i].y*vvv[i].y+vvv[i].z*vvv[i].z)+0.01f;
+			int i1=Poly[i][0];
+			int i2=Poly[i][1];
+			int i3=Poly[i][2];
+			XYZ n=CalcNormal2(Pnt[i1],Pnt[i2],Pnt[i3]);
+			vvv[i].x/=a;
+			vvv[i].y/=a;
+			vvv[i].z/=a;
+			XYZ light={100,-120,200};
+			double lenlight=sqrt((double)light.x*light.x+light.y*light.y+light.z*light.z);
+			light.x/=lenlight;
+			light.y/=lenlight;
+			light.z/=lenlight;
+			double spe1=vvv[i1].x*light.x+vvv[i1].y*light.y+vvv[i1].z*light.z;
+			double spe2=vvv[i2].x*light.x+vvv[i2].y*light.y+vvv[i2].z*light.z;
+			double spe3=vvv[i3].x*light.x+vvv[i3].y*light.y+vvv[i3].z*light.z;
+            double vv=sqrt((vpx*vpx+vpy*vpy+vpz*vpz));//視線ベクトルの長さを求める
+			double v1=vpx/(vv+0.01f);//視線ベクトルの正規化
+			double v2=vpy/(vv+0.01f);//
+			double v3=vpz/(vv+0.01f);//
+			ViewCos[i]=(v1*n.x+v2*n.y+v3*n.z);
+			//視線ベクトルと面ごとの法線ベクトルの内積
+			if(ViewCos[i]<0)
+				TRIANGLE2(Pnt2[Poly[i][0]].x,Pnt2[Poly[i][0]].y,
+					Pnt2[Poly[i][1]].x,Pnt2[Poly[i][1]].y,
+					Pnt2[Poly[i][2]].x,Pnt2[Poly[i][2]].y,
+					RGB(spe1*255,spe1*255,spe1*255),
+					RGB(spe2*255,spe2*255,spe2*255),
+					RGB(spe3*255,spe3*255,spe3*255)
+				);
+		}
+        //反時計回りの座標軸
+        //三角形を表示する
+			//DrawTextureMapping(320,0,0,480,640,0);
+			//FillGraRect(0,0,320,20,255,0,0,0,255,0);
+			CalcFPS(hdc);
+		/* DIBSectionをDIBとして描画 */
-		int vpx=200,vpy=-120,vpz=100;
+		StretchDIBits(hdc, 0, 0, WIDTH, HEIGHT,
+		     0, 0, WIDTH, HEIGHT, lpPixel, &biInfo,
+		     DIB_RGB_COLORS,SRCCOPY);
-		for(int i=0;i<12;i++)
+		EndPaint(hwnd, &ps);
-		{
-		double a=sqrt(vvv[i].x*vvv[i].x+vvv[i].y*vvv[i].y+vvv[i].z*vvv[i].z)+0.01f;
-			vvv[i].x/=a;
-			vvv[i].y/=a;
-			vvv[i].z/=a;
-            		double vv=sqrt((vpx*vpx+vpy*vpy+vpz*vpz));//視線ベクトルの長さを求める
+		DeleteObject(SelectObject(hdc , GetStockObject(WHITE_BRUSH)));
-			double v1=vpx/(vv+0.01f);//視線ベクトルの正規化
-			double v2=vpy/(vv+0.01f);//
-			double v3=vpz/(vv+0.01f);//
-			ViewCos[i]=(v1*vvv[i].x+v2*vvv[i].y+v3*vvv[i].z);
-			//視線ベクトルと面ごとの法線ベクトルの内積
-			if(ViewCos[i]<0)
-				TRIANGLE2(Pnt2[Poly[i][0]].x,Pnt2[Poly[i][0]].y,
-					Pnt2[Poly[i][1]].x,Pnt2[Poly[i][1]].y,
-					Pnt2[Poly[i][2]].x,Pnt2[Poly[i][2]].y,
-					RGB(255,255,255),
-					RGB(255,255,255),
-					RGB(255,255,255)
-				);
-			//とりあえず真っ白にしたが、頂点カラーを求めた後にどうすれば
-			//３点の色を指定する処理になるか分からないです。
 		}
-		CalcFPS(hdc);
-		/* DIBSectionをDIBとして描画 */
-		StretchDIBits(hdc, 0, 0, WIDTH, HEIGHT,
-		     0, 0, WIDTH, HEIGHT, lpPixel, &biInfo,
-		     DIB_RGB_COLORS,SRCCOPY);
-		EndPaint(hwnd, &ps);
-		DeleteObject(SelectObject(hdc , GetStockObject(WHITE_BRUSH)));
-		}
 		return 0;
@@ -823,3 +897,5 @@
 }
 ```
+![イメージ説明](6af0066a11228dbaa3ab7964da8a4eb2.jpeg)

インデントが上手く行かなかったので<code>を使います。

2019/11/01 12:42

投稿

littlestream

スコア37

test CHANGED Viewed

File without changes

test CHANGED Viewed

@@ -146,6 +146,8 @@
+```C言語
 void PSET(int x,int y,int col)
@@ -819,3 +821,5 @@
 }
+```

ソースを公開します。

2019/11/01 12:16

投稿

littlestream

スコア37

test CHANGED Viewed

File without changes

test CHANGED Viewed

@@ -139,3 +139,683 @@
 5.うまく行かないので、困っています。どうすれば求まるのか
 誰か教えてください。
+一部のソースコードを公開します。
+void PSET(int x,int y,int col)
+{
+    //点を打つ
+    if(x<0 || x>WIDTH) return;
+    if(y<0 || y>HEIGHT) return;
+    lpPixel[x+(y*WIDTH)]=col;
+}
+inline void ScanEdge(int x1,int y1,int x2,int y2,int c1,int c2)
+{
+	int b1=(c1>>16) & 0xFF;
+	int g1=(c1>>8) & 0xFF;
+	int r1=c1 & 0xFF;
+	int b2=(c2>>16) & 0xFF;
+	int g2=(c2>>8) & 0xFF;
+	int r2=c2 & 0xFF;
+    for(double i=0.0f;i<=1.00f;i+=0.0001f)//0.0～1.0を、0～1024に直す
+    {
+        int ax=x1*(1.0-i)+x2*i;//線形補間の固定小数点処理化(X成分
+        int ay=y1*(1.0-i)+y2*i;//(固定小数点(Y成分
+        double k=(ay-y1)/(double)(y2-y1);//0.0～1.0を0～256に直す
+        if(ay<0 || ay>=HEIGHT) continue;
+		if(minX2[ay]>ax)//minなX(左側の線分)
+		{
+			minX2[ay]=ax;
+			int b=(1.0-k)*b1+k*b2;
+			int g=(1.0-k)*g1+k*g2;
+			int r=(1.0-k)*r1+k*r2;
+			if(r<0) r=0; if(r>255) r=255;
+			if(g<0) g=0; if(g>255) g=255;
+			if(b<0) b=0; if(b>255) b=255;
+			cmin[ay]=RGB(b,g,r);
+			PSET(ax,ay,cmin[ay]);
+		}
+		if(maxX2[ay]<ax)//maxなX（右側の線分)
+		{
+			maxX2[ay]=ax;
+			int b=(1.0-k)*b1+k*b2;
+			int g=(1.0-k)*g1+k*g2;
+			int r=(1.0-k)*r1+k*r2;
+			if(r<0) r=0; if(r>255) r=255;
+			if(g<0) g=0; if(g>255) g=255;
+			if(b<0) b=0; if(b>255) b=255;
+			cmax[ay]=RGB(b,g,r);
+			PSET(ax,ay,cmax[ay]);
+		}
+    }
+}
+inline void TRIANGLE2(int x1,int y1,int x2,int y2,int x3,int y3,int c1,int c2,int c3)
+{
+	for(int i=0;i<HEIGHT;i++)
+	{
+		minX2[i]=65536;
+		maxX2[i]=-65536;
+		cmin[i]=RGB(0,0,0);
+		cmax[i]=RGB(0,0,0);
+	}
+    ScanEdge(x1,y1,x2,y2,c1,c2);
+    ScanEdge(x2,y2,x3,y3,c2,c3);
+    ScanEdge(x3,y3,x1,y1,c3,c1);
+    for(int y=0;y<HEIGHT;y++)
+    {
+        if(minX2[y]==65536 || maxX2[y]==-65536) continue;
+        for(int x=minX2[y];x<=maxX2[y];x++)
+	{
+	int r1=cmin[y] & 0xFF;
+	int g1=cmin[y]>>8 & 0xFF;
+	int b1=cmin[y]>>16 & 0xFF;
+	int r2=cmax[y] & 0xFF;
+	int g2=cmax[y]>>8 & 0xFF;
+	int b2=cmax[y]>>16 & 0xFF;
+	int j=((x-minX2[y])<<8)/(maxX2[y]-minX2[y]+1);
+	int r=r1*((1<<8)-j)+(r2)*j>>8;
+	int g=g1*((1<<8)-j)+(g2)*j>>8;
+	int b=b1*((1<<8)-j)+(b2)*j>>8;
+	int NewR=r;
+	int NewG=g;
+	int NewB=b;
+	if(NewR<0) NewR=0; if(NewR>255) NewR=255;
+	if(NewG<0) NewG=0; if(NewG>255) NewG=255;
+	if(NewB<0) NewB=0; if(NewB>255) NewB=255;
+            PSET(x,y,RGB(NewB,NewG,NewR));
+	}
+    }
+}
+XYZ CalcNormal2(XYZ a,XYZ b,XYZ c)
+{
+	XYZ d={0.0f};
+	XYZ a2,b2;
+	a2.x=a.x-b.x;
+	a2.y=a.y-b.y;
+	a2.z=a.z-b.z;
+	b2.x=a.x-c.x;
+	b2.y=a.y-c.y;
+	b2.z=a.z-c.z;
+	double len=sqrt(a2.x*a2.x+a2.y*a2.y+a2.z*a2.z);
+	a2.x/=len;
+	a2.y/=len;
+	a2.z/=len;
+	double len2=sqrt(b2.x*b2.x+b2.y*b2.y+b2.z*b2.z);
+	b2.x/=len2;
+	b2.y/=len2;
+	b2.z/=len2;
+        //23-32 31-13 12-21
+        //aY*bZ-aZ*bY aZ*bX-aX*bZ aX*bY-aY*bX
+	d.x=a2.y*b2.z-a2.z*b2.y;
+	d.y=a2.z*b2.x-a2.x*b2.z;
+	d.z=a2.x*b2.y-a2.y*b2.x;
+	return d;
+}
+LRESULT CALLBACK WndProc(HWND hwnd, UINT iMsg, WPARAM wParam, LPARAM lParam) {
+	HDC hdc;
+	PAINTSTRUCT ps;
+	switch (iMsg) {
+    case WM_TIMER:
+        {
+            if(GetAsyncKeyState(VK_ESCAPE)<0)
+                SendMessage(hwnd,WM_DESTROY,NULL,NULL);
+            double X=0.0f;
+            double Y=0.0f;
+            double Z=0.0f;
+            //Z+=0.1f;
+            //Y+=0.1f;
+            if(GetAsyncKeyState(VK_RIGHT)<0) Y=-2.5f;
+            else if(GetAsyncKeyState(VK_LEFT)<0) Y=2.5f;
+            else Y=0.0f;
+            if(GetAsyncKeyState(VK_DOWN)<0) X=-2.5f;
+            else if(GetAsyncKeyState(VK_UP)<0) X=2.5f;
+            else X=0.0f;
+			Rotate(X,Y,0);
+            static int X1=0,X2=0;
+			int n=8;
+			InvalidateRect(hwnd,NULL,NULL);
+        }
+        return 0;
+	case WM_CREATE:
+		MyOutputDebugString("起動しました。");
+        Init();//3次元座標を代入
+        SetTimer(hwnd,100,5,NULL);
+		/* BITMAPINFOをゼロクリア */
+		ZeroMemory(&biInfo, sizeof(BITMAPINFO));
+		/* BITMAPINFO設定 */
+		biInfo.bmiHeader.biSize = sizeof(BITMAPINFOHEADER);
+		biInfo.bmiHeader.biWidth = WIDTH;
+		biInfo.bmiHeader.biHeight = -HEIGHT;
+		biInfo.bmiHeader.biPlanes = 1;
+		biInfo.bmiHeader.biBitCount = 32;
+		biInfo.bmiHeader.biCompression = BI_RGB;
+		/* ウインドウのDCを取得 */
+		hdc = GetDC(hwnd);
+		/* biInfoの形式でDIBSectionを作成 */
+		hBMP =CreateDIBSection(hdc, &biInfo, DIB_RGB_COLORS, (LPVOID*)(&lpPixel), NULL, 0);
+		/* DIBSection用のメモリDCを作成 */
+		hdcBMP = CreateCompatibleDC(hdc);
+		/* メモリDCにDIBSectionを選択 */
+		hBMPOLD = (HBITMAP)SelectObject(hdcBMP, hBMP);
+		/* 不要になったウインドウのDCを解放 */
+		ReleaseDC(hwnd, hdc);
+		return 0;
+	case WM_PAINT:
+		{
+		hdc = BeginPaint(hwnd, &ps);
+		ClearScreen();//画面を消す
+		Projection2DAxisFrom3DAxis(hdc);//透視投影
+		XYZ vvv[12];//全ての頂点の法線をクリアする
+		for(int i=0;i<12;i++)
+		{
+			vvv[i].x=0.0f;//Clear
+			vvv[i].y=0.0f;//Clear
+			vvv[i].z=0.0f;//Clear
+			//これが無いとどんどん加算されてしまう
+		}
+		for(int i=0;i<12;i++)
+		{
+			int i1=Poly[i][0];
+			int i2=Poly[i][1];
+			int i3=Poly[i][2];
+				XYZ n=CalcNormal2(Pnt[i1],Pnt[i2],Pnt[i3]);
+				vvv[i1].x+=n.x;
+				vvv[i1].y+=n.y;
+				vvv[i1].z+=n.z;
+				vvv[i2].x+=n.x;
+				vvv[i2].y+=n.y;
+				vvv[i2].z+=n.z;
+				vvv[i3].x+=n.x;
+				vvv[i3].y+=n.y;
+				vvv[i3].z+=n.z;
+				vvv[i1].count++;
+				vvv[i2].count++;
+				vvv[i3].count++;
+					//全ての頂点毎に、全ての面の頂点座標と比較して、
+		}
+		for(int i=0;i<8;i++)
+		{
+			if(vvv[i].count>0)
+			{
+				vvv[i].x/=vvv[i].count ;
+				vvv[i].y/=vvv[i].count ;
+				vvv[i].z/=vvv[i].count ;
+			}
+		}
+		//頂点ベクトルは求まった(はず...)から、頂点カラーを面にある全ての頂点で求める
+		double ViewCos[12];
+		int vpx=200,vpy=-120,vpz=100;
+		for(int i=0;i<12;i++)
+		{
+		double a=sqrt(vvv[i].x*vvv[i].x+vvv[i].y*vvv[i].y+vvv[i].z*vvv[i].z)+0.01f;
+			vvv[i].x/=a;
+			vvv[i].y/=a;
+			vvv[i].z/=a;
+            		double vv=sqrt((vpx*vpx+vpy*vpy+vpz*vpz));//視線ベクトルの長さを求める
+			double v1=vpx/(vv+0.01f);//視線ベクトルの正規化
+			double v2=vpy/(vv+0.01f);//
+			double v3=vpz/(vv+0.01f);//
+			ViewCos[i]=(v1*vvv[i].x+v2*vvv[i].y+v3*vvv[i].z);
+			//視線ベクトルと面ごとの法線ベクトルの内積
+			if(ViewCos[i]<0)
+				TRIANGLE2(Pnt2[Poly[i][0]].x,Pnt2[Poly[i][0]].y,
+					Pnt2[Poly[i][1]].x,Pnt2[Poly[i][1]].y,
+					Pnt2[Poly[i][2]].x,Pnt2[Poly[i][2]].y,
+					RGB(255,255,255),
+					RGB(255,255,255),
+					RGB(255,255,255)
+				);
+			//とりあえず真っ白にしたが、頂点カラーを求めた後にどうすれば
+			//３点の色を指定する処理になるか分からないです。
+		}
+		CalcFPS(hdc);
+		/* DIBSectionをDIBとして描画 */
+		StretchDIBits(hdc, 0, 0, WIDTH, HEIGHT,
+		     0, 0, WIDTH, HEIGHT, lpPixel, &biInfo,
+		     DIB_RGB_COLORS,SRCCOPY);
+		EndPaint(hwnd, &ps);
+		DeleteObject(SelectObject(hdc , GetStockObject(WHITE_BRUSH)));
+		}
+		return 0;
+	case WM_DESTROY: /* ウインドウ破棄時 */
+		SelectObject(hdcBMP, hBMPOLD);
+		DeleteObject(hdcBMP);
+		DeleteObject(hBMP);
+		PostQuitMessage(0);
+		return 0;
+	}
+	return DefWindowProc(hwnd, iMsg, wParam, lParam);
+}