Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  782 1485 (N36E04 L=175)   775    145        3 ( 720/1441- 859/1517)  -5.7 E09-W01 N33-N37
      2 1706 1390 (N07W54 L=233)   465     71        2 (1675/1367-1735/1409) -21.8 W51-W56 N06-N08
      3 1066  860 (S07E01 L=178)  3546    798        7 ( 877/ 609-1194/1049)  72.3 E12-W05 S21-N04
      4  148  910 (N14E66 L=114)   899    172        3 ( 113/ 867- 184/ 951)  70.0 E70-E61 N10-N16
      5  205  887 (N12E59 L=120)   598     73        2 ( 190/ 840- 221/ 932)  80.5 E62-E57 N08-N14
      6  866  708 (S11E16 L=163)   232     80        0 ( 837/ 671- 892/ 749)  38.7 E18-E13 S13-S09
      7 1255  634 (S25W05 L=184)   545    144        0 (1184/ 590-1333/ 676)   5.7 W00-W10 S26-S23
      8  618  347 (S28E46 L=133)  1011    234        1 ( 577/ 272- 652/ 449)  61.7 E50-E42 S34-S21
  

    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