基因資料處理121之SSW的score matrix調整，使得與SparkSW評分一致

阿新 • • 發佈：2018-12-25

基因資料處理系列

1.解釋

SSW的評分矩陣是128*128的，是按char的int值來進行計算的。而blosum50是蛋白質的，而且不是按ABC順序來的，所以需要轉換。

java中無法執行128*128的String矩陣，所以需要提供blosum的轉換

2.程式碼：

    public static int[][] blosum50() {
        int[][] a = {

        };
        return a;
    }

    public static int[][] blosum50(String[][] strArr) {
        int[][] result = new int[128][128];

        int m = strArr.length;
        int n = strArr[0].length;

        for (int i = 0;i< result.length;i++) {
            for (int j = 0;j<  result[i].length;j++) {
                result[i][j] = Integer.valueOf(strArr[m - 1][n - 1]);
            }
        }
        for (int i = 0 ;i< m;i++) {
             char c = strArr[i][0].charAt(0);
            for (int j = 0 ;j<  result.length;j++) {
                result[Integer.valueOf(c)][j] = Integer.valueOf(strArr[i][strArr[i].length - 1]);
                result[j][Integer.valueOf(c)] = Integer.valueOf(strArr[i][strArr[i].length - 1]);
            }
        }

        //    println(result(0)(0))
        for (int i = 0 ;i< m;i++) {
            char strI = strArr[i][0].charAt(0);
            //       strArr(i)(0).toCharArray

            int valueI = Integer.valueOf(strArr[i][0].charAt(0));
            //      println(strI + "=" + valueI + "=" + strArr(i)(0).charAt(0).toInt)

            for (int j = 1 ;j< n;j++) {

                int valueJ =Integer.valueOf(strArr[j - 1][0].charAt(0));
                int tmp = Integer.valueOf(strArr[i][j]);
                result[valueI][valueJ] = tmp;
                //        arr(value)(j)

            }
        }

        //println(result(0)(0))
        return result;
    }
    public static String[][] blosum50Base() {
        String[][] a = {
                {"A","5","-2","-1","-2","-1","-1","-1","0","-2","-1","-2","-1","-1","-3","-1","1","0","-3","-2","0","-2","-2","-1","-1","-5"},
                {"R","-2","7","-1","-2","-4","1","0","-3","0","-4","-3","3","-2","-3","-3","-1","-1","-3","-1","-3","-1","-3","0","-1","-5"},
                {"N","-1","-1","7","2","-2","0","0","0","1","-3","-4","0","-2","-4","-2","1","0","-4","-2","-3","5","-4","0","-1","-5"},
                {"D","-2","-2","2","8","-4","0","2","-1","-1","-4","-4","-1","-4","-5","-1","0","-1","-5","-3","-4","6","-4","1","-1","-5"},
                {"C","-1","-4","-2","-4","13","-3","-3","-3","-3","-2","-2","-3","-2","-2","-4","-1","-1","-5","-3","-1","-3","-2","-3","-1","-5"},
                {"Q","-1","1","0","0","-3","7","2","-2","1","-3","-2","2","0","-4","-1","0","-1","-1","-1","-3","0","-3","4","-1","-5"},
                {"E","-1","0","0","2","-3","2","6","-3","0","-4","-3","1","-2","-3","-1","-1","-1","-3","-2","-3","1","-3","5","-1","-5"},
                {"G","0","-3","0","-1","-3","-2","-3","8","-2","-4","-4","-2","-3","-4","-2","0","-2","-3","-3","-4","-1","-4","-2","-1","-5"},
                {"H","-2","0","1","-1","-3","1","0","-2","10","-4","-3","0","-1","-1","-2","-1","-2","-3","2","-4","0","-3","0","-1","-5"},
                {"I","-1","-4","-3","-4","-2","-3","-4","-4","-4","5","2","-3","2","0","-3","-3","-1","-3","-1","4","-4","4","-3","-1","-5"},
                {"L","-2","-3","-4","-4","-2","-2","-3","-4","-3","2","5","-3","3","1","-4","-3","-1","-2","-1","1","-4","4","-3","-1","-5"},
                {"K","-1","3","0","-1","-3","2","1","-2","0","-3","-3","6","-2","-4","-1","0","-1","-3","-2","-3","0","-3","1","-1","-5"},
                {"M","-1","-2","-2","-4","-2","0","-2","-3","-1","2","3","-2","7","0","-3","-2","-1","-1","0","1","-3","2","-1","-1","-5"},
                {"F","-3","-3","-4","-5","-2","-4","-3","-4","-1","0","1","-4","0","8","-4","-3","-2","1","4","-1","-4","1","-4","-1","-5"},
                {"P","-1","-3","-2","-1","-4","-1","-1","-2","-2","-3","-4","-1","-3","-4","10","-1","-1","-4","-3","-3","-2","-3","-1","-1","-5"},
                {"S","1","-1","1","0","-1","0","-1","0","-1","-3","-3","0","-2","-3","-1","5","2","-4","-2","-2","0","-3","0","-1","-5"},
                {"T","0","-1","0","-1","-1","-1","-1","-2","-2","-1","-1","-1","-1","-2","-1","2","5","-3","-2","0","0","-1","-1","-1","-5"},
                {"W","-3","-3","-4","-5","-5","-1","-3","-3","-3","-3","-2","-3","-1","1","-4","-4","-3","15","2","-3","-5","-2","-2","-1","-5"},
                {"Y","-2","-1","-2","-3","-3","-1","-2","-3","2","-1","-1","-2","0","4","-3","-2","-2","2","8","-1","-3","-1","-2","-1","-5"},
                {"V","0","-3","-3","-4","-1","-3","-3","-4","-4","4","1","-3","1","-1","-3","-2","0","-3","-1","5","-3","2","-3","-1","-5"},
                {"B","-2","-1","5","6","-3","0","1","-1","0","-4","-4","0","-3","-4","-2","0","0","-5","-3","-3","6","-4","1","-1","-5"},
                {"J","-2","-3","-4","-4","-2","-3","-3","-4","-3","4","4","-3","2","1","-3","-3","-1","-2","-1","2","-4","4","-3","-1","-5"},
                {"Z","-1","0","0","1","-3","4","5","-2","0","-3","-3","1","-1","-4","-1","0","-1","-2","-2","-3","1","-3","5","-1","-5"},
                {"X","-1","-1","-1","-1","-1","-1","-1","-1","-1","-1","-1","-1","-1","-1","-1","-1","-1","-1","-1","-1","-1","-1","-1","-1","-5"},
                {"*","-5","-5","-5","-5","-5","-5","-5","-5","-5","-5","-5","-5","-5","-5","-5","-5","-5","-5","-5","-5","-5","-5","-5","-5","1"}
        };
        return a;
    }

