Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1275 1843 (N46W55 L=254)   432     54        0 (1256/1817-1297/1860)   6.3 W50-W59 N44-N47
      2 1250 1604 (N31W32 L=230)   732    325        0 (1096/1583-1373/1626) -30.5 W21-W38 N26-N36
      3  886 1548 (N38W05 L=204)   689    204        0 ( 799/1519- 984/1577)  -8.7 E01-W12 N36-N39
      4 1625 1447 (N12W50 L=249)   988    158        2 (1566/1386-1678/1512) -77.5 W47-W52 N07-N16
      5  165 1200 (N32E65 L=134)   584     61        2 ( 130/1176- 204/1222) -11.3 E73-E58 N30-N34
      6  503 1029 (N15E31 L=167)   235     64        1 ( 468/1010- 537/1049) -36.9 E33-E29 N13-N16
      7 1571  835 (S21W29 L=228)   715    182        0 (1530/ 769-1591/ 937)  52.1 W25-W32 S23-S14
      8 1258  789 (S15W08 L=207)   196     81        0 (1218/ 775-1298/ 808) -31.0 W05-W10 S16-S13
      9 1795  616 (S43W53 L=252)  1618    274        1 (1699/ 544-1879/ 721)  17.7 W40-W64 S44-S37
     10  861  522 (S22E23 L=176)   348    139        0 ( 829/ 461- 904/ 572)  36.9 E26-E18 S25-S19
     11  672  447 (S22E38 L=161)   890    220        2 ( 639/ 345- 716/ 540)  58.4 E41-E32 S29-S16
     12 1376  345 (S50W06 L=205)  1516    303        1 (1250/ 292-1515/ 413)  -2.3 E05-W21 S51-S46
  

    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