Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1591 1455 (N28W47 L=301)  1472    242        3 (1529/1385-1669/1552) -67.4 W41-W52 N22-N34
      2  343 1420 (N32E55 L=200)   717    111        2 ( 295/1392- 386/1451) -15.9 E63-E49 N29-N33
      3  276 1340 (N26E60 L=195)  1509    187        2 ( 232/1310- 347/1407) -15.5 E68-E50 N24-N30
      4  648 1240 (N20E24 L=230)   354    169        1 ( 582/1206- 712/1279)  12.5 E28-E19 N17-N21
      5  989  916 (S02E03 L=252)   198     78        0 ( 962/ 894-1027/ 935)  26.6 E04-W00 S03-N00
      6 1640  873 (S09W41 L=295)   627    128        1 (1616/ 822-1662/ 933)  63.4 W38-W42 S11-S05
      7 1749  858 (S11W51 L=305)   492     32        1 (1725/ 835-1773/ 884)  11.3 W48-W53 S12-S09
      8  563  560 (S23E36 L=219)   677    134        3 ( 533/ 500- 598/ 629)  71.6 E37-E33 S27-S18
      9  377  549 (S24E54 L=201)   604     98        1 ( 343/ 523- 410/ 588)  16.7 E58-E48 S26-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