Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1183 1863 (N53W58 L=262)   601    117        0 (1131/1841-1242/1876) -12.1 W49-W65 N50-N55
      2 1181 1676 (N41W34 L=239)   559     94        1 (1136/1660-1230/1692) -15.9 W30-W37 N38-N43
      3 1443 1518 (N23W42 L=247)   478     92        2 (1401/1482-1482/1553) -68.2 W40-W44 N20-N26
      4 1309 1400 (N21W28 L=232)   939    333        1 (1168/1323-1441/1458) -47.1 W20-W33 N12-N27
      5 1006 1236 (N20W04 L=209)  2609    611        4 ( 792/1094-1218/1278) -42.0 E09-W14 N06-N25
      6  407 1183 (N30E38 L=166)   571    157        2 ( 371/1119- 450/1248) -66.0 E41-E35 N25-N34
      7 1761 1010 (S14W47 L=252)   349     71        0 (1743/ 972-1790/1031)  26.6 W45-W50 S16-S12
      8 1806  853 (S25W50 L=255)   840     94        0 (1766/ 824-1849/ 889)   5.7 W45-W55 S26-S24
      9 1663  745 (S27W35 L=239)   275    103        0 (1642/ 693-1686/ 784)  45.0 W32-W37 S29-S24
     10 1748  507 (S50W50 L=254)   608     46        0 (1726/ 489-1774/ 528) -24.8 W44-W57 S52-S46
     11 1494  456 (S42W17 L=222)  1304    233        3 (1399/ 378-1574/ 521)  12.5 W07-W25 S45-S38
     12 1564  373 (S53W25 L=230)   545     51        0 (1529/ 358-1594/ 392) -12.5 W20-W29 S55-S51
  

    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