Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  978 1519 (N38W10 L=126)  1236    327        1 ( 943/1361-1019/1646)  48.6 W02-W17 N28-N45
      2 1205 1610 (N37W29 L=146)   353     55        1 (1175/1587-1240/1627)   0.0 W26-W32 N36-N38
      3 1433 1588 (N30W44 L=160)   715    137        2 (1379/1560-1506/1608)  -4.8 W39-W51 N27-N31
      4 1636 1474 (N17W54 L=170)   576    168        1 (1568/1442-1703/1504) -31.0 W49-W59 N14-N20
      5 1042 1292 (N22W07 L=123)   650    183        3 ( 996/1235-1096/1366) -78.7 W04-W09 N17-N27
      6  172 1123 (N27E66 L= 51)  1058    126        0 ( 148/1069- 194/1183)  60.3 E69-E62 N23-N30
      7  866  994 (N09E10 L=106)   481    213        2 ( 784/ 962- 976/1018)  -9.5 E15-E03 N06-N09
      8 1782  810 (S24W48 L=164)  2334    350        2 (1701/ 665-1847/ 963)  45.0 W39-W56 S32-S16
      9  249  884 (N12E57 L= 59)   949    251        0 ( 199/ 815- 294/ 957)  37.9 E63-E51 N08-N16
     10 1366  779 (S15W15 L=131)  1088    375        7 (1303/ 647-1414/ 929)  51.3 W08-W20 S21-S06
     11  313  863 (N10E51 L= 65)   400     88        3 ( 294/ 819- 330/ 912)  68.2 E52-E48 N07-N12
     12  689  767 (S01E25 L= 91)   203     88        0 ( 677/ 723- 700/ 811)  53.1 E26-E23 S03-N01
     13 1524  734 (S21W25 L=141)   505    169        2 (1456/ 692-1576/ 767)   0.0 W19-W28 S23-S19
     14  890  638 (S13E16 L=100)   392    187        0 ( 868/ 571- 926/ 728)  84.8 E17-E14 S17-S06
     15 1082  634 (S17E05 L=111)   800    307        3 (1021/ 559-1130/ 691) -84.3 E08-E02 S22-S12
     16 1226  668 (S18W04 L=120)   295     67        1 (1191/ 645-1252/ 689) -45.0 W02-W06 S19-S15
     17  491  450 (S18E51 L= 66)  1242    236        1 ( 388/ 405- 567/ 512)   8.1 E62-E41 S20-S14
     18 1071  335 (S37E15 L=101)  4430    711        7 ( 795/ 216-1320/ 477)   5.1 E41-W06 S42-S33
     19 1680  409 (S53W46 L=162)   579     60        1 (1654/ 397-1710/ 428) -18.4 W40-W55 S56-S50
  

    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