Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1425 1643 (N44W35 L=132)  1294    229        4 (1383/1573-1460/1691)  45.0 W31-W41 N39-N48
      2  947 1477 (N32E06 L= 90)   318     80        1 ( 933/1436- 965/1521) -80.5 E07-E04 N29-N35
      3 1732 1466 (N30W57 L=153)   543     89        1 (1701/1441-1759/1487) -31.0 W54-W59 N28-N31
      4 1574 1311 (N21W38 L=135)   490     92        0 (1531/1267-1624/1346) -39.8 W35-W42 N17-N23
      5 1569 1225 (N15W37 L=133)   546    156        2 (1502/1202-1635/1254)   5.2 W31-W42 N14-N17
      6 1120 1166 (N12W06 L=102)   397    108        0 (1091/1114-1150/1228) -66.8 W04-W07 N09-N16
      7  608 1131 (N10E27 L= 69)   213     73        2 ( 574/1106- 642/1152) -21.8 E29-E24 N08-N10
      8 1316 1100 (N08W18 L=115)   263    110        0 (1263/1062-1363/1143) -39.8 W15-W21 N05-N10
      9 1581 1096 (N07W37 L=133)   431    130        0 (1523/1073-1631/1129) -20.6 W32-W40 N05-N09
     10  510  991 (N01E33 L= 63)   418    141        0 ( 457/ 947- 570/1028)  32.0 E37-E29 S02-N03
     11  404  699 (S18E44 L= 53)   376     83        1 ( 388/ 653- 417/ 733)  80.5 E44-E42 S21-S15
     12 1592  620 (S23W41 L=137)   641     98        2 (1556/ 577-1627/ 667)  71.6 W38-W43 S25-S19
     13  725  626 (S22E20 L= 76)   635    115        0 ( 669/ 599- 789/ 661) -20.6 E23-E15 S23-S19
     14  498  589 (S25E38 L= 59)   587    133        2 ( 447/ 571- 559/ 609)   6.3 E42-E32 S26-S23
     15 1137  413 (S37W09 L=105)   412     82        2 (1111/ 377-1160/ 451) -71.6 W06-W11 S40-S34
     16  736  375 (S41E24 L= 73)   432     82        1 ( 708/ 350- 758/ 409)  68.2 E26-E21 S43-S38
     17 1014  344 (S43E01 L= 96)   796    117        2 ( 977/ 310-1064/ 393)  26.6 E03-W03 S46-S38
  

    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