3.結果：

調整後與SparkSW執行基本一致，但是ref Data比較大時會不同，比如：待分析原因解決

SparkSSW
time:10
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
SSW
time:1
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)

記錄

[email protected]:~/disk2/xubo/project/alignment/SparkSW/SparkSW20161114/alluxio-1.3.0$ tail -f sparkSWSSWtime201611270135.txt 
SparkSW:
time:1
(1734,UniRef100_Q6GZX4)                                                         
time:2
(1734,UniRef100_Q6GZX4)                                                         
time:3
(1734,UniRef100_Q6GZX4)                                                         
time:4
(1734,UniRef100_Q6GZX4)                                                         
time:5
(1734,UniRef100_Q6GZX4)                                                         
time:6
(1734,UniRef100_Q6GZX4)                                                         
time:7
(1734,UniRef100_Q6GZX4)                                                         
time:8
(1734,UniRef100_Q6GZX4)                                                         
time:9
(1734,UniRef100_Q6GZX4)                                                         
time:10
(1734,UniRef100_Q6GZX4)                                                         
SSW
time:1
(1734,UniRef100_Q6GZX4)                                                         
16/11/27 01:49:22 WARN QueuedThreadPool: 1 threads could not be stopped
time:2
(1734,UniRef100_Q6GZX4)                                                         
time:3
(1734,UniRef100_Q6GZX4)                                                         
time:4
(1734,UniRef100_Q6GZX4)                                                         
time:5
(1734,UniRef100_Q6GZX4)                                                         
16/11/27 01:51:26 WARN QueuedThreadPool: 4 threads could not be stopped
time:6
(1734,UniRef100_Q6GZX4)                                                         
time:7
(1734,UniRef100_Q6GZX4)                                                         
time:8
(1734,UniRef100_Q6GZX4)                                                         
time:9
(1734,UniRef100_Q6GZX4)                                                         
time:10
(1734,UniRef100_Q6GZX4)                                                         
SparkSW:
time:1
(1734,UniRef100_Q6GZX4)                                                         
time:2
(1734,UniRef100_Q6GZX4)                                                         
time:3
(1734,UniRef100_Q6GZX4)                                                         
time:4
(1734,UniRef100_Q6GZX4)                                                         
time:5
(1734,UniRef100_Q6GZX4)                                                         
time:6
(1734,UniRef100_Q6GZX4)                                                         
time:7
(1734,UniRef100_Q6GZX4)                                                         
time:8
(1734,UniRef100_Q6GZX4)                                                         
time:9
(1734,UniRef100_Q6GZX4)                                                         
time:10
(1734,UniRef100_Q6GZX4)                                                         
SSW
time:1
(1734,UniRef100_Q6GZX4)                                                         
16/11/27 02:19:43 WARN QueuedThreadPool: 2 threads could not be stopped
time:2
(1734,UniRef100_Q6GZX4)                                                         
time:3
(1734,UniRef100_Q6GZX4)                                                         
time:4
(1734,UniRef100_Q6GZX4)                                                         
time:5
(1734,UniRef100_Q6GZX4)                                                         
time:6
(1734,UniRef100_Q6GZX4)                                                         
time:7
(1734,UniRef100_Q6GZX4)                                                         
time:8
(1734,UniRef100_Q6GZX4)                                                         
16/11/27 02:25:02 WARN QueuedThreadPool: 2 threads could not be stopped
time:9
(1734,UniRef100_Q6GZX4)                                                         
16/11/27 02:25:48 WARN QueuedThreadPool: 5 threads could not be stopped
time:10
(1734,UniRef100_Q6GZX4)                                                         
SparkSW:
time:1
(1734,UniRef100_Q6GZX4)                                                         
time:2
(1734,UniRef100_Q6GZX4)                                                         
time:3
(1734,UniRef100_Q6GZX4)                                                         
time:4
(1734,UniRef100_Q6GZX4)                                                         
time:5
(1734,UniRef100_Q6GZX4)                                                         
time:6
(1734,UniRef100_Q6GZX4)                                                         
time:7
(1734,UniRef100_Q6GZX4)                                                         
time:8
(1734,UniRef100_Q6GZX4)                                                         
time:9
(1734,UniRef100_Q6GZX4)                                                         
time:10
(1734,UniRef100_Q6GZX4)                                                         
SSW
time:1
(1734,UniRef100_Q6GZX4)                                                         
16/11/27 03:12:40 WARN QueuedThreadPool: 2 threads could not be stopped
time:2
(1734,UniRef100_Q6GZX4)                                                         
time:3
(1734,UniRef100_Q6GZX4)                                                         
16/11/27 03:15:08 WARN QueuedThreadPool: 4 threads could not be stopped
time:4
(1734,UniRef100_Q6GZX4)                                                         
16/11/27 03:16:24 WARN QueuedThreadPool: 1 threads could not be stopped
time:5
(1734,UniRef100_Q6GZX4)                                                         
time:6
(1734,UniRef100_Q6GZX4)                                                         
time:7
(1734,UniRef100_Q6GZX4)                                                         
16/11/27 03:20:13 WARN QueuedThreadPool: 2 threads could not be stopped
time:8
(1734,UniRef100_Q6GZX4)                                                         
time:9
(1734,UniRef100_Q6GZX4)                                                         
16/11/27 03:22:45 WARN QueuedThreadPool: 4 threads could not be stopped
time:10
(1734,UniRef100_Q6GZX4)                                                         
SparkSW:
time:1
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:2
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:3
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:4
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:5
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:6
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:7
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:8
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:9
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:10
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
SSW
time:1
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:2
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:3
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:4
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:5
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:6
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:7
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:8
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:9
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:10
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
SparkSW:
time:1
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:2
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:3
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:4
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:5
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:6
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:7
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:8
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:9
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:10
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
SSW
time:1
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:2
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:3
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:4
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:5
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:6
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:7
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:8
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:9
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:10
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
SparkSW:
time:1
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:2
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:3
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:4
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:5
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:6
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:7
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:8
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:9
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:10
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
SSW
time:1
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:2
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:3
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:4
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:5
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:6
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:7
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:8
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:9
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
16/11/27 12:57:11 WARN QueuedThreadPool: 1 threads could not be stopped
time:10
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
SparkSW:
time:1
[Stage 1:====================================>                  (86 + 16) / 128]Exception in thread "main" org.apache.spark.SparkException: Job cancelled because SparkContext was shut down

