Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  643 1572 (N34E35 L=178)   967    159        2 ( 578/1528- 710/1607)  -8.7 E42-E28 N31-N35
      2  339 1532 (N28E64 L=149)  1919    158        1 ( 290/1461- 384/1591) -14.0 E82-E54 N23-N30
      3  154  903 (S13E70 L=144)   625     71        3 ( 136/ 869- 174/ 931)  26.6 E74-E65 S15-S11
      4  576  835 (S13E28 L=186)   418    138        0 ( 548/ 785- 607/ 894)  60.3 E29-E25 S16-S09
      5  959  769 (S14E02 L=211)   567    208        0 ( 866/ 745-1054/ 793)   4.8 E08-W04 S15-S13
      6 1546  657 (S17W39 L=252)   465     83        2 (1504/ 635-1601/ 678)  -6.3 W35-W44 S18-S16
      7 1468  595 (S22W35 L=248)   358     95        1 (1448/ 557-1489/ 640) -51.3 W32-W36 S24-S19
      8 1541  436 (S33W48 L=262)   400     49        1 (1515/ 421-1572/ 457)  -8.1 W45-W52 S33-S31
  

    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