Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1065 1679 (N48W18 L=327)   565    113        0 (1032/1641-1115/1701)   0.0 W14-W22 N45-N50
      2 1157 1585 (N39W21 L=330)   828    123        2 (1122/1542-1200/1637) -74.1 W18-W23 N35-N43
      3 1559 1516 (N26W50 L=359)   412     82        0 (1520/1487-1598/1541) -51.3 W47-W51 N23-N28
      4 1377 1375 (N21W29 L=339)   381     99        0 (1318/1359-1424/1389) -23.2 W25-W32 N19-N23
      5  160 1147 (N23E70 L=239)   984    147        1 ( 142/1069- 176/1214) -52.1 E76-E66 N17-N26
      6  477 1108 (N18E35 L=274)   705    229        0 ( 385/1052- 581/1160) -28.1 E42-E27 N13-N22
      7 1474  851 (S11W26 L=335)   689    145        2 (1439/ 787-1517/ 922)  45.0 W22-W29 S14-S07
      8 1692  846 (S16W43 L=352)   628    237        0 (1595/ 788-1760/ 921)  14.9 W34-W49 S18-S12
      9 1786  826 (S20W52 L=  1)  1380    182        4 (1766/ 736-1805/ 919) -72.9 W49-W55 S26-S13
     10  972  788 (S07E07 L=302)   789    185        3 ( 928/ 691-1005/ 865)  55.0 E10-E03 S12-S02
     11 1477  643 (S25W24 L=333)   288    110        0 (1449/ 600-1501/ 699)  56.3 W22-W26 S27-S21
  

    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