Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  486 1553 (N48E29 L=140)   467    129        3 ( 460/1501- 507/1595) -63.4 E32-E26 N44-N52
      2  928 1552 (N36W07 L=175)   881    177        1 ( 844/1507-1013/1585)  -8.7 W00-W13 N34-N37
      3  339 1455 (N45E43 L=125)   414    105        3 ( 307/1438- 372/1475)   0.0 E46-E39 N42-N46
      4 1814 1399 (N05W66 L=234)  1277    155        3 (1767/1355-1863/1451) -60.3 W61-W71 N01-N08
      5 1622 1316 (N06W44 L=212)   427     87        1 (1583/1276-1656/1345) -51.3 W41-W45 N02-N07
      6  215 1180 (N29E53 L=115)   582    156        1 ( 168/1154- 253/1200) -12.5 E58-E49 N26-N30
      7 1310  935 (S08W14 L=182)  1968    465        4 (1206/ 732-1370/1110)  82.2 W08-W17 S19-N03
      8  301 1070 (N20E46 L=122)   286    100        1 ( 266/1047- 336/1093) -38.7 E49-E43 N18-N22
      9 1093  926 (S04W02 L=170)  1175    366        4 (1005/ 732-1163/1090)  48.8 E06-W08 S13-N03
     10 1594 1027 (S10W34 L=202)   246     71        1 (1561/1009-1626/1050) -45.0 W32-W36 S11-S07
     11  232  956 (N15E55 L=113)   453    127        3 ( 188/ 929- 258/ 980) -26.6 E59-E52 N13-N16
     12 1412  686 (S25W16 L=184)   797    230        4 (1299/ 600-1512/ 743)   7.1 W07-W23 S28-S23
     13 1195  666 (S22W02 L=170)   192     58        0 (1173/ 635-1225/ 695)  26.6 W00-W04 S23-S20
     14  744  364 (S30E36 L=132)   749    158        2 ( 719/ 305- 767/ 439)  56.3 E39-E32 S34-S25
     15 1395  340 (S51W10 L=178)   449     81        0 (1367/ 310-1418/ 379)  45.0 W06-W11 S53-S47
     16  437  320 (S27E64 L=105)   799    150        0 ( 417/ 280- 459/ 357)  63.4 E65-E61 S30-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