Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  868 1747 (N53E18 L=158)   524     69        1 ( 836/1723- 905/1767)  26.6 E21-E14 N50-N54
      2  358 1537 (N33E60 L=116)  2081    208        3 ( 296/1473- 432/1616) -11.9 E88-E49 N29-N37
      3 1346 1552 (N38W24 L=200)   398     99        0 (1308/1536-1391/1571)  -9.5 W20-W27 N36-N39
      4 1424 1511 (N35W30 L=206)   476     71        0 (1395/1481-1454/1539)  18.4 W26-W32 N32-N36
      5  625 1459 (N29E30 L=146)  1396    324        2 ( 482/1425- 765/1492)   0.0 E42-E19 N27-N31
      6 1455 1460 (N31W31 L=207)   287     56        0 (1430/1443-1475/1484)  26.6 W28-W33 N29-N32
      7 1469 1334 (N23W30 L=206)   358     74        0 (1454/1287-1480/1381)  80.5 W29-W31 N19-N25
      8 1793 1161 (N11W56 L=232)   665    140        0 (1735/1113-1832/1198) -29.1 W50-W59 N08-N13
      9  459 1131 (N08E37 L=138)   541    156        0 ( 431/1059- 494/1197) -53.1 E40-E34 N03-N11
     10  807 1045 (N04E13 L=162)   240     65        0 ( 788/1007- 833/1075)  76.0 E14-E11 N01-N05
     11 1847  972 (S00W61 L=237)   541     84        0 (1830/ 934-1863/1014)  63.4 W59-W63 S02-N02
     12 1466  861 (S06W28 L=204)   313     98        0 (1444/ 818-1489/ 917) -66.8 W26-W30 S09-S02
     13  418  780 (S14E41 L=135)   326    114        0 ( 409/ 731- 434/ 838) -90.0 E41-E39 S16-S10
     14 1765  751 (S14W55 L=231)  1102    107        3 (1722/ 703-1811/ 797) -20.0 W50-W61 S16-S10
     15 1080  565 (S26W05 L=180)  1135    227        5 (1009/ 518-1165/ 651) -38.7 W00-W10 S29-S20
     16  878  566 (S26E09 L=167)   387     95        0 ( 828/ 544- 921/ 598)  26.6 E12-E06 S27-S23
     17 1375  216 (S56W45 L=221)   530     48        1 (1361/ 189-1388/ 248) -35.0 W40-W51 S59-S52
  

    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