Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1593 1521 (N18W52 L= 15)   812    122        2 (1531/1488-1639/1545) -45.0 W48-W55 N14-N20
      2  608 1197 (N27E21 L=301)   763    223        0 ( 531/1114- 686/1261) -56.3 E26-E17 N20-N32
      3 1256 1239 (N13W19 L=341)   202     42        1 (1233/1217-1278/1258) -90.0 W17-W20 N11-N14
      4 1138  876 (S04W03 L=325)   490    182        1 (1085/ 813-1218/ 912)   0.0 E01-W07 S06-S01
      5  413  874 (N12E42 L=280)  1012    339        2 ( 293/ 841- 547/ 901) -19.3 E52-E32 N07-N14
      6 1332  802 (S13W12 L=334)   362    154        1 (1260/ 786-1413/ 819) -29.1 W07-W16 S16-S10
      7 1338  658 (S22W09 L=331)   296     70        1 (1304/ 631-1378/ 693)   9.5 W06-W12 S23-S21
      8  969  596 (S16E14 L=308)   509     99        3 ( 934/ 557-1012/ 646) -81.9 E16-E13 S19-S12
      9  423  502 (S10E55 L=267)  1032    171        5 ( 382/ 444- 461/ 556)  63.4 E58-E50 S13-S05
     10  851  395 (S26E29 L=293)   326     79        1 ( 806/ 384- 891/ 409)  -9.5 E32-E26 S27-S24
     11  521  296 (S26E60 L=263)  1202    237        3 ( 500/ 252- 550/ 336)  35.0 E64-E54 S29-S22
  

    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