Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  758 1690 (N38E25 L= 38)   385     77        0 ( 719/1668- 799/1718) -23.2 E28-E21 N36-N40
      2  374 1622 (N33E59 L=  4)  1401    159        5 ( 344/1548- 407/1682) -42.3 E74-E52 N27-N37
      3  459 1382 (N16E40 L= 23)  1289    426        3 ( 404/1212- 530/1548)  90.0 E43-E37 N05-N26
      4 1498 1448 (N25W30 L= 93)   294    116        0 (1442/1422-1548/1471) -14.0 W26-W34 N23-N26
      5 1149 1416 (N20W06 L= 69)   212    118        0 (1099/1372-1182/1455) -45.0 W02-W08 N17-N22
      6  631 1330 (N13E26 L= 37)   744    293        0 ( 539/1258- 741/1385) -12.1 E32-E18 N08-N16
      7  558 1259 (N08E31 L= 33)   232    103        0 ( 525/1213- 572/1289) -68.2 E33-E29 N05-N10
      8 1933 1096 (N07W71 L=134)   730     55        1 (1914/1071-1954/1118)  -9.5 W68-W76 N05-N08
      9  662  924 (S12E22 L= 41)   196     48        1 ( 636/ 909- 691/ 940)  18.4 E23-E20 S12-S10
     10 1233  879 (S12W14 L= 77)   455    211        2 (1170/ 810-1278/ 936)  58.0 W09-W16 S16-S08
     11  313  847 (S17E49 L= 14)  1461    320        4 ( 206/ 795- 471/ 883)   2.2 E61-E35 S19-S14
     12  553  775 (S21E30 L= 33)   222     37        2 ( 533/ 753- 573/ 797)  71.6 E31-E28 S22-S19
     13 1332  734 (S20W21 L= 85)   683    170        3 (1286/ 675-1373/ 786) -50.2 W18-W24 S23-S16
     14 1809  709 (S17W62 L=125)  1154    156        1 (1749/ 668-1852/ 775)  45.0 W56-W67 S20-S13
     15 1116  718 (S22W07 L= 71)   391    168        0 (1045/ 688-1187/ 750) -16.7 W02-W12 S23-S20
     16  290  693 (S26E55 L=  8)   328     53        1 ( 261/ 677- 317/ 713)  -9.5 E58-E52 S27-S25
     17 1468  500 (S34W37 L=101)   598     72        2 (1429/ 472-1501/ 525) -20.6 W33-W41 S35-S32
     18  948  408 (S45E02 L= 61)   246     67        1 ( 926/ 390- 982/ 436) -26.6 E04-W00 S45-S42
  

    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