Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  807 1809 (N37E40 L=272)   669     67        1 ( 772/1786- 848/1835)   0.0 E43-E35 N35-N39
      2  672 1437 (N09E31 L=281)  2069    386        5 ( 626/1220- 728/1654) -61.7 E36-E22 S03-N23
      3  300 1493 (N06E63 L=248)  1133    318        0 ( 253/1456- 352/1541) -16.4 E73-E56 N04-N09
      4 1492 1211 (N15W23 L=334)   438    114        0 (1461/1152-1517/1278) -54.5 W19-W25 N11-N18
      5  359 1145 (S14E45 L=266)  1022    331        2 ( 244/1119- 477/1165)  11.9 E55-E36 S15-S11
      6  113 1018 (S24E68 L=244)  1345    210        1 (  81/ 972- 133/1057)  45.0 E70-E65 S26-S21
      7 1833  925 (N09W56 L=  8)   729     91        0 (1810/ 877-1853/ 972) -45.0 W52-W59 N06-N11
      8 1589  641 (S13W46 L=357)   798    177        1 (1513/ 589-1659/ 722)  67.4 W42-W48 S18-S06
      9 1641  374 (S25W66 L= 17)  1276    195        2 (1589/ 348-1679/ 414)  35.5 W61-W70 S27-S22
  

    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