Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  914 1902 (N72W46 L= 76)  1005     84        0 ( 891/1865- 939/1944)   8.0 W33-W84 N68-N75
      2 1764 1476 (N13W69 L= 99)  1315    176        5 (1701/1410-1810/1525) -31.0 W61-W74 N08-N15
      3 1245 1505 (N30W27 L= 57)   295     88        0 (1199/1492-1291/1519) -26.6 W24-W30 N28-N31
      4 1125 1421 (N28W16 L= 46)   275     65        1 (1105/1386-1145/1458)  68.2 W14-W17 N25-N30
      5 1734 1424 (N11W61 L= 91)   682     85        2 (1694/1398-1768/1456) -51.3 W56-W64 N09-N14
      6 1067 1313 (N23W09 L= 39)   504    209        0 (1000/1223-1103/1381)  35.0 W03-W13 N19-N26
      7  701 1351 (N34E14 L= 16)   246     84        1 ( 666/1313- 728/1374) -51.3 E16-E12 N30-N35
      8 1751 1227 (N00W53 L= 84)   980    172        1 (1711/1155-1788/1304) -90.0 W51-W55 S04-N05
      9  766 1145 (N20E13 L= 17)   220     52        0 ( 743/1125- 795/1169) -45.0 E14-E11 N18-N21
     10  597 1078 (N19E26 L=  5)   269     77        0 ( 558/1055- 631/1105) -45.0 E27-E23 N16-N21
     11  231  947 (N16E59 L=331)   418     75        0 ( 199/ 926- 262/ 964)   0.0 E63-E55 N15-N17
     12 1648  886 (S16W36 L= 67)   586    133        1 (1591/ 853-1696/ 934)   0.0 W31-W40 S17-S14
     13 1529  665 (S27W25 L= 55)   620     90        2 (1504/ 625-1556/ 716) -90.0 W22-W26 S29-S23
     14  636  659 (S07E32 L=358)   233     75        0 ( 621/ 624- 665/ 687)  36.9 E33-E29 S08-S05
     15  593  449 (S19E43 L=347)   378     58        0 ( 561/ 433- 629/ 463) -21.8 E45-E40 S20-S18
     16  876  420 (S27E25 L=  5)   761    154        1 ( 806/ 385- 948/ 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