Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1333 1848 (N56W39 L=286)  1782    144        1 (1287/1804-1397/1902)  84.8 W34-W45 N50-N62
      2 1484 1717 (N43W44 L=290)   370     36        1 (1463/1696-1511/1731)   9.5 W41-W47 N40-N43
      3 1535 1575 (N31W41 L=287)  3214    485        3 (1328/1503-1735/1629)  -7.9 W23-W59 N26-N35
      4 1837 1423 (N21W68 L=315)   834     51        2 (1807/1393-1859/1445) -31.0 W64-W72 N19-N23
      5 1052 1366 (N18W03 L=249)   788    248        0 ( 955/1302-1162/1430) -33.7 E03-W09 N13-N22
      6 1744 1313 (N14W52 L=298)  1605    319        4 (1670/1207-1826/1423) -52.4 W47-W58 N07-N21
      7 1501 1319 (N14W32 L=278)   651    173        0 (1443/1260-1552/1398) -56.3 W28-W34 N10-N19
      8 1806 1336 (N16W59 L=306)   441     80        1 (1781/1303-1830/1366) -45.0 W57-W62 N13-N17
      9  589 1220 (N10E27 L=220)   249     85        0 ( 549/1202- 630/1240)  21.8 E29-E24 N09-N11
     10  223  782 (S14E61 L=186)  1621    268        2 ( 193/ 655- 252/ 909)  71.6 E65-E57 S21-S06
     11 1553  876 (S12W34 L=280)   219     41        0 (1537/ 852-1568/ 905)  71.6 W32-W34 S13-S10
     12 1733  732 (S21W51 L=298)   626    174        1 (1706/ 658-1770/ 818)  78.7 W50-W53 S25-S15
     13  583  783 (S16E29 L=217)   385    182        0 ( 499/ 758- 664/ 807) -10.3 E35-E23 S17-S14
     14 1793  776 (S18W57 L=303)   372     53        1 (1772/ 747-1816/ 804)  63.4 W54-W58 S20-S16
     15 1355  490 (S37W24 L=270)  1562    365        3 (1230/ 428-1487/ 620)  35.2 W14-W32 S41-S28
     16 1600  461 (S39W48 L=295)   479     75        1 (1584/ 425-1619/ 502) -45.0 W44-W51 S42-S36
     17  366  408 (S39E67 L=180)   845     83        0 ( 331/ 374- 395/ 439)   7.4 E82-E59 S41-S37
  

    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