Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  863 1726 (N55W15 L=171)   301     51        1 ( 841/1702- 887/1754)  23.2 W11-W18 N53-N57
      2 1264 1549 (N30W32 L=188)   772    153        2 (1203/1507-1339/1581)   0.0 W25-W38 N28-N31
      3  732 1416 (N37E08 L=148)   665    200        2 ( 668/1361- 823/1460)   4.4 E13-W00 N34-N38
      4  601 1313 (N34E20 L=136)   444    114        1 ( 568/1252- 632/1357)  39.8 E23-E17 N31-N36
      5 1730 1185 (S04W49 L=205)   846    183        2 (1690/1100-1766/1265) -84.3 W46-W50 S09-N01
      6 1693  882 (S20W39 L=195)  1290    359        7 (1597/ 774-1773/1038)  32.9 W29-W47 S24-S12
      7 1526  761 (S21W24 L=180)   520    149        1 (1490/ 737-1583/ 819) -54.5 W21-W27 S24-S17
      8 1726  713 (S32W40 L=196)   903    251        2 (1645/ 619-1789/ 800)  20.6 W31-W47 S35-S27
      9 1296  547 (S28W04 L=160)  1069    256        4 (1189/ 509-1407/ 614)  -3.8 E03-W12 S30-S26
     10  998  537 (S20E14 L=142)   328    147        1 ( 935/ 525-1070/ 561) -29.7 E17-E10 S22-S17
     11 1148  407 (S34E10 L=146)   507    143        1 (1113/ 350-1179/ 467)  33.7 E14-E05 S36-S29
  

    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