Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  659 1731 (N41E45 L= 88)   480     71        1 ( 622/1710- 697/1746) -11.3 E50-E40 N39-N42
      2  831 1663 (N39E26 L=107)   401    106        1 ( 785/1641- 884/1682)  14.0 E29-E21 N36-N40
      3 1274 1436 (N29W11 L=144)  2022    465        3 (1070/1340-1463/1518)  -8.1 E03-W25 N24-N32
      4  455 1130 (S01E38 L= 96)   797    188        3 ( 403/1080- 494/1221)  66.0 E41-E35 S04-N05
      5 1463 1171 (N15W26 L=159)   254    101        0 (1450/1110-1488/1205)  90.0 W24-W27 N11-N17
      6 1782  940 (N05W54 L=187)   323     67        0 (1768/ 893-1794/ 974) -68.2 W52-W55 N02-N07
      7 1572  747 (S09W40 L=173)  1154    413        1 (1439/ 663-1704/ 869)  41.2 W31-W47 S15-N00
      8  737  787 (S18E15 L=119)   257    136        0 ( 680/ 763- 770/ 839)  45.0 E18-E12 S19-S14
      9 1073  760 (S15W07 L=140)   195     78        0 (1029/ 752-1110/ 770)  21.8 W03-W08 S15-S13
     10 1677  601 (S17W54 L=187)   563    112        0 (1652/ 559-1699/ 661) -81.9 W51-W55 S19-S12
     11 1513  517 (S24W43 L=176)  1391    323        2 (1365/ 487-1652/ 551)  16.9 W31-W54 S27-S20
  

    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