Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  433 1458 (N31E38 L=191)   785    189        1 ( 371/1421- 508/1496) -21.0 E44-E31 N28-N34
      2 1057 1348 (N18W05 L=234)   320     81        2 (1007/1324-1104/1367)  -9.5 W01-W07 N16-N19
      3 1841 1291 (N08W62 L=292)   810    142        3 (1819/1241-1862/1347) -90.0 W58-W65 N05-N11
      4  307 1305 (N22E48 L=182)   489    144        0 ( 266/1266- 350/1343) -36.9 E51-E43 N19-N25
      5  311 1128 (N11E47 L=183)  1010    211        5 ( 268/1056- 354/1223) -65.6 E50-E43 N06-N17
      6 1275 1146 (N04W16 L=245)   334     93        0 (1260/1095-1291/1197)  80.5 W15-W16 N00-N06
      7 1578 1092 (S02W35 L=265)   599    135        0 (1536/1051-1638/1127) -29.7 W32-W39 S04-N00
      8 1625  974 (S09W38 L=267)   261     81        1 (1599/ 936-1649/1005)  53.1 W35-W39 S11-S07
      9  212  936 (N01E58 L=171)   638    159        1 ( 171/ 893- 247/ 992)  24.0 E63-E53 S01-N04
     10 1415  891 (S13W23 L=252)   447     93        2 (1377/ 850-1441/ 926) -45.0 W20-W24 S14-S10
     11  216  719 (S12E63 L=167)  5377    765       15 ( 111/ 545- 303/ 872) -81.5 E76-E55 S22-S02
     12  497  822 (S09E35 L=194)   220     47        1 ( 482/ 794- 517/ 846)  56.3 E36-E33 S10-S07
     13 1390  721 (S23W21 L=250)   640    197        1 (1321/ 630-1470/ 808)  50.7 W16-W26 S28-S17
     14 1288  729 (S21W14 L=243)   470    154        0 (1218/ 686-1339/ 797)  32.0 W09-W17 S23-S17
     15  509  536 (S26E41 L=188)   830    167        2 ( 463/ 496- 564/ 576)  24.0 E45-E36 S28-S23
     16 1401  389 (S46W25 L=254)  1738    297        4 (1266/ 346-1503/ 453)   8.1 W12-W33 S49-S41
  

    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