Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  970 1686 (N31E25 L=327)   323     76        1 ( 937/1675-1014/1706)  21.8 E27-E21 N29-N31
      2 1499 1622 (N45W14 L=  5)   967    127        2 (1458/1579-1533/1656)  58.0 W10-W16 N40-N49
      3 1566 1595 (N46W20 L= 12)   448     74        2 (1544/1562-1591/1619)  71.6 W18-W23 N42-N48
      4 1724 1385 (N37W37 L= 28)  3396    509        7 (1605/1190-1804/1536) -43.4 W25-W45 N25-N43
      5  416 1482 (N04E52 L=299)   383    119        0 ( 373/1449- 446/1501)   9.5 E55-E48 N02-N05
      6  184 1261 (S12E66 L=285)   933    182        2 ( 141/1220- 217/1291)  11.3 E73-E61 S14-S10
      7 1082 1247 (N07E03 L=348)   211     72        0 (1067/1209-1099/1285) -63.4 E04-W00 N05-N09
      8 1579 1148 (N16W29 L= 20)   439    123        0 (1534/1085-1607/1209) -35.5 W24-W32 N13-N18
      9  300 1110 (S19E49 L=302)   292     65        0 ( 272/1073- 325/1135)  63.4 E51-E47 S21-S17
     10 1497  996 (N05W27 L= 19)   782    168        2 (1438/ 951-1576/1059) -10.3 W22-W33 N03-N07
     11 1881  975 (N18W59 L= 50)   496     62        1 (1860/ 947-1894/1002) -71.6 W56-W61 N15-N19
     12 1742  926 (N09W47 L= 38)  1156    235        1 (1662/ 891-1819/ 945)  19.7 W40-W54 N06-N12
     13 1742  740 (S00W54 L= 45)   921    157        1 (1676/ 681-1792/ 816)  68.2 W50-W56 S05-N05
     14 1156  701 (S21W17 L=  8)   993    334        1 (1033/ 637-1321/ 782)  40.6 W10-W24 S25-S12
     15  767  740 (S30E07 L=344)   205     79        0 ( 730/ 717- 796/ 775)  68.2 E09-E06 S31-S26
     16 1246  556 (S26W28 L= 19)   519    131        0 (1208/ 509-1291/ 614) -26.6 W24-W32 S27-S23
     17 1511  553 (S17W45 L= 36)   295     71        0 (1479/ 527-1550/ 579)  59.0 W43-W46 S19-S14
  

    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