Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  668 1532 (N34E21 L= 10)  4500    692       10 ( 415/1326- 920/1679) -32.0 E52-E03 N17-N47
      2 1368 1625 (N34W32 L= 63)   538    162        2 (1296/1596-1446/1650) -24.4 W26-W37 N31-N36
      3 1106 1564 (N32W11 L= 42)   275     61        0 (1083/1537-1135/1596) -68.2 W09-W12 N29-N34
      4 1694 1519 (N24W56 L= 88)   412     22        1 (1673/1500-1716/1538)   9.5 W53-W59 N22-N25
      5 1243 1438 (N22W18 L= 49)   284    123        0 (1196/1404-1301/1473) -39.8 W15-W21 N19-N24
      6  415 1319 (N22E39 L=352)   962    217        2 ( 389/1258- 462/1425) -61.4 E42-E35 N17-N28
      7  189 1166 (N14E61 L=330)  3885    437       11 ( 101/ 988- 283/1334)  57.0 E74-E50 N04-N24
      8  317 1204 (N15E47 L=344)   301     52        0 ( 293/1179- 343/1225) -21.8 E49-E44 N13-N16
      9  394 1115 (N09E40 L=351)   968    236        2 ( 304/1078- 500/1164)   7.1 E48-E32 N06-N12
     10  806  764 (S15E16 L= 16)   748    324        1 ( 711/ 660- 898/ 878)  42.7 E23-E09 S20-S08
     11 1485  737 (S22W28 L= 59)  1100    282        3 (1383/ 684-1595/ 834) -32.0 W20-W37 S26-S15
     12  525  742 (S14E35 L=356)  1217    256        1 ( 499/ 616- 571/ 828)  59.7 E37-E30 S21-S08
     13 1681  792 (S20W44 L= 75)   719    126        1 (1630/ 761-1756/ 813)   5.2 W39-W50 S21-S19
     14  448  635 (S20E44 L=348)   271     47        0 ( 417/ 617- 474/ 650) -26.6 E46-E41 S21-S18
     15 1690  478 (S42W56 L= 87)   518     60        0 (1662/ 456-1721/ 504) -18.4 W50-W62 S43-S39
     16  367  457 (S31E62 L=330)   475     44        0 ( 348/ 428- 386/ 482)  23.2 E65-E57 S32-S28
     17  632  436 (S35E36 L=355)   305     61        0 ( 611/ 406- 655/ 469)  33.7 E39-E33 S36-S32
  

    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