Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1635 1189 (N01W41 L=241)   898    212        1 (1580/1096-1711/1230)   8.1 W34-W48 S02-N02
      2  413 1182 (N29E38 L=161)   515    148        0 ( 338/1153- 473/1207)   0.0 E44-E32 N27-N30
      3 1874 1010 (S17W60 L=259)  4502    278        8 (1826/ 872-1930/1107)  47.5 W52-W69 S23-S10
      4  160 1045 (N24E67 L=132)  1218     83        3 ( 106/1018- 198/1073)   3.4 E77-E60 N22-N25
      5 1313  757 (S15W11 L=210)   310     95        2 (1276/ 710-1341/ 804)  51.3 W07-W12 S17-S12
      6  804  674 (S08E21 L=178)   248     24        1 ( 780/ 649- 822/ 694) -71.6 E21-E19 S09-S06
      7  683  481 (S18E35 L=164)   827    239        1 ( 585/ 445- 782/ 514)  -8.1 E42-E28 S21-S14
      8  583  325 (S27E52 L=147)  3767    277        7 ( 527/ 212- 646/ 424)  48.4 E66-E46 S36-S18
      9  868  351 (S30E29 L=170)   458    192        0 ( 798/ 293- 934/ 409)  20.6 E37-E21 S33-S27
     10  498  337 (S24E60 L=140)  1090    123        5 ( 471/ 296- 539/ 393)  54.5 E64-E54 S27-S20
  

    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