Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1193 1677 (N40W35 L=239)   573     82        0 (1151/1662-1248/1695) -38.7 W31-W39 N38-N42
      2 1706 1578 (N17W70 L=274)   745     69        2 (1683/1551-1735/1606)  76.0 W66-W77 N15-N19
      3 1449 1519 (N23W43 L=247)   420    104        2 (1407/1485-1487/1553) -71.6 W40-W44 N20-N26
      4 1323 1400 (N20W29 L=233)   943    348        0 (1177/1321-1461/1462) -49.1 W21-W34 N12-N27
      5 1011 1245 (N20W05 L=209)  2429    540        6 ( 786/1176-1227/1285) -33.1 E08-W16 N10-N25
      6  417 1169 (N29E38 L=167)   406    124        2 ( 377/1119- 455/1222) -58.0 E41-E34 N25-N33
      7 1810  856 (S25W50 L=254)   861    105        1 (1769/ 824-1855/ 892)   5.7 W45-W56 S26-S23
      8 1666  742 (S27W35 L=239)   296    104        0 (1647/ 693-1690/ 787)  38.7 W32-W37 S30-S24
      9 1507  464 (S42W18 L=223)  1412    261        4 (1403/ 385-1596/ 541)  11.9 W08-W27 S45-S38
     10 1753  509 (S50W50 L=254)   577     42        1 (1730/ 494-1777/ 528) -25.0 W44-W59 S53-S46
     11 1567  375 (S53W25 L=230)   516     48        0 (1531/ 361-1597/ 393) -18.4 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