Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1061 1734 (N38E22 L=133)  1209    250        2 ( 961/1667-1176/1798)   0.0 E32-E12 N35-N40
      2  691 1726 (N26E46 L=109)   706    194        1 ( 664/1668- 741/1792) -45.0 E51-E39 N22-N31
      3  555 1678 (N20E53 L=102)   717    104        1 ( 518/1650- 609/1711)  -4.8 E58-E46 N18-N21
      4 1125 1607 (N31E12 L=143)   308     16        0 (1108/1579-1141/1634) -63.4 E13-E10 N28-N32
      5  485 1408 (N02E43 L=112)  2257    376        7 ( 414/1268- 605/1535)  86.4 E46-E38 S05-N11
      6 1285 1456 (N26W03 L=158)   389    105        1 (1258/1407-1315/1501) -45.0 W00-W05 N23-N28
      7  817 1282 (N03E18 L=136)   714    270        3 ( 774/1178- 860/1383) -48.4 E22-E13 S01-N08
      8 1703  995 (N13W42 L=197)   364    105        0 (1693/ 942-1713/1053) -71.6 W39-W43 N09-N15
      9 1914  762 (N08W76 L=231)  1137     75        0 (1895/ 729-1938/ 795)  45.0 W72-W80 N06-N11
     10 1524  375 (S25W58 L=213)  1124    141        1 (1461/ 338-1593/ 411)  53.1 W54-W62 S29-S21
     11 1234  152 (S49W67 L=221)  1276    132        0 (1180/ 115-1290/ 188)   4.2 W55-W85 S50-S46
  

    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