Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  930 1624 (N45W14 L= 74)   735    151        1 ( 870/1589- 993/1647) -21.8 W09-W19 N41-N48
      2 1599 1483 (N15W52 L=112)   938    217        1 (1577/1388-1638/1563)  28.3 W47-W60 N09-N18
      3 1481 1383 (N14W38 L= 98)   555    116        2 (1419/1358-1531/1404) -20.6 W33-W42 N12-N15
      4  800 1204 (N23E08 L= 51)   357     61        1 ( 768/1171- 827/1231) -63.4 E10-E06 N20-N24
      5 1778 1126 (S09W51 L=111)   396     88        0 (1745/1097-1808/1152)   0.0 W47-W55 S10-S08
      6  492  993 (N18E34 L= 26)   530    161        0 ( 444/ 948- 556/1054)  16.7 E37-E27 N15-N19
      7  193 1011 (N24E62 L=358)   492     53        0 ( 159/ 988- 220/1032)   6.3 E67-E58 N23-N25
      8  809  938 (N07E14 L= 46)   244     77        0 ( 788/ 899- 827/ 976)  45.0 E16-E12 N05-N09
      9  267  925 (N18E55 L=  5)   615    168        1 ( 238/ 859- 311/ 963)  71.6 E57-E49 N14-N20
     10  177  854 (N15E68 L=352)   975    175        0 ( 149/ 770- 205/ 924) -78.7 E70-E64 N09-N19
     11 1349  840 (S12W14 L= 73)   295     52        0 (1325/ 808-1371/ 861)  45.0 W11-W15 S13-S10
     12 1625  792 (S23W31 L= 91)   758    184        0 (1566/ 721-1681/ 855)  24.4 W25-W36 S26-S21
     13 1717  767 (S29W39 L= 99)   591    130        0 (1683/ 724-1740/ 840)  45.0 W35-W41 S30-S24
     14 1214  672 (S18W02 L= 61)   959    373        3 (1041/ 652-1392/ 706) -20.9 E08-W13 S22-S13
     15 1724  599 (S41W40 L=100)   573     76        1 (1689/ 577-1752/ 633)  31.0 W35-W43 S42-S38
     16  422  433 (S14E61 L=359)  2108    295        6 ( 383/ 339- 467/ 575)  43.3 E73-E49 S21-S05
     17  705  431 (S20E38 L= 22)   290     56        1 ( 691/ 401- 722/ 462)  68.2 E38-E36 S22-S17
  

    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