Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  855 1866 (N66W25 L= 41)   750     69        1 ( 836/1824- 875/1899)  28.1 W19-W34 N60-N70
      2 1159 1831 (N49W44 L= 60)  1007    151        1 (1109/1780-1211/1889)   9.9 W35-W58 N46-N52
      3  762 1822 (N65W09 L= 25)   590     65        1 ( 735/1788- 783/1851)  63.4 W06-W15 N60-N68
      4  891 1759 (N54W16 L= 32)   521     99        0 ( 845/1732- 933/1791)   5.2 W11-W22 N52-N55
      5 1066 1709 (N44W26 L= 42)   610    109        1 (1040/1662-1093/1763)  37.9 W21-W30 N40-N47
      6  830 1704 (N51W08 L= 24)   339     65        1 ( 794/1687- 868/1719) -45.0 W05-W10 N48-N53
      7 1420 1655 (N29W47 L= 63)   296     44        1 (1404/1632-1437/1680)  53.1 W45-W48 N27-N31
      8 1340 1593 (N28W37 L= 53)   538     98        0 (1297/1558-1385/1637) -71.6 W35-W39 N25-N31
      9 1423 1413 (N16W34 L= 50)   223     40        0 (1400/1391-1442/1441)  33.7 W32-W36 N14-N17
     10 1472 1350 (N11W35 L= 51)   704    132        1 (1441/1288-1509/1405)  36.9 W31-W39 N07-N13
     11  794 1345 (N28E06 L= 10)   487    174        0 ( 726/1275- 828/1383) -69.4 E09-E04 N23-N31
     12 1144 1345 (N19W15 L= 31)   534    201        0 (1046/1305-1228/1378) -37.9 W09-W19 N14-N23
     13  988 1293 (N20W04 L= 20)   240    111        0 ( 936/1273-1036/1310) -38.7 W01-W06 N17-N21
     14  248 1073 (N24E53 L=323)   592    104        0 ( 221/1019- 275/1114) -54.5 E55-E50 N19-N26
     15  204  932 (N16E60 L=316)  1196    240        3 ( 174/ 830- 236/1011)  84.3 E64-E55 N10-N20
     16 1363  903 (S11W16 L= 32)   295     92        0 (1317/ 887-1417/ 922) -45.0 W13-W19 S13-S08
     17  668  578 (S14E33 L=343)   467     78        0 ( 649/ 530- 684/ 626)  56.3 E35-E30 S16-S10
     18 1150  368 (S39E11 L=  5)  1416    299        0 (1049/ 278-1254/ 465)  20.0 E22-W00 S44-S35
     19 1397  244 (S61W05 L= 21)   530     43        0 (1371/ 227-1425/ 264)   0.0 W01-W09 S63-S58
  

    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