Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  646 1796 (N30E52 L= 32)  1497    290        2 ( 615/1730- 672/1870) -47.5 E58-E47 N25-N37
      2 1624 1451 (N36W27 L=111)  2377    433        6 (1498/1298-1747/1582) -24.6 W15-W39 N30-N42
      3  282 1531 (N06E70 L= 14)   778    240        2 ( 259/1492- 309/1570) -21.0 E77-E63 N03-N08
      4  462 1532 (N09E50 L= 33)  1206    209        2 ( 382/1498- 542/1557)  18.4 E56-E44 N06-N10
      5  423 1284 (S07E44 L= 40)  1567    233        6 ( 314/1115- 538/1439)  72.3 E49-E39 S18-N04
      6  931 1415 (N12E16 L= 68)   376    122        0 ( 871/1393- 994/1439)  29.7 E19-E12 N09-N15
      7 1330 1242 (N13W11 L= 95)  1546    345        5 (1184/1074-1422/1378) -28.8 W00-W20 N06-N17
      8  600 1097 (S13E26 L= 57)   238    133        0 ( 547/1030- 651/1140)  63.4 E28-E24 S18-S09
      9  148 1042 (S25E63 L= 21)  1143    313        1 ( 100/1011- 221/1086)  -8.1 E70-E55 S26-S22
     10 1571  434 (S22W56 L=140)   691    149        0 (1518/ 417-1626/ 451)  24.0 W51-W60 S23-S19
  

    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