Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1269 1717 (N45W39 L=353)  1227    173        1 (1216/1660-1349/1753)   0.0 W31-W47 N41-N46
      2  633 1702 (N60E24 L=289)   519     84        1 ( 605/1658- 655/1740) -58.0 E28-E20 N55-N63
      3 1475 1579 (N31W46 L=359)   414     87        2 (1438/1558-1521/1595) -26.6 W42-W49 N29-N32
      4 1708 1533 (N22W65 L= 19)   434     92        0 (1667/1502-1745/1565) -45.0 W62-W68 N19-N24
      5  994 1455 (N33W06 L=320)  2054    547        2 ( 761/1370-1233/1534) -34.2 E09-W19 N23-N42
      6  257 1405 (N41E63 L=251)   810     81        3 ( 221/1366- 292/1437) -13.2 E72-E55 N37-N42
      7 1597 1341 (N14W44 L=358)   451    156        0 (1533/1282-1635/1373) -63.4 W39-W46 N10-N16
      8  183 1088 (N21E64 L=250)   830     83        1 ( 166/1036- 203/1132)  71.6 E66-E61 N17-N23
      9  274 1071 (N19E53 L=261)   526    142        0 ( 251/1018- 300/1130)  74.1 E55-E50 N15-N22
     10  628 1007 (N12E25 L=289)  1239    373        2 ( 494/ 884- 763/1105) -43.2 E33-E17 N02-N18
     11  892 1060 (N11E07 L=306)   184     72        0 ( 853/1051- 931/1069) -11.3 E09-E04 N10-N11
     12 1588  903 (S11W33 L=347)   451    118        0 (1571/ 840-1603/ 959)  60.3 W31-W35 S14-S07
     13  888  774 (S06E12 L=302)   320    127        0 ( 865/ 721- 908/ 834)  81.9 E13-E11 S09-S02
     14 1067  685 (S14E03 L=311)   483    212        0 (1031/ 636-1122/ 749) -82.9 E04-W00 S17-S09
     15  933  673 (S13E11 L=302)   199     79        0 ( 894/ 658- 976/ 688) -31.0 E13-E08 S14-S11
     16 1243  478 (S32W05 L=318)  1375    328        1 (1129/ 410-1354/ 596)  27.9 E04-W13 S34-S25
     17  949  341 (S36E20 L=294)   288     59        0 ( 915/ 328- 980/ 357)   0.0 E22-E17 S37-S35
     18  748  223 (S43E44 L=270)  1203    154        2 ( 686/ 184- 823/ 266)   8.5 E54-E34 S46-S40
  

    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