Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1134 1927 (N60W58 L=265)   671     47        1 (1113/1897-1153/1954)  10.3 W49-W72 N57-N62
      2  905 1834 (N60W19 L=226)  1344    264        2 ( 818/1766-1005/1900)  12.3 W06-W38 N55-N63
      3 1644 1061 (S09W39 L=246)   465    123        1 (1609/1003-1680/1108)  29.7 W35-W42 S11-S06
      4 1780  679 (S36W50 L=258)   864    112        0 (1744/ 652-1831/ 719) -21.0 W45-W58 S37-S32
      5  623  552 (S17E36 L=171)   313    136        1 ( 593/ 502- 642/ 614)  54.5 E39-E33 S20-S13
      6  862  554 (S22E20 L=187)   284     93        0 ( 822/ 526- 909/ 572)   9.5 E23-E17 S22-S20
      7  798  428 (S29E29 L=178)   236     62        1 ( 775/ 404- 819/ 464)  31.0 E31-E26 S30-S26
      8  684  303 (S35E46 L=162)  1047    210        2 ( 594/ 280- 787/ 340) -28.3 E51-E37 S38-S30
      9 1403  280 (S58W12 L=220)   434     67        0 (1363/ 267-1440/ 298)  -5.7 W07-W17 S59-S56
  

    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