參考

    【1】https://github.com/xubo245
    【2】http://blog.csdn.net/xubo245/

基因資料處理121之SSW的score matrix調整，使得與SparkSW評分一致

基因資料處理系列 1.解釋 SSW的評分矩陣是128*128的，是按char的int值來進行計算的。而blosum50是蛋白質的，而且不是按ABC順序來的，所以需要轉換。 java中無法執行128*128的String矩陣，所以需要提供blosum

基因資料處理122之SSW和SparkSW評分不一致，query為Q9

基因資料處理系列 1.解釋 RT,但是順序一致 2.程式碼： [email protected]:~/disk2/xubo/project/alignment/SparkSW/SparkSW20161114/alluxio-1.3.

基因資料處理123之SSW程式碼不正確，到時比SparkSW時間長

基因資料處理系列 1.解釋由於要生成新的score matrix：blosum50，第一次使用靜態方法，直接傳給align，到時每次執行都需要進行一次score matrix的計算，而這個是將blosum50的矩陣轉換成128*128的矩陣，當計算

基因資料處理114之BWA建立全基因組索引成功

執行記錄 [email protected]:~/disk2/home/hadoop/xubo/ref/buildIndex$ bwa index GCA_000001405.15_GRCh38_full_analysis_set.fna [bwa

基因資料處理119之java呼叫SSW在linux下執行

基因資料處理系列 1.解釋測試自帶Example： [email protected]:~/xubo/tools/Complete-Striped-Smith-Waterman

基因資料處理120之scala呼叫SSW在linux下執行

基因資料處理系列 1.解釋先有java提供轉換，使用jni呼叫c 然後scala呼叫java 2.程式碼： 2.1 java： package ssw; /** * Created by xubo on 2016/11/25.

基因資料處理117之重新多次執行SparkBWA Yarn叢集

基因資料處理系列 1.解釋重新執行，跟換了檔案地址後來終端，需要執行其他的 2.程式碼： [email protected]:~/disk2/xubo/project/alignment/sparkBWA$ vi sparkB

基因資料處理115之重新執行SparkBWA本地模式

基因資料處理系列 1.記錄 [email protected]:~/disk2/xubo/tools/SparkBWA/build$ ./pairedGRCH38L1Local.sh [Stage 3:>

基因資料處理113之對avocado識別的SparkBWA變異資料進行疾病分析_

（一）問題問題1：avocado資料讀取： avocado存到磁碟是：RDD[Genotype] val processedGenotypes: RDD[Genotype] = postProcessVariants(calledVariants, st

基因資料處理118之SSW執行

基因資料處理系列 1.解釋 SSW是一個更快的SW演算法，並且提供了c語言lib和java的呼叫程式碼： https://github.com/mengyao/Complete-Str

基因資料處理116之重新執行SparkBWA Yarn叢集模式

基因資料處理系列 1.解釋很久沒執行SparkBWA了，系統檔案有點多，重新執行。 2.程式碼： end [email protected]:~/disk2/xubo/

基因資料處理5之GRCH38資料來源和檢視資訊

資料來源： ftp://ftp.ncbi.nlm.nih.gov/genomes/all/GCA_000001405.15_GRCh38/seqs_for_alignment_pipelines.ucsc_ids/GCA_000001405.15_GRCh38_full_

基因資料處理6之BWA_MEM無法分配記憶體

基因資料處理之BWA_MEM無法分配記憶體：建立BWA索引的時候記憶體不足，現在用BWA-MEM又記憶體不足，真耗記憶體 [email protected]:~/cloud/

大資料處理框架之:Storm + Kafka + zookeeper 叢集

Storm kafka zookeeper 叢集我們知道storm的作用主要是進行流式計算，對於源源不斷的均勻資料流流入處理是非常有效的，而現實生活中大部分場景並不是均勻的資料流，而是時而多時而少的資料流入，這種情況下顯然用批量處理是不合適的，如果使用storm做實時計算的話可能因為資

大資料處理基礎之利用hadoop寫的簡單mapreduce案例

案例：需要處理的資料： 13877779999 bj zs 2145 13766668888 sh ls 1028 13766668888 sh ls 9987 13877779999 bj zs 5678 13544445555 sz ww 10577 1387777999

大資料處理過程之核心技術ETL詳解

核心技術架構挑戰： 1、對現有資料庫管理技術的挑戰。 2、經典資料庫技術並沒有考慮資料的多類別（variety）、SQL（結構化資料查詢語言），在設計的一開始是沒有考慮到非結構化資料的儲存問題。 3、實時性技術的挑戰：一般而言，傳統資料倉庫系統，BI應用，對處理時間的要求

數學建模資料處理模型之變數相關性類（灰色相關聯、相關性分析）

相關類灰色關聯 1作用：系統分析主要因素；次要因素，因素對系統發展的影響，以便對各因素強化發展或者抑制發展。 2 灰色關聯分析的基本思想：根據序列曲線的幾何形狀的相似程度判斷其聯絡緊密性 3 具體操作步驟：（1）繪圖：各指標，各系統的發展趨勢（2）確定分析數列：母序列：能反映系統行為特徵的資料序列。（

數字影象處理成長之路4： C語言與離散傅立葉變換(DFT)

這幾天一直學習傅立葉變換，看了很多國內外資料，網上講原理的很多，到了程式實現這塊大多是Matlab，opencv等，這些軟體的api對於我們理解DFT在計算機中的實現並沒有多大幫助。於是想用C/C++實現DFT，經過不斷的閱讀與程式設計實驗，最終程式有了還算滿意

Matlab之山東大學數字影象處理實驗一（影象載入，顯示與合成）

本文屬於原創，轉載請註明出處。實驗1.1：影象載入、顯示利用影象庫的功能，實現從檔案載入影象，並在視窗中進行顯示的功能；利用常見的影象檔案格式（.jpg; .png; .bmp; .gif）進行測試；實驗1.2：影象合成現有一張4通道透明影象a.p

資料結構基礎之圖的鄰接矩陣實現與鄰接表實現

【鄰接矩陣】鄰接矩陣，就是一個反應邊與邊之間聯絡的二維陣列。這個二維陣列我們用matrix[numV][numV]表示，其中numV是頂點數。對於無權圖若頂點Vi和Vj之間有邊，則matrix[Vi][Vj]=1;否則matrix[Vi][Vj]=0。對

基因資料處理121之SSW的score matrix調整，使得與SparkSW評分一致

1.解釋

2.程式碼：

3.結果：

記錄

相關推薦