Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  901 1704 (N40E14 L= 35)   735    222        2 ( 801/1662- 986/1754) -16.4 E23-E06 N36-N43
      2  456 1682 (N37E54 L=356)  2077    148        4 ( 401/1636- 512/1732) -17.7 E71-E45 N33-N41
      3  514 1520 (N25E39 L= 10)   526    109        4 ( 480/1467- 551/1564) -60.3 E42-E35 N20-N27
      4  647 1391 (N16E26 L= 23)   569    173        0 ( 614/1316- 692/1464)  78.7 E27-E24 N11-N21
      5 1673 1422 (N25W45 L= 94)   651    120        0 (1624/1386-1728/1455) -18.4 W40-W49 N22-N26
      6 1377 1398 (N21W21 L= 70)   248     96        0 (1331/1364-1412/1433) -38.7 W17-W23 N18-N22
      7  190 1382 (N16E69 L=341)   952     76        2 ( 160/1344- 220/1410) -21.0 E78-E64 N13-N18
      8  809 1357 (N15E15 L= 34)   413    140        1 ( 745/1322- 865/1383) -14.0 E19-E11 N12-N16
      9  195 1136 (N00E61 L=349)  2296    380        8 ( 159/ 989- 236/1270)  86.4 E66-E55 S08-N08
     10  622 1169 (N03E25 L= 24)   300     80        0 ( 603/1131- 653/1210) -53.1 E26-E23 N00-N05
     11  312  850 (S17E49 L=  0)   384    144        0 ( 278/ 797- 341/ 915)  49.4 E53-E46 S20-S13
     12  553  860 (S16E30 L= 20)   693    263        1 ( 433/ 823- 664/ 887)  13.2 E38-E21 S18-S14
     13 1556  813 (S14W36 L= 85)   429    121        4 (1492/ 794-1620/ 836)  12.5 W31-W40 S15-S11
     14 1514  739 (S18W34 L= 84)   353     66        0 (1482/ 717-1548/ 766)  36.9 W32-W36 S19-S16
     15 1290  709 (S22W19 L= 69)   346    167        0 (1219/ 682-1362/ 747) -11.3 W14-W24 S23-S19
     16 1517  671 (S22W36 L= 86)   416     91        1 (1480/ 634-1542/ 709)  68.2 W33-W37 S24-S19
     17 1605  482 (S34W51 L=101)   608     89        0 (1570/ 457-1629/ 524) -68.2 W47-W55 S35-S30
     18 1230  419 (S42W21 L= 71)   270     48        1 (1208/ 401-1258/ 447)  14.0 W19-W23 S43-S39
  

    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