Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1340 1378 (N21W28 L=195)   629    194        1 (1244/1353-1427/1401) -24.4 W21-W33 N18-N24
      2 1573 1072 (S01W35 L=202)  1330    440        3 (1526/ 865-1640/1234)  82.2 W32-W39 S15-N07
      3 1693 1023 (S07W44 L=210)   339     73        0 (1679/ 988-1707/1062)  76.0 W42-W45 S09-S05
      4 1812  931 (S16W54 L=221)  1302    211        5 (1764/ 841-1840/1021)  74.7 W48-W57 S21-S10
      5 1457  776 (S17W22 L=189)   483    213        1 (1414/ 692-1480/ 871)  52.1 W18-W25 S21-S11
      6 1054  815 (S06E02 L=165)   216    117        2 (1008/ 772-1104/ 844)  21.8 E05-W01 S07-S04
      7 1247  698 (S17W07 L=174)   522    200        3 (1160/ 656-1305/ 752)  11.3 W01-W11 S18-S14
      8  989  704 (S12E08 L=159)   201    135        0 ( 934/ 683-1038/ 731)   0.0 E11-E04 S13-S09
      9 1677  601 (S35W41 L=207)  2511    287        5 (1607/ 457-1725/ 711)  66.4 W34-W45 S45-S26
     10  353  520 (S15E61 L=106)  1646    236        4 ( 326/ 422- 395/ 644)  47.0 E70-E55 S21-S06
     11 1148  404 (S35E05 L=162)  3377    602       18 ( 938/ 236-1325/ 593)  24.0 E25-W11 S46-S24
     12  563  545 (S15E40 L=127)   582    160        3 ( 495/ 515- 647/ 562) -10.5 E44-E33 S16-S13
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y