Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  844 1435 (N26E16 L= 86)  1775    417        4 ( 684/1301- 960/1552) -26.6 E30-E06 N18-N32
      2  283 1447 (N21E62 L= 40)   715    212        0 ( 240/1425- 324/1469)   0.0 E68-E57 N19-N21
      3  671  885 (S10E21 L= 81)   234     51        2 ( 651/ 858- 689/ 918)  76.0 E22-E19 S11-S07
      4  456  823 (S16E36 L= 66)   311    109        0 ( 431/ 783- 501/ 850)  38.7 E38-E32 S18-S14
      5  303  689 (S27E51 L= 51)   563    194        0 ( 257/ 658- 349/ 739)  28.6 E56-E45 S29-S22
      6 1101  649 (S21W08 L=111)   277     86        0 (1082/ 603-1121/ 687)  80.5 W07-W09 S24-S18
      7 1617  538 (S24W49 L=152)   419    101        2 (1585/ 513-1657/ 568)  38.7 W46-W52 S26-S22
      8 1036  506 (S32W06 L=109)   226     55        0 (1001/ 487-1063/ 521) -11.3 W03-W08 S32-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