Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  914 1902 (N73W46 L= 76)  1010     95        0 ( 892/1865- 939/1945)   6.5 W32-W85 N68-N75
      2 1740 1498 (N15W68 L= 98)   687     85        0 (1699/1475-1772/1526)  -4.8 W61-W73 N13-N15
      3 1246 1501 (N30W27 L= 58)   311    117        0 (1199/1479-1303/1516) -39.8 W24-W30 N26-N31
      4 1788 1449 (N11W70 L=100)   728     44        0 (1766/1419-1809/1490)  38.7 W66-W74 N09-N13
      5 1729 1421 (N11W60 L= 91)   730     92        2 (1690/1395-1766/1454) -45.0 W56-W64 N08-N14
      6 1121 1415 (N28W16 L= 46)   276     67        1 (1099/1381-1140/1448)  76.0 W14-W16 N25-N29
      7  731 1251 (N27E14 L= 17)   816    331        3 ( 657/1124- 791/1378) -80.5 E17-E11 N18-N36
      8 1060 1305 (N23W09 L= 39)   537    196        0 ( 999/1219-1097/1377)  41.6 W02-W12 N18-N26
      9 1747 1225 (N00W53 L= 83)   905    152        1 (1708/1154-1786/1285) -82.9 W50-W55 S04-N04
     10  601 1070 (N18E25 L=  5)   311    125        0 ( 549/1029- 640/1106) -50.2 E28-E23 N15-N21
     11  229  946 (N16E59 L=331)   388     71        0 ( 199/ 925- 259/ 963)   0.0 E63-E55 N15-N17
     12 1642  881 (S16W36 L= 66)   567    117        2 (1586/ 852-1691/ 924)  -6.3 W31-W40 S17-S14
     13 1526  661 (S27W24 L= 55)   663     91        3 (1499/ 618-1557/ 713) -90.0 W22-W26 S30-S23
     14  634  661 (S07E32 L=359)   259     74        0 ( 616/ 623- 662/ 698)  45.0 E33-E29 S08-S04
     15  592  448 (S19E43 L=347)   352     62        0 ( 559/ 432- 628/ 461) -21.8 E45-E40 S20-S18
     16  874  419 (S27E25 L=  6)   797    171        1 ( 804/ 384- 953/ 454)   4.8 E30-E18 S28-S25
  

    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