Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  711 1789 (N39E51 L=138)  1383    177        1 ( 647/1727- 778/1846) -10.5 E66-E39 N36-N43
      2  756 1630 (N29E34 L=154)   925    222        2 ( 641/1591- 836/1671)  48.0 E40-E30 N23-N33
      3 1595 1512 (N43W32 L=221)  1576    226        1 (1524/1416-1636/1584) -45.0 W25-W36 N36-N47
      4 1532 1319 (N27W27 L=215)  2474    717        1 (1351/1031-1661/1534) -39.8 W10-W40 N12-N37
      5  897 1510 (N24E20 L=169)   242     75        0 ( 865/1485- 931/1530)  -9.5 E22-E16 N23-N25
      6  432 1481 (N12E52 L=137)   724    116        1 ( 399/1429- 470/1527)  90.0 E53-E49 N08-N15
      7 1106 1424 (N23E04 L=184)   250    106        0 (1068/1390-1152/1464) -15.9 E07-W00 N22-N25
      8 1038 1307 (N15E05 L=183)   216    138        0 ( 993/1273-1097/1337) -15.9 E08-E01 N13-N15
      9 1834 1174 (N26W57 L=246)  1252    120        0 (1809/1115-1854/1233) -74.1 W54-W60 N22-N29
     10  287 1027 (S17E51 L=138)  1847    390        5 ( 185/ 876- 368/1132)  43.5 E63-E44 S28-S10
     11  539  895 (S20E28 L=161)   343    123        0 ( 489/ 856- 570/ 944)  54.5 E31-E26 S23-S16
     12 1300  832 (S07W20 L=209)   238     87        1 (1258/ 814-1339/ 852)  36.9 W18-W22 S08-S05
     13 1516  786 (S05W36 L=224)   276     97        0 (1478/ 757-1567/ 807)  33.7 W33-W39 S07-S02
     14 1469  713 (S10W35 L=224)   273     42        0 (1446/ 688-1490/ 735) -26.6 W32-W36 S11-S08
     15 1397  622 (S17W33 L=222)   674    214        2 (1298/ 603-1485/ 665)  38.7 W27-W39 S20-S12
     16  374  609 (S43E43 L=146)  1189    158        1 ( 297/ 589- 445/ 632)  11.3 E54-E34 S45-S40
     17 1532  504 (S21W49 L=238)   325     90        0 (1489/ 477-1573/ 524)  38.7 W46-W52 S23-S18
  

    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