Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  948 1641 (N48W06 L=326)   417     58        0 ( 919/1618- 982/1659)   9.5 W03-W09 N46-N49
      2 1038 1539 (N38W11 L=331)   548    103        1 (1007/1494-1079/1582) -66.8 W08-W12 N34-N41
      3  937 1441 (N33W01 L=321)   249     87        1 ( 913/1401- 951/1477) -90.0 W00-W02 N30-N36
      4  343 1095 (N19E47 L=273)   798    214        0 ( 264/1037- 431/1155) -26.6 E55-E39 N14-N22
      5 1565  797 (S16W32 L=352)   486    206        0 (1484/ 743-1624/ 874)  22.6 W25-W37 S18-S12
      6  284  842 (N03E55 L=265)   372     52        1 ( 261/ 816- 311/ 871) -76.0 E57-E52 N01-N05
      7  796  766 (S06E18 L=302)   742    184        1 ( 756/ 668- 833/ 835)  56.3 E22-E14 S11-S02
      8 1289  784 (S12W13 L=333)   295     56        1 (1266/ 750-1311/ 814)  53.1 W11-W14 S14-S10
      9 1675  731 (S23W41 L=  1)   734    156        1 (1651/ 679-1720/ 794) -76.0 W39-W46 S27-S19
     10 1298  540 (S29W11 L=331)   431    203        0 (1269/ 453-1338/ 629)  61.4 W08-W14 S34-S23
     11 1012  379 (S36E12 L=308)   265     41        1 ( 994/ 351-1031/ 401)  63.4 E14-E10 S38-S33
  

    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