Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  622 1372 (N38E17 L=159)   891    196        3 ( 552/1319- 702/1406)   4.4 E23-E10 N35-N39
      2  160 1208 (N37E61 L=115)   788    258        1 ( 123/1174- 198/1249) -32.0 E65-E56 N34-N39
      3 1809  981 (S19W50 L=226)  2295    431        8 (1744/ 874-1896/1077)  26.6 W42-W61 S23-S15
      4  351 1002 (N21E44 L=132)   982    197        2 ( 294/ 971- 425/1057) -35.0 E49-E38 N18-N25
      5 1428  781 (S18W17 L=193)   591    197        1 (1359/ 712-1504/ 843)   9.5 W11-W23 S20-S15
      6  189  779 (N10E67 L=108)  2527    340        4 ( 141/ 714- 227/ 838) -82.9 E75-E62 N06-N14
      7 1035  765 (S08E06 L=170)   390     96        2 (1019/ 710-1059/ 810)  51.3 E07-E03 S11-S05
      8  632  551 (S11E37 L=139)  2086    389        1 ( 576/ 382- 723/ 722)  73.1 E39-E31 S23-N00
  

    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