Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1646 1657 (N20W66 L= 48)  1453    221        3 (1589/1624-1712/1686) -30.3 W59-W72 N16-N23
      2 1839 1320 (S03W65 L= 48)  1495    153        3 (1812/1263-1869/1387)  71.6 W61-W69 S06-N01
      3  384 1300 (N38E37 L=306)   287     60        0 ( 355/1282- 403/1327)  11.3 E39-E34 N37-N39
      4 1729 1158 (S07W47 L= 30)   391     99        1 (1704/1109-1744/1208)  80.5 W44-W49 S09-S03
      5  704 1088 (N17E17 L=326)   222     88        0 ( 660/1074- 755/1108)  -9.5 E19-E13 N16-N18
      6  333  812 (N10E50 L=293)  2631    706        5 ( 152/ 757- 550/ 860) -14.4 E70-E34 N03-N14
      7  880  679 (S10E17 L=325)   194     65        0 ( 844/ 664- 918/ 701)   0.0 E19-E14 S10-S08
      8 1125  621 (S20E05 L=337)   260    112        2 (1106/ 584-1145/ 655)  36.9 E07-E03 S21-S17
      9  674  367 (S23E43 L=300)  1728    387        5 ( 637/ 245- 712/ 517)  51.7 E50-E35 S32-S13
     10 1021  253 (S42E28 L=314)   338     85        1 (1005/ 224-1040/ 280)  78.7 E30-E25 S44-S39
  

    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