Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  637 1801 (N59E31 L=191)   643     82        0 ( 602/1765- 675/1820)  38.7 E36-E25 N55-N61
      2 1637 1615 (N30W57 L=280)  1453    232        0 (1550/1528-1696/1672) -55.0 W51-W61 N23-N34
      3  949 1614 (N37W01 L=224)   539    110        0 ( 915/1580- 999/1639)   9.5 E01-W05 N35-N39
      4 1244 1600 (N33W23 L=245)   229     50        1 (1216/1582-1269/1618) -56.3 W20-W24 N31-N34
      5  346  886 (S02E46 L=176)   886    254        0 ( 257/ 785- 429/ 948) -36.9 E53-E40 S09-N02
      6  259  636 (S16E63 L=160)   842     93        3 ( 222/ 601- 292/ 677)  15.3 E69-E58 S18-S13
      7 1715  654 (S30W49 L=271)   308     38        2 (1695/ 632-1733/ 677) -31.0 W46-W51 S31-S28
      8  331  581 (S20E57 L=165)   362     75        2 ( 313/ 548- 354/ 609)  76.0 E59-E53 S22-S18
      9  492  534 (S25E44 L=178)   810    270        0 ( 426/ 443- 597/ 588)  20.9 E54-E33 S30-S22
     10 1208  302 (S52W06 L=229)   474     85        1 (1166/ 284-1252/ 318)   7.1 W02-W10 S52-S49
     11  955  187 (S59E24 L=198)  1374    167        2 ( 893/ 157-1014/ 229)   3.4 E33-E16 S62-S55
  

    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