Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1302 1677 (N35W38 L=133)  1081    174        1 (1228/1650-1392/1704) -17.1 W31-W44 N31-N37
      2 1563 1599 (N23W53 L=149)   482     69        0 (1526/1568-1599/1622)   6.3 W48-W57 N22-N24
      3 1110 1339 (N19W12 L=107)  1404    422        2 ( 923/1277-1276/1448)   0.0 W00-W25 N16-N21
      4  579 1386 (N34E22 L= 73)  1771    317        3 ( 451/1301- 702/1445) -32.3 E32-E13 N25-N38
      5  149 1269 (N35E70 L= 25)   783     75        1 ( 131/1236- 178/1304) -24.4 E77-E63 N32-N37
      6  345 1210 (N27E43 L= 53)   356     86        1 ( 331/1163- 362/1247) -90.0 E43-E40 N24-N29
      7 1849 1218 (S05W62 L=157)   406     45        0 (1834/1190-1869/1240)  31.0 W59-W65 S06-S03
      8 1736  882 (S21W44 L=139)  1749    363        1 (1672/ 753-1812/1047)  45.0 W37-W53 S28-S13
      9 1884  833 (S29W62 L=158)   522     59        2 (1863/ 810-1904/ 862)  21.8 W59-W65 S29-S26
     10  869  727 (S11E15 L= 80)   634    172        1 ( 849/ 647- 898/ 826)  55.0 E18-E11 S15-S05
     11 1393  657 (S27W15 L=110)   477     90        0 (1359/ 607-1421/ 706)  51.3 W12-W17 S30-S24
     12 1330  523 (S34W08 L=104)   228     61        1 (1300/ 506-1359/ 539) -14.0 W06-W10 S35-S32
  

    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