Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  928 1618 (N40W10 L=  8)   273     61        1 ( 897/1601- 965/1638) -21.8 W07-W12 N38-N42
      2  574 1272 (N27E22 L=335)   212     58        0 ( 549/1252- 604/1298) -45.0 E24-E20 N25-N28
      3 1903 1201 (S08W68 L= 65)   567    119        1 (1892/1150-1922/1256)  56.3 W64-W72 S11-S05
      4 1747 1189 (S05W50 L= 47)   637    138        2 (1692/1151-1811/1229)  -4.8 W44-W56 S06-S03
      5 1493 1111 (S03W29 L= 27)  1591    366        2 (1385/ 947-1540/1224) -83.3 W24-W32 S12-N05
      6  551 1165 (N21E26 L=332)   642    217        3 ( 433/1140- 634/1190)  -7.6 E34-E19 N18-N22
      7 1802 1086 (S12W52 L= 50)  1421    227        4 (1752/ 998-1842/1176)  41.6 W45-W57 S15-S07
      8  889 1118 (N11E06 L=352)   431    134        1 ( 832/1099- 964/1143) -29.7 E08-E01 N08-N13
      9 1693 1044 (S11W42 L= 39)   438     87        2 (1670/1015-1727/1088)  36.9 W39-W44 S13-S08
     10 1855  852 (S27W56 L= 53)  1007    141        1 (1798/ 825-1911/ 899) -21.8 W49-W64 S30-S22
     11 1090  789 (S12E01 L=356)   838    184        5 (1014/ 746-1156/ 855)  29.7 E06-W02 S13-S08
     12 1145  648 (S21E02 L=356)   910    242        1 (1049/ 577-1199/ 751)  35.0 E08-W03 S24-S15
     13 1683  648 (S35W37 L= 35)   905    141        2 (1652/ 582-1717/ 698) -77.5 W34-W41 S40-S31
     14  649  645 (S10E31 L=327)   393    131        1 ( 586/ 610- 699/ 686)   7.1 E35-E27 S11-S08
     15 1732  559 (S44W44 L= 42)  3990    249       10 (1666/ 442-1811/ 646)  69.0 W36-W56 S52-S37
     16  768  449 (S25E31 L=327)  2537    578        4 ( 609/ 252- 967/ 600)  20.6 E53-E13 S35-S19
     17  412  422 (S19E60 L=298)   454     43        0 ( 395/ 392- 428/ 450)  53.1 E62-E57 S20-S16
  

    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