Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1414 1585 (N26W44 L= 37)   516     99        1 (1387/1538-1442/1630)  74.1 W42-W46 N22-N29
      2 1385 1437 (N19W34 L= 27)   237     87        0 (1354/1397-1418/1457)   8.1 W30-W37 N17-N20
      3 1154 1423 (N25W20 L= 12)   317    111        1 (1100/1402-1209/1449) -50.2 W16-W22 N22-N28
      4 1035 1317 (N23W09 L=  2)   560    203        0 ( 951/1262-1104/1350) -53.1 W04-W11 N17-N26
      5  757 1185 (N23E11 L=341)   526    203        3 ( 721/1116- 799/1261) -84.3 E13-E10 N17-N27
      6  298 1208 (N35E48 L=305)   706    131        2 ( 247/1169- 364/1250) -23.2 E55-E41 N31-N38
      7 1902 1104 (S16W66 L= 59)  1099     95        3 (1859/1075-1949/1131) -22.4 W59-W76 S20-S13
      8 1673 1052 (S10W40 L= 33)   736    186        1 (1628/ 975-1736/1120)  22.6 W35-W47 S13-S08
      9 1866 1004 (S20W58 L= 51)   550    103        3 (1848/ 957-1893/1048)  53.1 W56-W62 S22-S17
     10  121 1001 (N25E75 L=278)  1057    118        3 (  98/ 956- 149/1044) -40.6 E81-E69 N21-N27
     11 1285  786 (S14W08 L=  1)   894    310        1 (1217/ 666-1342/ 923)  38.2 W01-W15 S18-S07
     12  531  847 (N09E34 L=319)   648    176        1 ( 477/ 789- 593/ 914) -60.9 E36-E30 N04-N13
     13 1602  827 (S21W30 L= 23)   353     75        0 (1585/ 792-1617/ 868)  63.4 W28-W31 S22-S18
     14 1044  542 (S21E12 L=341)   316    126        1 ( 999/ 522-1105/ 573) -26.6 E14-E08 S22-S18
     15 1303  537 (S30W03 L=356)   411     38        1 (1283/ 508-1331/ 565)  63.4 W01-W05 S31-S27
     16  861  451 (S22E27 L=326)   348    113        1 ( 831/ 414- 879/ 499)  33.7 E30-E24 S24-S19
     17 1015  252 (S42E29 L=324)   583    107        0 ( 978/ 213-1052/ 306)  28.6 E33-E22 S44-S37
     18  831  229 (S38E44 L=309)   411     66        0 ( 807/ 196- 852/ 260)  66.8 E46-E40 S41-S34
  

    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