Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  476 1684 (N40E47 L=144)   735    138        2 ( 432/1658- 512/1707) -12.5 E51-E42 N37-N42
      2  785 1687 (N40E20 L=170)   360     67        1 ( 752/1667- 820/1705) -18.4 E22-E16 N37-N41
      3  218 1212 (N08E58 L=132)   660    129        3 ( 188/1171- 244/1268)  80.5 E60-E54 N05-N11
      4  164 1127 (N03E64 L=126)  2386    374        7 ( 116/1031- 207/1212) -73.3 E70-E58 S02-N08
      5 1850 1073 (N02W59 L=249)   682    101        1 (1822/1031-1879/1121) -38.7 W55-W62 N00-N04
      6  988  858 (S14E02 L=188)   425    140        0 ( 958/ 784-1032/ 927)  71.6 E03-W00 S18-S09
      7  489  558 (S33E39 L=151)   444    119        0 ( 447/ 530- 543/ 581)   6.3 E43-E34 S34-S30
  

    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