Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1174 1854 (N65E08 L= 83)  4520    441        6 ( 977/1777-1350/1919)  -1.1 E41-W17 N58-N70
      2 1119 1684 (N45E06 L= 85)  1488    231        3 (1020/1628-1205/1724)  -6.7 E15-W02 N41-N47
      3  516 1665 (N32E55 L= 35)  1038    192        2 ( 464/1589- 588/1709)  78.7 E59-E50 N26-N36
      4  339 1508 (N19E63 L= 28)   574     91        1 ( 321/1472- 373/1548) -15.3 E67-E56 N16-N21
      5  328 1361 (N10E56 L= 35)  1108    210        4 ( 286/1287- 382/1427) -38.7 E62-E49 N05-N13
      6 1492 1206 (N18W28 L=118)  1378    383        3 (1374/1080-1583/1330) -35.2 W18-W35 N11-N24
      7  436 1153 (S01E40 L= 50)   270     70        0 ( 421/1132- 466/1185) -14.0 E41-E37 S02-N00
      8 1853 1073 (N15W62 L=153)  1019    156        0 (1836/1014-1883/1148) -49.4 W59-W66 N12-N19
      9  341  992 (S12E46 L= 44)   346    117        1 ( 303/ 953- 383/1039)  45.0 E49-E42 S14-S08
     10  151  863 (S23E71 L= 19)   756     68        2 ( 133/ 832- 168/ 892) -45.0 E76-E67 S25-S21
     11  257  791 (S26E56 L= 35)   396     88        0 ( 237/ 747- 276/ 832)  45.0 E59-E52 S29-S23
     12  359  771 (S26E44 L= 46)   608    125        1 ( 329/ 740- 406/ 811)  35.5 E46-E39 S28-S23
     13 1191  734 (S15W15 L=105)   258     82        0 (1163/ 703-1226/ 757)  45.0 W13-W16 S17-S13
     14  497  701 (S28E32 L= 59)  1374    316        2 ( 343/ 658- 617/ 737)  17.7 E47-E22 S33-S23
     15 1197  612 (S23W18 L=109)   322    135        1 (1156/ 555-1245/ 665) -32.0 W14-W22 S25-S19
     16 1022  585 (S27W07 L= 98)   346    122        0 ( 978/ 551-1070/ 628)  54.5 W04-W09 S30-S23
     17 1431  510 (S26W38 L=129)   279     42        1 (1414/ 486-1453/ 532) -56.3 W36-W40 S27-S24
     18  670  219 (S69E27 L= 64)   878     87        2 ( 643/ 178- 699/ 262)  58.4 E33-E22 S75-S62
  

    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