Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1047 1813 (N43E29 L=198)   426    115        0 (1013/1757-1073/1856) -28.3 E36-E23 N39-N46
      2  939 1774 (N36E33 L=195)   572    173        2 ( 922/1695- 968/1837) -47.5 E38-E27 N30-N42
      3  989 1455 (N17E15 L=213)   679    187        0 ( 929/1378-1053/1508) -24.4 E19-E08 N13-N18
      4 1852 1056 (N22W55 L=283)  1041    176        0 (1795/ 974-1882/1109) -18.4 W48-W61 N17-N24
      5  615 1017 (S18E24 L=204)   554    254        0 ( 512/ 968- 735/1055)  39.8 E29-E16 S22-S12
      6  309  908 (S32E46 L=181)   835    271        3 ( 251/ 801- 374/1017)  53.1 E53-E41 S39-S23
      7  715  846 (S25E13 L=215)   330    116        1 ( 660/ 824- 775/ 872)  40.6 E16-E09 S27-S21
      8  793  697 (S31E04 L=224)  1277    316        2 ( 637/ 657- 909/ 728)  18.4 E15-W03 S35-S27
      9 1112  601 (S27W19 L=246)   796    330        0 ( 977/ 535-1247/ 665)   5.4 W08-W29 S28-S25
     10 1438  546 (S20W41 L=269)   490    137        1 (1384/ 504-1495/ 584)  76.0 W39-W43 S23-S15
     11  329  483 (S62E57 L=171)   747     53        1 ( 295/ 467- 361/ 500)   7.8 E69-E47 S63-S59
     12 1502  318 (S29W64 L=292)   880    145        0 (1441/ 279-1561/ 368)  82.9 W62-W66 S32-S24
  

    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