Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  650 1745 (N57E37 L=276)   569     96        0 ( 603/1723- 705/1763)   4.4 E43-E30 N54-N58
      2  826 1522 (N38E12 L=301)   494    123        1 ( 766/1484- 889/1555)  24.0 E16-E06 N35-N39
      3 1628 1444 (N26W49 L=  2)   392     62        0 (1608/1412-1645/1482) -90.0 W47-W50 N23-N28
      4 1755 1375 (N20W59 L= 12)   962    166        1 (1719/1300-1810/1438) -73.3 W56-W62 N14-N24
      5 1797 1385 (N20W65 L= 18)   639     81        1 (1767/1347-1830/1413) -51.3 W62-W68 N17-N22
      6  231 1227 (N19E60 L=253)  1628    216        2 ( 195/1135- 273/1321)  71.6 E63-E57 N13-N25
      7  932 1132 (N12E05 L=308)   421    187        1 ( 857/1085-1001/1159) -26.6 E09-E01 N08-N13
      8  790 1081 (N09E14 L=299)   232    111        1 ( 775/1022- 809/1125) -90.0 E15-E12 N05-N12
      9  345 1059 (N09E46 L=266)   553    140        1 ( 308/1007- 398/1091)  33.7 E49-E41 N05-N11
     10  610  887 (S02E27 L=286)  1653    442        1 ( 546/ 762- 741/1068) -58.6 E30-E19 S09-N09
     11 1508  954 (S03W30 L=343)   415    135        2 (1454/ 904-1552/1004)  35.5 W26-W33 S05-N00
     12 1407  896 (S05W23 L=336)   227     77        0 (1371/ 872-1437/ 913) -26.6 W20-W25 S07-S04
     13 1606  784 (S14W38 L=351)   306     88        1 (1588/ 737-1635/ 813) -78.7 W36-W40 S17-S12
     14  286  734 (S11E55 L=257)   453     81        1 ( 263/ 702- 322/ 770) -36.9 E58-E51 S13-S08
     15  261  656 (S16E61 L=251)  2257    202        4 ( 205/ 589- 308/ 751) -59.0 E67-E57 S20-S10
     16  860  570 (S23E14 L=299)  1084    339        2 ( 737/ 540-1023/ 619) -11.9 E21-E02 S25-S19
     17 1019  339 (S42E05 L=308)  2453    346        6 ( 917/ 239-1122/ 463)  41.6 E15-W03 S50-S31
     18  597  260 (S48E51 L=262)  1568    129        4 ( 525/ 231- 666/ 287)   7.9 E69-E40 S50-S44
  

    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