Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1371 1635 (N46W22 L=150)   412     76        1 (1343/1610-1398/1655)  33.7 W19-W25 N43-N47
      2  521 1492 (N25E42 L= 87)  2362    495        8 ( 391/1357- 604/1582) -20.6 E55-E31 N17-N30
      3 1205 1512 (N34W08 L=136)   399     95        1 (1178/1473-1237/1562) -56.3 W05-W10 N31-N37
      4 1647 1339 (N27W42 L=170)   289     90        2 (1614/1319-1684/1354)  -9.5 W39-W45 N25-N27
      5  911 1164 (N09E08 L=120)   419    168        0 ( 875/1099- 933/1252)  90.0 E09-E07 N04-N14
      6 1697 1195 (N18W45 L=174)   967    136        4 (1653/1146-1736/1245) -23.2 W41-W48 N14-N20
      7 1814 1217 (N20W57 L=185)   477     94        2 (1791/1188-1832/1244) -56.3 W54-W59 N18-N21
      8  414 1101 (S00E41 L= 87)  1942    365        2 ( 312/1023- 550/1162)   6.0 E49-E30 S05-N03
      9  544  756 (S20E30 L= 99)   326    162        0 ( 512/ 724- 587/ 802)  38.7 E31-E26 S22-S16
     10  368  676 (S28E45 L= 84)   486    147        0 ( 330/ 649- 400/ 721)  35.5 E48-E41 S29-S24
     11  700  561 (S32E18 L=110)   337     77        3 ( 671/ 540- 731/ 586)  21.8 E20-E15 S33-S29
     12  578  398 (S46E31 L= 97)   378     76        1 ( 559/ 377- 598/ 423)  33.7 E34-E28 S48-S44
  

    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