Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  536 1639 (N21E49 L=219)   538    111        2 ( 505/1602- 563/1691) -50.2 E52-E45 N18-N24
      2  530 1554 (N15E45 L=223)   930    211        3 ( 456/1496- 611/1615) -16.4 E53-E36 N12-N18
      3 1511 1462 (N31W22 L=290)  1232    391        0 (1396/1309-1609/1592) -34.4 W11-W30 N23-N37
      4 1348 1522 (N30W09 L=277)   407    150        0 (1323/1456-1378/1591)  90.0 W07-W10 N25-N36
      5  946 1441 (N16E14 L=254)   669    154        2 ( 882/1414-1025/1483)   0.0 E17-E08 N13-N17
      6 1599 1053 (N08W34 L=302)  1919    411        2 (1536/ 878-1678/1211) -63.4 W30-W43 N00-N18
      7 1705 1109 (N14W41 L=309)   259     97        1 (1690/1067-1720/1153) -78.7 W39-W42 N11-N16
      8 1822  941 (N08W55 L=323)  1339    238        1 (1744/ 896-1874/1009)  45.0 W47-W60 N04-N13
      9 1456  760 (S12W32 L=300)   288    113        0 (1412/ 738-1509/ 777)  26.6 W29-W35 S14-S10
     10 1593  594 (S17W48 L=316)   771    121        1 (1534/ 557-1646/ 636)  54.5 W45-W51 S20-S13
     11 1737  596 (S13W61 L=329)   548     77        0 (1713/ 554-1760/ 630) -78.7 W59-W62 S15-S10
     12 1089  562 (S32W16 L=284)   208     68        0 (1046/ 553-1124/ 573)  26.6 W12-W18 S33-S31
     13 1734  517 (S17W67 L=335)  1246    156        2 (1658/ 486-1795/ 560)  36.0 W60-W73 S20-S12
     14 1386  451 (S31W41 L=309)   758    129        0 (1330/ 416-1443/ 494)  66.8 W37-W44 S34-S27
     15 1263  413 (S37W35 L=303)   312     35        0 (1243/ 393-1288/ 435) -11.3 W32-W37 S38-S35
     16 1474  295 (S37W63 L=330)   511     50        1 (1444/ 278-1505/ 314)   6.3 W57-W66 S38-S35
  

    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